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ABSTRACT 


This  project  is  concerned  with  the  relation  of  meteorological  conditions 

to  parameters  and  processes  that  influence  the  optical  propagation  properties 

of  the  turbulent  atmosphere.  The  approach  is  centered  around  the  use  of 

aircraft  and  surface  based  remote  sensors  to  investigate  the  relationship  of 

2 

the  refractive  index  structure  function  parameter,  Cn  ,  to  various 

meteorological  variables  and  to  generally  evaluate  the  accuracy  of 

2 

measurements  of  by  in  situ  and  remote  means.  The  relationship  of  the 
of  the  atmospheric  turbulence  profile  to  the  synoptic  context  and  the  use  of 
physical  models  to  predict  the  profile  using  standard  meteorological  profile 
data  was  also  investigated.  The  study  features  two  modes  of  data  acquisition: 
(1)  continuous  archiving  of  1  hr  average  wind  profiles  and  turbulence  levels 
with  VHF  doppler  radar  profilers,  and  (2)  high  time  resolution  aircraft 
measurements  in  association  with  other  measurements  (groundbased  optical 
scintillometers,  thermosondes,  radar  profilers,  and  doppler  sodar  profilers). 
In  one  case  all  of  the  surface  based  systems  were  assembled  in  an  intensive 
'campaign'  given  the  acronym  EWAK. 

The  atmospheric  turbulence  profiles  and  resultant  optical  propagation 
parameters  have  been  found  to  be  strongly  influenced  by  synoptic  conditions. 

In  particular,  the  turbulence  is  substantially  affected  by  to  strength  and 
location  of  the  jet  stream.  A  very  strong  correlation  between  wind  shear 
(which  is  maximum  above  and  below  the  core  of  the  jet)  and  pilot  reports  of 
turbulence  was  found.  Richardson  number  gave  a  much  weaker  indication, 
possibly  because  of  the  poorer  quality  of  the  vertical  temperature  gradient 
data.  A  comparison  of  five  different  methods  (four  measurement  and  one  model) 
of  obtaining  optical  Cn  showed  average  disagreements  as  large  as  a  factor 
of  three.  The  planetary  boundary  layer  was  found  to  contribute  significantly 
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to  the  correlation  length,  rQ.  A  study  of  the  ratio  of  temperature  to 

velocity  microturbulence  showed  that  the  assumption  of  a  constant  mixing 

efficiency  (used  in  the  Van  Zandt  model)  may  not  be  valid  for  very  weak 

o 

turbulence.  Variations  of  C  at  fixed  altitude  were  found  to  be  as 

n 

large  as  one  order  of  magnitude  for  one  hour  time  or  50  km  space  scales .  An 

evaluation  of  a  Xondar  III  doppler  sodar  was  conducted  to  assess  its  utility 

2 

for  quantitative  monitoring  of  Cn  in  the  boundary  layer.  Profiles  of 
2 

Cn  were  obtained  to  a  maximum  height  of  250  to  900  m  depending  on 
atmospheric  conditions.  We  believe  that  the  boundary  layer  performance  of 
this  system  can  be  improved  significantly  with  minor  modifications  to  the 
hardware  and  software. 
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1 .  Introduction 

1.1  Background 

This  report  describes  work  done  in  a  program  of  atmospheric  measurements 

and  analysis  to  verify  and  improve  optical  measurements  and  models  of 

2 

atmospheric  refractive  index  structure  parameter  (Cn  )  and  other  derived 

atmospheric  optical  parameters  such  as  the  coherence  length  (rQ)  and  the 

isoplanatic  angle  (0Q).  The  primary  in  situ  standard  for  this  work  was  the 

Airborne  Research  Associates'  (ARA)  instrumented  aircraft.  Two  remote  sensors 

were  also  used,  a  VHF  doppler  radar  wind  profiler  and  a  monostatic  doppler 

sodar  wind  profiler.  The  intensity  of  the  backscattered  signal  measured  with 

2  2 

these  remote  sensors  is  proportional  to  radar  Cn  and  C^,  , 
respectively. 

The  basic  goals  of  the  study  were: 

a.  Direct  comparison  of  aircraft  and  stellar  scintillometer  measurements 

2 

of  optical  Cn  vertical  profiles  from  2-10  km  altitude. 

2 

b.  Aircraft  measurements  of  C  below  2  km  altitude  to  assess  the 

n 

relative  importance  of  the  boundary  layer  region  on  rQ  and  0 Q. 

2 

c.  Constant  altitude  aircraft  measurements  of  C^,  to  provide  an 
assessment  of  the  horizontal  variability. 

2 

d.  An  evaluation  of  the  Van  Zandt  et  al.  (1981)  model  for  Cn  . 

e.  A  study  of  the  usefulness  of  third  generation  acoustic  radar  (sodar) 

2 

for  continuous  groundbased  monitoring  of  Cn  profiles  (via 
2 

Cj  )  in  the  boundary  layer. 


The  work  described  here,  which  emphasizes  atmospheric  optical 


-  7  - 


propagation,  was  funded  by  and  done  in  conjunction  with  the  Rome  Air 
Development  Center  (RADC) .  An  adjunct  activity  at  Penn  State  University, 
basic  research  on  the  dynamics  of  free  atmospheric  turbulence,  was  carried  out 
with  funding  from  the  Air  Force  Office  of  Scientific  Research  (AFOSR) .  The 
aircraft  measurements  were  carried  out  in  two  modes:  (1)  A  major  intensive 
field  experiment  involving  a  variety  of  groundbased  optical  measurement 
systems.  (2)  More  limited  studies  involving  aircraft  turbulence  profiles  with 
a  single  groundbased  optical  system  (at  Rome,  NY)  and/or  surface  based  doppler 
systems  (usually  at  Penn  State).  A  total  of  22  flights  were  made. 

The  intensive  field  program  (which  was  given  the  acronym  EWAK)  was  held 
at  Penn  State  in  early  May,  1986.  Given  the  scale  of  the  experiment,  optical 

propagation  scientists  from  several  other  laboratories  were  also  invited  to 

2 

participate.  Four  different  optical  systems  plus  a  balloon  born 

system  (thermosonde)  were  operated  in  conjunction  with  the  aircraft  and 

doppler  systems  (see  Table  1) .  All  optical  equipment  was  operated  between  30 

April  and  6  May.  Fortunately,  during  this  week  skies  were  clear  on  every 

night  but  one.  During  this  period  four  aircraft  flights  were  made  (another 

four  flights  were  made  during  the  optical  setup  period)  and  35  thermosondes 

2 

were  launched.  C  was  also  measured  on  a  short  tower  at  the 
n 

optical/sodar  site  at  Penn  State. 

Following  EWAK,  an  additional  14  aircraft  flights  were  made  over  the  next 
year  and  a  half.  Since  EWAK  took  place  in  the  spring,  the  remaining  flights 
were  made  in  summer,  fall,  and  winter.  Only  the  last  two  flights  were  not 
done  in  conjunction  with  a  groundbased  remote  turbulence  sensing  system.  A 
chronology  is  given  in  Table  2. 


Tc.>'-*»  1.  Summary  of  participants  and  measurements  for  the  EWAK  intensive 
field  experiment  held  at  Penn  State  in  May,  1986. 


Measurement 

Institution 

Contact 

Surface  micrometeorology 

PSU 

C.  Fairall 

Sodar 

PSU 

D .  Thomson 

Radar  (VHF1.VHF2) 

PSU 

D.  Thomson 

2 

Thermosonde  (C  ) 
v  n  ' 

AFGL 

J .  Brown 

Aircraft 

ARA 

R.  Marks on 

2 

Optical  Cn  profile#l 

RADC 

D.  Stebbins 

2 

Optical  CR  profile#2 

AFGL 

E .  Murphy 

Optical  scintillometer  (rQ) 

NPS 

D.  Walter 

Optical  scintillometer  (0Q) 

AFWL 

J .  Davidson 
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Table  2.  Chronology  of  aircraft  flights  for  the  entire  program. 


Flight# 

Site 

Date 

GMT 

Remote  Sensors 

Comments 

1 

PSU 

4/13/86 

1400 

*** 

Test  flight 

2 

PSU 

4/14/86 

1500 

Radar ,  sodar 

3 

PSU 

4/18/86 

2100 

Radar,  sodar 

4 

PSU 

4/19/86 

1700 

Radar ,  sodar 

5 

PSU 

5/3/86 

1400 

Recording  failure 

6 

PSU 

5/4/86 

0400 

EWAK 

7 

PSU 

5/4/86 

0300 

EWAK 

8 

PSU 

5/6/86 

1700 

EWAK 

9 

RADC 

8/11/86 

2250 

2 

Optical  Cn  , sodar 

10 

RADC 

8/12/86 

1500 

2 

Optical  Cn  .sodar 

11 

RADC 

8/13/86 

0200 

2 

Optical  Cn  , sodar 

11-b 

RADC 

8/13/86 

1400 

2 

Optical  Cn  , sodar 

12 

RADC 

10/10/86 

1700 

Optical  Cn2 

13 

PSU 

8/1/87 

2100 

Radar ,  sodar 

14 

PSU 

8/6/87 

2100 

Radar,  sodar 

15 

PSU 

8/11/87 

2100 

Radar ,  sodar 

16 

PSU 

8/12/87 

1700 

Radar,  sodar 

17 

PSU 

9/3/87 

2100 

Radar,  sodar 

Test  flight 

18 

PSU 

9/4/87 

1800 

Radar,  sodar 

19 

PSU 

10/14/87 

2100 

Radar,  sodar 

20 

PSU 

10/15/87 

1800 

Radar,  sodar 

21 

Boston 

12/21/87 

1800 

Horizontal 

22 

Boston 

12/23/87 

1600 

Horizontal 

10 


1 . 2  Scope  of  this  report 

The  analysis  of  the  results  has  been  a  very  large  job  involving  a  number 

of  people.  As  a  result,  certain  components  have  been  completed  and  written  up 

already.  Results  for  the  first  four  goals  outlined  above  have  been  described 

extensively  in  an  MS  thesis  (Beecher,  1987)  and  a  final  technical  report  to 

the  AFOSR  (Fairall  et  al.,  1988).  Both  of  these  documents  are  included  as 

appendices  in  this  report.  The  main  text  of  this  report  will  be  devoted  to  an 

evaluation  of  the  sodar  (item  e  above)  and  a  summary  of  our  overall 

conclusions.  A  discussion  and  evaluation  of  the  accuracy  of  microthermal 
2 

Cn  measurements  is  also  included.  In  the  near  future  we  anticipate 

2 

completion  of  the  in-depth  study  of  the  radar  Cn  applications  in  the  form 

of  a  Ph.D.  thesis  by  Capt.  Michael  Moss  (Florida  State  University).  Most  of 

2 

the  analysis  of  the  radar  Cn  data  presented  in  the  appendix  has  been 
performed  by  Capt.  Moss  as  a  part  of  his  Ph.D.  dissertation  research.  Capt. 
Moss,  a  graduate  student  at  Florida  State  University,  performed  most  of  this 
thesis  research  at  Penn  State  as  it  was  the  only  university  in  the  country  to 
have  VHF  and  UHF  profilers  available. 


2. 


Aircraft  measurement  of  C 

n 


2 


2.1  Theory 
2 

Cn  can  be  related  to  micrometeorological  variables  by  (Wesely,  1976) 


cl  -  A <4  +  2  B  CTq  +  B  C*  ) 


<1) 


where  A  and  B  are  functions  of  temperature  (T) ,  pressure  and  specific  humidity 
2  2 

(q) ;  Cip  ,C  ,  and  are  temperature  and  humidity  structure 

function  parameters.  The  values  of  A  and  B  are  radiation- type  and  wavelength 

2 

dependent.  For  optical  wavelengths,  is  almost  always  dominated  by  the 

temperature  structure  function  parameter  in  the  atmosphere.  Therefore,  we 

2 

will  focus  our  attention  on  Cp  . 

The  temperature  structure  function  parameter  can  be  estimated  from  the 
one -dimensional  temperature  variance  spectral  density,  ^T(f) ,  which  is 
computed  from  a  time  series  of  temperature  fluctuation  measurements . 


C2  -  4  (2*/u)2/3f5/3  *T(f) 


(2) 


where  u  is  the  mean  speed  of  the  air  relative  to  the  sensor  (i.e.,  the 
aircraft  true  air  speed)  and  f  the  frequency.  In  using  (2)  we  are  assuming 
that  the  measurements  are  in  a  valid  region  of  the  inertial  subrange  of 
locally  isotropic  turbulence  (Panofsky  and  Dutton,  1984) . 


2.2  Data  Processins 

The  temperature  fluctuations  necessary  to  use  (2)  are  measured  using 


standard  microthermal  temperature -resistance  sensors  and  very  sensitive 
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resistance  bridges.  Fluctuations  in  the  bridge  output  voltage  is  proportional 
to  atmospheric  temperature  fluctuations .  The  spectrum  of  the  bridge  voltage 
fluctuations,  4> v.  is  expected  to  be 

*v  -  <i  G  Gb  a)2  <*T  +  Nu  +  NT  +  N^iR^)2)  (3) 

where  i  is  the  current  to  the  sensor  with  ice  point  resistance  and 
temperature-resistance  coefficient  a;  Gg  is  the  bridge  gain  and  G  the 
combined  subsequent  signal  conditioning  and  recording  gains;  is  the 
velocity  sensitivity  contamination  noise,  N^,  is  temperature  noise  detected 
by  the  sensors  but  not  due  to  turbulence,  and  Nv  is  the  broadband  voltage 
noise  of  the  detector.  This  relation  and  a  number  of  other  factors  relevant 
to  the  measurement  problem  are  discussed  in  great  detail  by  Fairall  and 
Markson  (1984). 

The  fluctuation  data  are  recorded  in  flight  on  FM  tape  and  later 

processed  in  the  laboratory.  The  time  series  can  be  spectrally  analyzed 

either  with  a  computer  based  digitizer  and  FFT  processor  or  with  a  realtime 

computer  controlled  spectral  analyzer  (Hewlett  Packard  3562A) .  The  results 

given  here  were  all  processed  with  the  HP  device.  Because  of  the 

contaminating  influence  of  the  various  noise  sources  included  in  (3), 

2 

computation  of  Cn  is  not  straightforward.  For  reasonably  strong 

turbulence,  the  real  signal  (^)  overwhelms  the  noise  terms  in  (3)  and  there 

2  -17  -2/3 

are  no  problems.  However,  in  weak  turbulence  (Cn  <10  m  '  )  some 

additional  effort  is  required  to  extract  the  signal  from  of  the  noise. 

Beecher  (1987)  developed  a  useful  algorithm  to  improve  the  automated  analysis 
of  the  data.  Use  of  this  algorithm  is  a  process  of  several  steps: 
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(i)  Noise  subtraction.  Determine  the  broadband  noise  (Nv)  by  looking 

in  the  high  frequency  end  of  the  spectrum  (where  the  real  signal  is 
negligible).  Then  subtract  the  noise  from  the  spectrum. 

(ii)  Defiltering.  Remove  the  effects  of  the  highpass  analog  filter  used 
in  recording.  This  will  restore  the  low  frequency  content  of  the 
atmospheric  spectrum  which  is  most  likely  to  be  well  above  the 
noise. 

(iii)  Slope  checking.  The  slope  of  the  spectrum  is  computed  in  log(^)  vs 

log(f)  coordinates.  According  to  (2),  this  slope  should  be  -5/3  in 
the  inertial  subrange.  The  slope  was  checked  in  several  selected 
frequency  bands. 

Initially  the  slope  of  the  spectrum  was  checked  between  5  and  18  Hz.  If  the 
slope  was  not  within  20%  of  -1.667,  then  various  steps  were  taken  (noise 
subtraction,  defiltering)  and  the  slope  checked  again.  Spectra  that  failed  to 
meet  this  final  criterion  were  simply  rejected.  Further  detail  and  examples 
are  given  in  Chapter  6  of  Beecher  (1987) . 

2 . 3  Calibration 

Examination  of  expressions  (2)  and  (3)  suggests  that  a  number  of  factors 

2 

are  involved  in  computing  Cn  and  that  each  represents  a  potential  source 

of  error.  During  EWAK  two  coincident  aircraft  and  AFGL  thermosonde  flights 

were  made  and  at  some  altitudes  considerable  disagreement  was  observed  in 
2 

Cn  values.  This  disagreement  could  have  been  due  to  sampling  error, 
improper  temperature  bridge  calibration,  processing  errors,  etc.  Lacking  a 
suitable  standard  for  high  altitude  turbulence ,  we  have  checked  the 
consistency  of  the  aircraft  data  by  comparison  with  ground  based  measurements. 

Two  comparisons  have  been  examined.  First,  in  order  to  check  the  bridge 


calibrations  and  general  methodology,  an  AFGL  and  a  Penn  State  temperature 

bridge  were  run  side-by-side  on  the  micromet  tower  at  the  site  at  Rock 

Springs,  PA.  Completely  independent  systems  (sensors  and  bridges)  were 

mounted  on  the  tower  at  a  height  of  8  m.  A  paired  sensor,  rms  processing 

method  was  used  (Fairall  and  Markson,  1984).  The  rms  outputs  of  both  bridges 

were  recorded  as  half-hour  averages  on  the  same  data  acquisition  system. 

2 

About  5  days  data  were  taken.  Values  of  were  computed  and  converted 
2 

to  Cn  using  (1).  The  comparison  of  the  two  systems  is  shown  in  Fig.  1. 

The  average  disagreement  was  9%  with  a  point  to  point  scatter  of  about  25%. 

2 

Direct  comparisons  of  the  aircraft  Cn  measurements  with  tower 

measurements  were  not  feasible  at  State  College  because  of  restrictions  on  low 

altitude  flights  in  a  nearby  populated  area.  Following  EWAK,  the  aircraft  was 

deployed  to  an  experiment  at  Hobbs,  NM,  where  a  ground  station  was  setup  at  an 

2 

abandoned  B-17  airbase.  There  the  ground  station  had  measurements  of  CT 

on  a  tower  at  a  height  of  2  m  above  the  surface  (which  was  very  flat  and 

2 

smooth).  A  number  of  aircraft  profiles  of  C^,  were  made  in  the  area. 

Typically  the  lowest  aircraft  altitude  used  was  6m.  A  sample  aircraft 
2 

profile  of  C^,  and  e  (rate  of  dissipation  of  turbulent  kinetic  energy)  is 

shown  in  Fig.  2.  The  lowest  two  data  points  on  the  graph  (E  for  e  and  C  for 
2 

C^,  )  are  taken  from  the  ground  station.  The  solid  lines  depict  similarity 

model  fits  to  the  data  (Fairall,  1987).  Since  the  aircraft  data  are  not  from 

the  same  altitude  as  the  tower  data,  a  similarity  model  has  been  used  to 

extrapolate  the  aircraft  results  to  2  m.  In  the  example  shown  here,  the  solid 

line  passes  right  through  the  tower  values,  indicating  excellent  agreement. 

2 

Using  this  approach,  we  have  compared  aircraft  and  tower  Cn  for  10  flybys 

2 

of  the  surface  station  (Fig.  3) .  The  values  of  Cn  are  much  larger  at 
Hobbs  than  at  PSU  because  of  the  much  greater  sensible  heat  fluxes  and  the 
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closer  proximity  to  the  surface  (2  m  versus  8  m) .  The  agreement  here  is  quite 
good,  especially  when  we  consider  the  locally  inhomogeneous  surface  (tarmac 
and  grass) . 

2.4  Aircraft  Data 

Some  22  flights  were  made  during  this  project,  three  did  not  yield  usable 

data  (see  Table  2) .  The  data  were  analyzed  in  100  second  blocks  which 

typically  yielded  a  vertical  resolution  of  300  m  for  profiles  and  a  horizontal 

resolution  on  the  order  of  5  km.  Profiles  were  usually  from  0.1  to  10  km 

altitude  but  several  times  the  upper  altitude  was  restricted  due  to  air 

traffic  control.  The  last  two  flights  were  directed  towards  long  horizontal 

runs.  Above  the  boundary  layer,  the  values  were  generally  between 
-18  -16  -2/3 

10  and  10  m'  '  ,  although  higher  values  were  found  in  occasional 

2  2 

strong  CAT  regions.  Besides  the  CR  ,  CT  ,  and  e  data,  aircraft 

measurements  of  mean  potential  temperature  and  mean  water  vapor  mixing  ratio 

are  available.  For  the  profiles  near  State  College,  mean  wind  and  microwave 
2 

Cn  profiles  were  also  available  from  the  Penn  State  radar.  For  flights 
2-12,  all  of  the  available  data  were  summarized  in  the  MS  thesis  by  Beecher. 
Since  this  document  is  included  as  an  appendix,  that  data  has  not  been 
reproduced  in  the  body  of  this  report.  For  flights  13-22,  we  have  included 
the  aircraft  turbulence  and  mean  data  as  Appendix  A.  Radar  data  were  archived 
for  flights  13-20,  but  are  still  being  analyzed  by  Capt.  Moss  and  are  not 


available  at  this  time. 
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Fig.  1  Comparison  of  measurements  of  optical  from  tower  mounted 
sensors  using  a  Penn  State  and  an  AFGL  temperature  bridge. 
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An  example  of  aircraft  measured  profiles  of  C-,  (+  symbol)and  t 

(e  symbol)  from  Hobbs,  NM.  The  solid  curve  is  the  Monin-Obukhov 
similarity  model.  The  ground  stations  measurements  are  designated  by 
E  (for  «)  and  C  (for  C_Z) . 
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Fig.  3  Direct  comparison  of  aircraft  (extrapolated  to  z-2  m)  versus  ground 
station  C  values. 
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3  Sodar  Evaluation 

3 . 1  Sodar  measurements  of  turbulence 

Acoustic  radars  (sodar)  have  been  used  for  quantitative  measurements  of 

2 

C,p  for  some  time  (e.g.,  Neff,  1975;  Thomson  et  al.,  1978).  Doppler 
acoustic  systems  are  now  available  that,  in  principle,  should  allow 

measurement  of  mean  velocity  profiles  and  variability.  As  is  the  case  for 

2 

clear-air  radar,  the  sodar  backscatter  intensity  is  proportional  to  Cn  . 

2 

For  monostatic  sodars  the  operative  variable  is  acoustic  Cn  which,  as  in 
2 

optical  ,  is  dominated  by  the  temperature  structure  function  (Wesely, 

2 

1976).  The  largest  values  of  typically  occur  in  the  atmospheric 

boundary  layer  (ABL)  which,  in  the  daytime,  has  a  vertical  extent  on  the  order 

of  1  to  2  km.  Thus  the  low  level  portion  of  the  Cn  profile  is  often 

critical  in  determining  integrated  optical  path  variables.  In  this  regard, 

the  optical  correlation  length,  rQ,  is  particularly  sensitive.  Since  most 

remote  sensors  with  high  altitude  capabilities  (e.g.,  profiling 

scintillometers  or  clear-air  radar)  are  unable  to  provide  information  about 
2 

Cn  in  the  ABL,  computation  of  rQ  and  $Q  usual)/  requires 

extrapolation  of  the  lowest  level  of  data  down  to  the  surface.  The  final 

results  can  be  quite  sensitive  to  the  extrapolation  method,  so  a  low  level 
2 

Cn  profiling  capability  is  very  desirable.  Historically,  sodar  have  not 

been  able  to  suitably  meet  this  need  because  of  their  limited  vertical  range 

(about  300  m) .  However  there  are  now  sodars  available  that  can  provide 

doppler  wind  measurements  up  to  about  1  km  altitude.  In  this  section  we  will 

examine  the  use  of  one  such  system  (Xontech  Corporation’s  model  Xondar  III)  to 

2 

measure  profiles  of  Cn  . 


3.2  Xhsgry 
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The  basis  for  relating  the  received  power  for  a  sodar  to  the  relevant 
atmospheric  turbulence  parameter  is  the  monostatic  acoustic  radar  equation 
(Neff,  1975): 

pr  -  Pcexp(-2aR)  (cr/2)  [A*G/R2]  o(w)  (4) 

where  pr  is  the  received  power,  pt  the  transmitted  power,  and  R  the  range. 

The  exponential  factor  describes  the  round-trip  loss  of  power  resulting  from 
attenuation  by  air  where  a  is  the  average  absorption  coefficient  for  sound  at 
the  frequency  of  the  sodar.  The  factor  (cr/2)  is  the  radial  extent  of  the 
range  gate  (at  range  R)  which  is  determined  by  the  speed  of  propagation,  c, 
and  the  acoustic  pulse  duration,  r.  The  factor  in  the  square  brackets  ([])  is 
the  solid  angle  subtended  by  the  antenna  aperture  area,  A,  at  range  R  from  the 
scattering  volume,  modified  by  an  effective  aperture  factor,  G,  arising  from 
the  antenna's  directivity.  The  quantity  o(n)  is  the  scattering  cross  section 
per  unit  volume  (for  scattering  180*  from  the  initial  direction  of 
propagation)  and  is  given  by 

a(ic)  -  0.00753  A'1/3  C^/T2  (5) 


where  A  is  the  acoustic  wavelength,  Tv  is  the  virtual  temperature,  and  we 
have  indicated  that  the  sodar  is  sensitive  to  the  virtual  temperature 
structure  function  parameter  (i.e.,  fluctuations  of  atmospheric  density). 
Normally  we  can  assume 


(6) 


to  within  about  20%  accuracy.  Combining  (4)  and  (5)  gives  us  an  expression 
for  the  temperature  structure  function  parameter 
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C2  -  132.8  T l  A1/3R2exp(2aR)/[PTAG(cr/2)Gd]  PR  (7) 

Note  in  (7)  we  have  included  an  additional  factor,  the  total  detector  gain, 

G^,  to  describe  various  stages  of  amplification  to  the  signal  after  it  is 
received  by  the  transducers. 

3.3  The  Xondar  III  system 

The  system  used  for  this  study  was  the  Xondar  III  manufactured  by  Xontech 
Corporation.  This  is  a  monostatic  sodar  using  a  phased  array  antenna  with  25 
transducers.  By  changing  the  phasing  of  the  transducers,  three  different  beam 
orientations  are  generated  to  obtain  doppler  estimates  of  the  three  wind 
components.  All  processing  an  control  is  done  with  an  IBM  PC  computer  with  a 
special  controller  board  and  operating  system.  Data  are  archived  on  5.25" 
floppy  discs.  The  system  specifications  are  summarized  in  Table  3. 

The  system  used  in  the  EWAK  study  was  serial  number  001  of  this  model . 

In  other  words,  Penn  State  acquired  the  first  unit  the  company  sold  and  EWAK 
was  also  our  first  experience  with  it.  Since  1986  we  have  acquired  two 
additional  units.  We  have  operated  these  units  at  State  College,  PA; 
Huntsville,  AL;  Savannah  River  Laboratory,  SC;  San  Nicolas  Island,  CA; 
Monterey,  CA;  and  Harrisburg,  PA  (at  the  Three  Mile  Island  nuclear  power 
facility) .  Unit  001  was  also  operated  briefly  at  Rome  Air  as  a  follow-on  to 
EWAK,  but  because  of  hardware  problems  the  results  were  not  very  useful.  As 
is  typical  with  technically  advanced  instrumentation,  there  were  many  startup 
problems  with  the  systems  and  our  own  learning  curve  in  understanding  their 

Table  3.  Xondar  III  acoustic  radar  system  specifications. 


Number  of  array  elements 


25 


Wavelength 
Transmitted  power: 
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0.21  m 


Electrical  704  W 

Acoustical  176  W 

Operating  frequency  1600  Hz 

Pulse  duration  200  ms 

Pulse  half  length  (cr/2)  34  m 

2 

Antenna  aperture  area  0.5  m 

•fc 

Antenna  effective  aperture  factor  (G  )  0.4 

Total  detector  gain: 

Wind  table  78  dB 

Sounder  output  108  dB 


*  See  Hall  and  Wescott,  1974. 


behavior  and  diagnosing  hardware/software  problems.  Xontech  has  made  several 

improvements  to  the  hardware  based  on  problems  that  we  have  uncovered.  They 

are  also  now  offering  several  different  versions  of  the  system  software  to 

meet  different  needs.  While  this  provides  greater  flexibility,  it  also  has 

added  to  our  confusion  since  it  seems  that  almost  every  experiment  for  which 

we  have  used  the  sodars  has  had  a  different  version  of  the  software. 

The  software  used  during  EWAK  provided  only  the  mean  received  power  as 

part  of  the  wind  table  at  the  end  of  an  averaging  period  (say  10  minutes) . 

2 

Thus  the  information  was  available  with  the  same  temporal  and  spatial 
averaging  as  was  the  wind  information.  A  later  version  of  the  software 
provided  an  additional  output  of  the  received  power  with  one  minute  time 
resolution  and  8.7  m  vertical  resolution.  This  is  the  same  information 


(referred  to  as  'sounder  output')  that  is  used  to  drive  a  facsimile  printer, 
but  is  provided  as  a  serial  stream  that  can  be  archived  digitally  on  a  second 
computer.  This  second  method,  which  was  first  used  on  the  San  Nicolas  Island 
deployment,  was  used  for  the  aircraft  comparisons  in  State  College  in  August 
and  September,  1987.  Archiving  the  received  power  separately  greatly 
complicates  the  analysis  but  it  allows  a  much  closer  look  at  atmospheric  fine 
structure  (both  temporal  and  spatial).  Furthermore,  the  sounder  output  has 
additional  electronic  stages  that  increase  the  total  gain  by  30  dB. 

3.4  Sodar/aircraft  comparison 

2 

Aircraft  Cn  measurements  were  taken  in  the  vicinity  of  the  sodar  at 
Rock  Springs,  PA,  during  two  periods:  EWAK  in  1986  and  August  to  October, 
1987.  The  comparisons  were  made  in  a  series  of  level  runs  by  the  aircraft  at 
fixed  altitude.  Usually,  the  horizontal  runs  were  2  minutes  in  duration 
(about  6  km)  over  agricultural  land  and  parallel  to  Tussey  Ridge.  The  sodar 
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site  was  typically  located  about  the  middle  of  the  run. 

The  EWAK  period  comparisons  are  shown  in  Figs.  4a-4d.  The  solid  curves 
2 

indicate  the  sodar  Cp  computed  from  the  average  power  given  in  the  wind 
table.  The  raw  powers  have  had  a  mean  background  removed.  This  background  is 
believed  to  the  noise  level  (in  some  cases  greater  levels  are  due  to 
environmental  noise  and  wind  noise  of  the  antenna)  of  the  detector.  Removing 
the  background  value  improves  the  signal  to  noise  significantly  at  the  upper 
altitudes.  Since  there  is  a  limit  to  the  accuracy  of  this  noise  subtraction 
process,  a  residual  background  noise  remains,  which  is  crudely  indicated  in 
these  figures  by  the  dashed  straight  line.  The  aircraft  values  fcr  the  2 
minute  level  runs  are  indicated  by  the  circles .  The  comparison  is  quite  good 
although  there  is  a  clear  indication  that  the  sodar  overestimates  Cp^ 
below  100  m.  This  is  probably  due  to  the  acoustic  energy  of  local  ground 
reflections . 

The  1987  comparisons  are  shown  in  Figs.  5a- 5e.  Remember,  these  sodar 
values  are  computed  from  the  sounder  output  which  is  available  about  once  a 
minute.  The  curves  shown  represent  20  minute  time  averages  and  5  point  median 
smoothing  in  the  vertical.  The  noise  level  has  also  been  subtracted  from  this 
sodar  data.  The  residual  noise  level  is  apparent  as  the  smooth  line  above  200 
to  300  m.  For  these  data  the  ambient  Cp  values  are  much  lower  than 
during  the  EWAK  period,  partly  because  most  of  the  flights  were  in  the  late 
afternoon  and  approaching  the  minimum  near  sundown.  Considerably  better 
altitude  coverage  is  observed  at  most  other  times  of  the  day.  Again,  the 

profiles  indicate  that  the  sodar  Cp  are  not  valid  below  about  100  m. 

2 

Only  on  the  afternoon  of  August  12  do  we  see  reasonable  Cp  values  of  the 

2 

300  m.  However,  the  values  for  these  days  appear  to  be  on  the  order 
of  10*^,  which  are  quite  low  for  the  boundary  layer.  Notice  also  that  the 
'noise  curve'  gives  an  indication  of  the  sensitivity  of  the  sodar  as  a 
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EUJAK  UNCORRECTED  SQDAR  PROFILE 
OF 

TEMPERATURE  STRUCTURE  PARAMETER 
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Fig.  4b  Profiles  of  C_2  from  4-18,  16:55  EST  during  EWAK.  The  solid  line 
denotes  e  half  hour  average  of  sodar  data,  the  circles  denote 
2-mlnute  level  aircraft  values. 
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EUAK  UNCORRECTED  9QDAR  PROFILE 
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TEMPERATURE  STRUCTURE  PARAMETER 
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Fig.  4c  Profiles  of  C-2  from  4-19,  12:01  EST  during  EWAK.  The  solid  line 
denotes  a  half  hour  average  of  sodar  data,  the  circles  denote 
2 -minute  level  aircraft  values. 
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Fig.  Ad  Profiles  of  C-  from  5-04,  1:26  EST  during  EVAK.  The  solid  line 
denotes  a  half  hour  average  of  sodar  data,  the  circles  denote 
2 -minute  level  aircraft  values. 
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EM  CTA2  DATA  COMPARISON 


icKcnO  (CB5*2/r(2/3)) 


-  „  e  r  2  fY.nm  o.oi  17*50  EDT  one  year  after  EWAK.  The 

Fi8'  S  Psolideiine  denotes  a  20  minute  median  of  sodar  data,  the  circles 
denote  2 -minute  level  aircraft  values. 
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Fig.  5b  Profiles  of  CT  from  8-06,  18:45  EDT  one  year  after  EWAK.  The 

solid  line  denotes  a  20  minute  median  of  sodar  data,  the  circles 
denote  2 -minute  level  aircraft  values. 


LGflcrc)  (CEn/r(2/3)) 

Fig.  5c  Profiles  of  from  8-11,  18:00  EDT  one  year  after  EWAK.  The 

solid  line  denotes  a  20  minute  median  of  sodar  data,  the  circles 
denote  2 -minute  level  aircraft  values. 
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Fig.  5d  Profiles  of  from  8-12,  12:15  EDT  one  year  after  EWAK.  The 

solid  line  denotes  a  20  minute  median  of  sodar  data,  the  circles 
denote  2 -minute  level  aircraft  values. 
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Fig.  5e  Profiles  of  C^,  from  8-12,  14:28  EDT  one  year  after  EWAK.  The 

solid  line  denotes  a  20  minute  median  of  sodar  data,  the  circles 
denote  2-minute  level  aircraft  values. 
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function  of  altitude.  For  example,  at  z-1  km  a  Cn  of  2*10"  (i.e., 

2  -2 

CT  -2*10  )  would  be  just  detectable  above  the  noise  while  at  200  m  a 

value  of  1*10  ^  is  detectable.  Of  course,  this  height  dependence  of  the 

2 

detectability  is  a  result  of  the  1/r  dependence  of  the  strength  of  the 
returned  signal. 


In  July  of  1987  we  operated  one  of  our  Xondar  sodars  on  San  Nicolas 

Island,  located  about  60  miles  west  of  Los  Angeles.  This  operation  was  part 

of  a  large  field  program  to  study  marine  stratocumulus  clouds.  Since  the 

Eastern  Pacific  is  characterized  by  strong  subsidence  and  weak  convection,  the 

boundary  layer  heights  are  typically  lower  than  those  at  inland  locations; 

between  200  and  1000  m.  This  is  well  within  the  range  of  the  sodar.  The 

sodar  was  operated  for  three  weeks  and  winds  were  obtained  from  100  m  up  to 

the  top  of  the  boundary  layer  continuously  during  the  entire  period.  Thus, 

2 

winds  and  CT  profiles  were  routinely  obtained  to  heights  of  600  and  700 

m. 

Two  sample  profiles  are  shown  (Figs.  6a-6b),  Fig.  6a  shows  a  case  (July 

8  about  midnight)  where  the  inversion  is  at  an  altitude  of  400  m.  Again,  the 

first  100  m  of  the  profiles  are  almost  certainly  not  usable.  The  inversion  is 

2 

clearly  visible  as  a  two  order  of  magnitude  increase  in  CT  in  less  than 

2 

100  m.  Above  the  inversion,  the  C^,  decreases  rapidly  down  to  the  noise 
level.  Note  that  the  noise  curve  is  virtually  identical  to  that  obtained  in 
1987  at  Rock  Springs.  This  probably  indicates  that  the  source  of  the  noise  is 
internal  to  the  sodar.  Much  later  in  the  day  on  July  8  (Fig.  6b)  the 
inversion  had  moved  up  to  650  m.  Notice  that  now  the  noise  level  above  the 
Inversion  is  not  a  nice  smooth  line.  This  implies  that  in  this  case  the  noise 
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level  was  due  to  external  sources  associated  with  the  daytime  activity  at  the 
site  (the  system  was  being  operated  in  a  Western  Sea  Gull  rookery  which, 
combined  with  sea  lions  and  surf  produced  a  continuous  shoreline  cacaphony) . 


HEIGHT  (HUNDREDS  OP"  METERS) 
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2 

Fig.  6a  Profiles  of  Cj  from  8-08,  00:50  GMT  from  San  Nicolas  Island.  The 
solid  line  denotes  a  20  minute  median  of  sodar  data.  The  capping 
inversion  is  visible  as  the  maximum  at  400  m. 


height  (hundreds  OP-  meters) 
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JUtf  8, 1987  (1039) 
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Fig.  6b  Profiles  of  C_^  from  8-08,  10:59  GMT  from  San  Nicolas  Island.  The 
solid  line  denotes  a  20  minute  median  of  sodar  data.  The  capping 
inversion  is  visible  as  the  maximum  at  630  m. 
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4.  Conclusions 

As  we  discussed  in  the  introduction,  the  majority  of  the  work  performed 

on  this  project  has  been  described  in  a  technical  report  (Fairall  et  al., 

1988;  Appendix  B)  and  an  MS  thesis  (Beecher,  1987;  Appendix  C) .  The  main  body 

of  the  present  report  is  intended  to  deal  with  the  subject  of  the  sodar 

performance.  However,  before  moving  on  to  that  subject,  we  will  state  a  few 

conclusions  that  summarize  the  results  of  Appendices  B  and  C. 

2 

*  In  Fig.  4  in  Appendix  B  a  comparison  of  Cn  profiles  from  five 
different  methods  (scintillometer,  thermosonde,  radar,  aircraft,  and  model) 
is  shown  for  two  cases.  The  rms  disagreement  at  any  altitude  is  on  the  order 
of  a  factor  of  three.  This  is  quite  large  and  is  still  basically  unexplained, 
although  it  may  be  primarily  due  to  sampling  variability. 

*  For  the  few  cases  examined  (section  7.4  in  Appendix  C) ,  the  boundary 
layer  contribution  to  6q  was  only  about  1%  but  the  contribution  to  rQ  was 
about  10  to  30  %.  These  were  nighttime  cases.  The  contribution  is  expected 
to  be  higher  during  the  day. 

2 

*  The  variability  of  CR  is  indicated  in  the  difference  between  up 

and  down  aircraft  profiles,  level  aircraft  data  (Appendix  A),  and  the  temporal 
variability  at  fixed  location  (e.g.,  from  the  radar).  Order  of  magnitude 
variations  are  observed  on  one  hour  time  scales  with  a  vertical  resolution  of 
300  to  1000  m.  Similar  variability  is  observed  for  spatial  scales  on  the 
order  of  30  km,  but  surprisingly  it  doesn't  increase  much  for  much  longer 
scales . 

*  The  van  Zandt  (1981)  model  was  found  to  agree  on  average  with  the 

aircraft  data,  but  with  considerable  rms  disagreement.  The  thermosonde  and 

scintillometer  were  about  a  factor  of  2  or  3  lower.  An  analysis  of  the  ratio 
2  2 

of  CT  /Cu  showed  that  the  assumption  of  constant  mixing  efficiency 
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(see  Appendix  B)  may  not  be  valid  for  weak  turbulence.  Thus  the  model  will 

2 

tend  to  underestimate  C  when  it  is  small. 

n 

*  Attempts  to  use  mean  meteorological  profiles  (wind  speed  and 

temperature)  to  compute  Richardson  numbers  in  order  to  predict  layers  of 

severe  CAT  were  not  successful.  A  much  higher  correlation  was  obtained  simply 

by  using  the  wind  shear  alone.  It  is  not  clear  if  this  implies  something 

significant  about  the  dynamics  or  if  it  is  simply  the  result  of  low  quality 

data  or  inadequate  sampling  by  rawinsondes .  In  other  words ,  rawinsonde  data 

do  not  seem  to  be  adequate  to  use  for  computing  Richardson  number. 

2 

Regarding  the  accuracy  of  the  aircraft  CR  measurements,  some 

progress  was  made.  Clearly  the  basic  sensors  and  bridges  used  by  PSU  and  AFGL 

are  in  close  agreement  (Fig.  1)  when  run  side  by  side  near  the  ground. 

Comparison  of  aircraft  measurements  with  a  small  tower  in  New  Mexico  also 

indicated  reasonable  agreement  (Fig.  3).  But,  since  both  of  these  comparisons 

are  for  very  strong  turbulence  (compared  to  the  free  troposphere),  they  cannot 

be  considered  as  definitive  verification  for  high  altitude  performance. 

The  evaluation  of  the  sodar  turned  out  to  be  much  more  complicated  than 

expected.  On  balance,  we  would  have  to  say  that  it  is  inconclusive.  The 

comparisons  with  the  aircraft  look  fairly  good  but  there  are  some 

inconsistencies  that  we  cannot  explain.  Only  recently  did  we  find  out  that 

the  power  given  in  the  wind  table  and  the  power  available  at  the  sounder 

output  do  not  have  the  same  gains.  This  was  established  by  comparing  the  dB 

given  simultaneously  from  each  source.  Thus  it  is  difficult  to  explain  the 

apparent  good  agreement  obtained  during  EWAK  (with  the  wind  table)  and  the 

2 

1987  period  (using  the  sounder  output).  Furthermore,  Cj  values  computed 
from  the  sounder  output  from  the  San  Nicolas  Island  measurements  are  very  high 
(as  large  as  3*10*^)  at  the  inversion.  While  we  don't  have  independent 
measurements  available  to  us  at  this  time,  such  values  are  at  about  an  order 
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of  magnitude  greater  than  those  predicted  by  models.  Finally,  the  behavior  of 
the  noise  background  of  the  sodar  suggests  that  its  sensitivity  has  not  been 
optimally  utilized.  Apparently  some  internal  electronic  noise  (such  as  the 
A/D  converter)  is  limiting  the  performance  rather  than  the  environmental 
acoustic  noise.  Since  the  doppler  system  uses  spectral  processing,  apparently 
much  more  sensitive  intensity  information  is  potentially  available.  This  is 
obvious  since  we  often  obtain  winds  when  the  backscatter  intensity  indicates 
nothing  but  the  background  noise.  In  conclusion  there  is  both  substantial 
room  for  improvement  in  the  hardware  and  much  more  work  will  be  required  to 
establish  methods  for  calibrating  these  systems. 
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APPENDIX  A 


Profiles  of  mean  meteorological  variables  ( 8  and  q) ,  profiles  of  turbulence 
parameters,  and  level  flight  turbulence  parameters  for  flights  13  through  22. 
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Part  AI 

Mean  Meteorological  Profiles 

The  series  of  graphs  in  this  section  depict  the  mean  profiles  of 
potential  temperature  ( 6 )  and  specific  humidity  (q)  on  the  days  and  times 
indicated  on  the  graph. 


ALTITUDE  (km) 
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Part  All 

Turbulence  Profiles 

The  series  of  graphs  in  this  section  depict  the  profiles  of  turbulence 

2 

quantities.  The  upper  panel  is  optical  computed  from  microthermal 

2 

measurements  of  Cj,  assuming  negligible  water  vapor  contribution.  The 

o 

middle  panel  depicts  the  structure  function  parameters  for  velocity  (cu  , 

2 

dotted  line)  and  temperature  (C^,  ,  solid  line) .  The  lower  panel  shows 
2 

C,j,  and  the  rate  of  dissipation  of  turbulent  kinetic  energy,  e  (dotted 
line).  The  flight  number,  date  and  time  are  indicated  on  the  graph. 
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Part  All I 

Horizontal  variability 

The  graphs  in  this  section  depict  turbulence  parameters  measured  during 

level  flight  and  are  intended  to  give  some  indication  of  the  horizontal 

2 

variability  of  the  parameters .  The  upper  panel  is  optical  CR  and  the 

2 

lower  panel  shows  e  (dotted  line)  and  C^,  (solid  line).  The  flight 
number,  date,  GMT,  and  altitude  (in  thousands  of  feet)  are  indicated  on  the 
graph . 
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appendix  b 

Reprint  of  a  final  report  "Long  term  studies  of  the  refractive  index 
structure  parameter  in  the  troposphere  and  stratosphere"  submitted  to  the  Air 
Force  Office  of  Scientific  Research.  Note,  this  report  also  contains 
appendices  A  and  B,  so  beware  of  confusion. 
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ABSTRACT 

This  project  is  concerned  with  the  relation  of  meteorological  conditions 
to  parameters  and  processes  that  influence  the  optical  propagation  properties 
of  the  turbulent  atmosphere .  The  approach  is  centered  around  the 
establishment  of  a  climatology  of  refractive  index  structure  function 
parameter  as  measured  with  a  network  of  Doppler  radars.  The  relation  of  the 
atmospheric  turbulence  profile  to  the  synoptic  context  and  the  use  of  physical 
models  to  predict  the  profile  using  standard  meteorological  profile  data  is 
also  being  investigated.  The  study  features  two  modes  of  data  archiving:  (1) 
continuous  archiving  of  1  hr  average  wind  profiles  and  turbulence  levels,  and 
(2)  high  time  resolution  measurements  in  association  with  other  measurements 
(ground-based  optical  scintillometers,  aircraft  or  radiosondes)  in  an 
intensive  'campaign'  given  the  acronym  EWAK. 

The  atmospheric  turbulence  profiles  and  resultant  optical  propagation 
parameters  have  been  found  to  be  strongly  influenced  by  synoptic  conditions. 

In  particular,  the  turbulence  is  substantially  affected  by  to  strength  and 
location  of  the  jet  stream.  A  very  strong  correlation  between  wind  shear 
(which  is  maximum  above  and  below  the  core  of  the  jet)  and  pilot  reports  of 
turbulence  was  found.  Richardson  number  gave  a  much  weaker  indication, 
possibly  because  of  the  poorer  quality  of  the  vertical  temperature  gradient 
data.  A  comparison  of  five  different  methods  (four  measurement  and  one  model) 
of  obtaining  optical  Cn  showed  average  disagreements  as  large  as  a  factor 
of  nhree.  A  study  of  the  ratio  of  temperature  to  velocity  microturbulence 
shoved  that  the  assumption  of  a  constant  mixing  efficiency  (used  in  the  Van 
Zandt  model)  may  not  be  valid  for  very  weak  turbulence.  The  potential  for 
using  operational  numerical  forecast  models  to  compute  turbulence  estimates 
froa  predicted  wind  and  temperature  profiles  was  examined  in  a  preliminary 


look  at  the  ability  of  the  NMC  Nested  Grid  Model  (NGM)  to  reproduce  the  wind 
speed  and  direction  directly  measured  by  the  radars.  The  standard  deviation 
between  the  radar  and  model  was  on  the  order  of  6  m/s  for  wind  speed  and  15 
degrees  for  wind  direction  at  the  initial  analysis  time.  The  uncertainty  in 
wind  direction  increased  to  about  25  to  30  degrees  for  a  48  hr  forecast  but 
the  uncertainty  in  wind  speed  did  not  change  significantly.  A  systematic 
difference  of  several  m/s  was  found  during  the  winter,  probably  due  to  over 
smoothing  of  mesoscale  features  by  the  model. 
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I.  Introduction 

The  importance  of  turbulence  in  the  atmospheric  boundary  layer,  which  is 
the  basic  living  environment  of  humankind,  has  long  been  recognized.  Only  in 
the  last  few  decades  has  the  importance  of  free  atmospheric  (i.e.,  above  the 
boundary  layer)  turbulence  been  recognized.  Even  in  the  absence  of  clouds, 
stably  stratified  fluids  (e.g.,  the  free  troposphere  and  stratosphere)  are 
observed  to  experience  intermittent  transitions  to  turbulent  regimes. 

Although  it  is  referred  to  as  Clear  Air  Turbulence  (CAT)  it  does  occur  in 
regions  of  stable  clouds  such  alto  or  cirrostratus .  The  beautiful  breaking 
wave  clouds  provide  visual  evidence  of  CAT.  CAT  is  important  in  aircraft 
performance  and  safety,  optical  propagation,  EM  propagation  and 
communications,  atmospheric  dispersion  of  pollutants,  and  as  a  source  of 
viscous  drag  in  the  'free'  atmosphere. 

2 

The  atmospheric  refractive  index  structure  function  parameter,  Cn  , 
is  important  in  a  broad  range  of  optical  propagation  applications  (see  Section 
II).  Examples  of  the  consequence  of  atmospheric  micro turbulence  on  optical 
systems  include:  reduced  intensity  (Yura,  1971),  beam  wander  (Fried,  1966), 

scintillation  and  coherence  (Fried  and  Schmeltzer , 1967) ,  and  anisoplanatism 

2 

(Fried,  1981).  The  consequences  of  spatially  and  temporally  varying  Cn 
can  be  computed  for  a  specific  optical  system,  at  least  In  principle,  by 
specifying  the  values  of  CR  along  the  optical  path.  While  an  extensive 
base  of  data  and  theory  is  available  for  planetary  boundary  layer  (PBL) 
properties  of  (see  Fairall  et  al.,  1982;  Fairall,  1987),  only  a  few 
case  studies  are  available  for  the  free  troposphere  and  the  lower  stratosphere 
(Walters  and  Kunkel,  1981).  The  ALLCAT  (i.e.,  HICAT,  MEDCAT,  etc.)  program  of 
the  late  1960 's  focused  primarily  on  large  scale  turbulence  severe  enough  to 
damage  aircraft  or  discomfit  passengers.  Now  the  structure  and  dynamical 
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2 

properties  of  both  CAT  and  Cn  in  the  free  troposphere  and  lower 

2 

stratosphere  are  receiving  renewed  attention.  For  Cn  this  is  because  of 

the  recognition  of  the  importance  of  the  isoplanatic  angle  in  the  performance 

of  a  variety  of  ground- to- space  systems  that  utilize  adaptive  optics. 

Concurrently  development  of  VHF  clear  air  Doppler  radars,  the  performance  of 

2 

which  are  directly  proportional  to  CR  ,  has  rekindled  the  interest  of  the 

meteorological  community  in  the  climatology  and  dynamics  of  CAT.  It  turns  out 

that  this  is  still  poorly  understood.  For  example,  a  recent  long  term  study 
2 

of  Cn  using  VHF  rcdar  (Nastrom  et  al.,  1981)  revealed  substantial  diurnal 
and  seasonal  variations  even  in  the  stratosphere.  This  study  suggested  that 
strong  tropospheric  convection  and  the  jet  stream  were  relevant  factors  in  the 
intensity  of  CAT  but,  at  present,  this  is  only  conjecture  based  on  physical 
plausibility. 

The  investigation  of  atmospheric  turbulence  involves  the  collection  of 
data,  the  development  of  theory  and  the  implementation  of  models.  The  data 
tells  us  what  is  'up  there'  at  the  time  of  the  measurement,  but  the  theory  and 
models  are  true  indices  of  our  understanding  of  the  phenomena.  With  models  we 
attempt  to  predict  an  important  but  difficult  to  measure  variable  (e.g., 

Cn  )  from  variables  that  we  expect  to  have  at  our  disposal  (e.g.,  regional 
scale  radiosonde  information) .  Dynamical  models  of  turbulence  based  on  the 
Navier-Stokes  and  turbulent  kinetic  energy  (TKE)  equations  have  a  rich 
tradition  in  PBL  research.  Both  second  order  closure  and  large  eddy 
simulations  have  provided  great  insights  into  the  structure  of  PBL  turbulence. 
In  contrast,  the  only  model  of  CAT  in  use  today  is  purely  statistical  in 
nature  (Van  Zandt  et  al.,  1981;  Moss,  1986). 

Remote  sensors  are  ideal  for  the  study  of  climatological  properties  of 
turbulence  in  the  free  troposphere  and  stratosphere.  In  this  regard,  the  VHF 
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and  UHF  Doppler  radars  are  particularly  well  suited  because  of  their  all 
weather  and  day /night  capabilities.  The  Penn  State  Department  of  Meteorology 
is  now  completing  a  mesoscale  triangular  network  of  three  VHP  radars  under 
funding  provided  by  the  DoD  (AFOSR)  University  Research  Instrumentation 
Program.  DoD  (AFOSR  and  ONR)  funding  has  also  been  obtained  for  the 
construction  of  a  multichannel  mm-wave  radiometer  system  for  continuous 
groundbased  measurements  of  temperature  profiles  and  integrated  water 
vapor/liquid.  Subsystems  of  this  radiometer  are  now  being  tested.  It  is 
expected  to  be  online  at  Penn  State  in  mid  1988.  Finally,  a  fourth  Doppler 
radar  has  been  constructed.  It  is  a  UHF  system  similar  in  concept  to  the  VHF 
systems  but  operating  at  0.75  m  rather  than  6  m  wavelength.  The  UHF  system 
has  a  smaller,  more  portable  antenna,  better  vertical  resolution  (100  m  vs. 

300  m  for  the  VHF  systems)  and  far  superior  low  altitude  capabilities  (minimum 
range  of  200  m) .  The  radar  systems  are  capable  of  producing  a  wind  and 
turbulence  profile  with  roughly  30  to  90  second  time  resolution.  Thus,  they 
are  also  ideal  for  highly  detailed  studies  of  turbulence  and  its  relationship 
to  mesoscale  phenomena.  The  combination  of  the  mesoscale  triangle  of  VHF 
profilers,  the  UHF  profiler  and  the  radiometers  constitutes  an  atmospheric 
observing  system  that  represents  a  quantum  jump  in  our  ability  to  study  the 
atmosphere . 

This  report  describes  a  two  year  project  to  use  these  systems  to  study 

2 

the  seasonal  and  diurnal  climatology  of  ,  the  atmospheric  dynamical 

2 

processes  responsible  for  the  variability  of  Cn  ,  and  to  investigate 

2 

measurement  methods  and  models  of  Cn  .  The  results  have  been  published  in 
the  open  literature  (see  section  VIII)  and  in  a  number  of  graduate  assistant 
M.S.  theses  (Beecher  (1987),  Carlson  (1987),  Knowlton  (1987),  Neiman  (1987), 
and  Syrett  (1987)].  The  work  combined  the  full  power  of  the  department's 
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observing  systems,  access  to  national  data  networks  and  weather/satellite 
information,  and  cooperation  with  a  number  of  government  research 
laboratories.  During  the  study  an  intensive  field  program  was  held  at  Penn 
State  involving  in  situ  measurements  with  aircraft  and  thermosonde  balloons 
and  a  variety  of  ground  based  electro-optical  systems  (e.g.,  scintillometers). 
This  study  was  given  the  acronym  EWAK  and  involved  scientists  from  a  number  of 
DoD  laboratories  (see  section  IV).  The  purpose  of  EWAK  was  to  compare 
measurement  methods  and  to  test/develop  models  of  clear  air  turbulence.  The 
emphasis  on  relating  turbulence  to  mesoscale  and  synoptic  scale  structures  was 
intended  to  promote  the  extension  of  the  results  to  different  climatological 
regions.  A  thorough  understanding  o.T  these  relationships  should  lead  to 
assessment  and  forecasting  of  atmospheric  turbulence  with  a  combination  of 
satellite  informational  and  numerical  global  weather  models. 


B-11 


II .  Refractive  Index  of  Air 

For  many  purposes  the  effects  of  atmospheric  gas  on  propagating 
electromagnetic  and  acoustic  radiation  may  be  conveniently  subdivided  into  the 
following  subcategories: 

a)  mean  density  gradients  result  in  beam  refraction  and,  consequently, 
tracking  or  pointing  errors, 

b)  density  and  velocity  fluctuation  caused  by  turbulence  produce 
refractive  index  variations  which,  in  turn,  degrade  system  performance. 

Regarding  the  bulk  radar  refractive  index,  n,  of  atmospheric  gas  it  is 
convenient  to  specify  it  in  terms  of  the  refractivity ,  N  -(n-l)*10^,  in 
whi(.h  case 


N-  77.6  (P  +  4810  e/T)/T  (1) 

where  N,  T,  P  and  e  are  in  units  of  ppm,  K,  mb  (total  pressure),  and  mb  (vapor 

pressure),  respectively.  The  second  term  including  e  specifies  the 

contribution  by  polar  molecules  (principally  water) .  It  Is  often  negligible 

at  optical  frequencies  where  the  humidity  coefficient  is  much  smaller. 

For  specification  of  the  effects  of  turbulence  the  relevant  atmospheric 

2 

properties  are  the  refractive  index  structure  function  parameter,  Cr  , 
(Tatarski,  1961)  and  the  inner  scale  or  Kolmogorov  scale  of  turbulence, 
(Livingston,  1972;  Hill  and  Clifford,  1978).  can  be  related  to 

micrometeorological  variables  by  (Wesely,  1976) 

C2  -  A(C2  +  2  B  CTq  +  B2C2)  (2) 


where  A  and  B  are  functions  of  temperature  (T) ,  pressure  and  specific  humidity 
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2  2 

(q);  ,  Cjq,  and  are  the  temperature  and  humidity  structure 

function  parameters.  The  values  of  A  and  B  are  also  radiation- type  and 
wavelength  dependent.  is  related  to  the  viscosity  and  density  of  air 

(both  functions  of  altitude)  and  the  rate  of  dissipation  of  turbulent  kinetic 
energy,  e  (Hinze,  1975). 


A.  Microturbulence 


In  general,  atmospheric  turbulence  is  anisotropic.  However,  it  is  known 

that  at  small  size  scales  (large  values  of  wavenumber,  k) ,  the  eddies  become 

increasingly  isotropic.  In  the  isotropic  limit  it  can  be  shown  that  (Hinze, 

2 

1975)  the  structure  function  parameter,  Cx  ,  for  the  unspecified 
variable, X,  (X  could  be  u,  T,  q,  or  n) 

c\  -  <(X(r)-X(r+d))2>/d2/3  (3) 

is  independent  of  the  spacing,  d.  In  this  equation  X(r)  denotes  the  value  of 

X  at  position  r  while  X(r+d)  denotes  the  value  at  a  position  a  distance  d  from 

r;  the  brackets  denote  an  average.  Using  dimensional  arguments,  Kolmogorov 

showed  that  in  the  isotropic  limit  the  one -dimensional  variance  spectrum  obeys 
-5/3 

a  k  '  wavenumber  dependence 

*>x< k)  -  Vxe'V3  k*5/3  (4) 

where  xx  is  the  rate  of  dissipation  of  one  half  of  the  variance  of  X  and 
0*  is  an  empirical  constant  (on  the  order  of  0.5  for  velocity  and  0.8  for 
scalars).  It  can  similarly  be  shown  that 

*x(k)  -  0.25  c\  k*5/3  (5) 

where  the  factor  0.25  represents  several  mathematical  constants.  Note  that 
(4)  and  (5)  imply  the  Corrsin  relation 
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Cx  “  4  Vx*  V3  (6) 

Velocity  is  a  special  case  (since  xu“«)  so  that 

C2  -  2.0  «2/3  (7) 


B.  TKE  and  Variance  Budget  Equations 

The  simplified,  horizontally  homogeneous  budget  equations  for  TKE  (symbol 

2 

E)  and  scalar  variance  (symbol  <x  >) ,  are 

DE/Dt  -  -<uw>3U/3z  +  (g /8)d8/dz  -«  +  transport  (8a) 

D(<x2>/2)/Dt  -  -<xw>3X/3z  -v  +  transport  (8b) 

where  the  small  symbols  denote  turbulent  fluctuations  and  the  capital  symbols 
denote  average  values,  U  is  the  streamwise  velocity,  W  the  vertical  velocity, 
g  the  acceleration  of  gravity,  8  the  potential  temperature  and  z  the  altitude. 
Assuming  a  state  of  dynamic  equilibrium,  neglecting  transport,  and  using  an 
eddy  diffusion  coefficient,  K,  to  express  the  covariance  term 

<uv>  -  -K  3U/3z  (9a) 

ul 

<xv>  -  -Kx3X/3z  (9b) 


we  can  use  (8)  to  express  the  dissipations 
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t  -  Km<3U/3zr  [1  -  Ri/Pr] 


(10a) 


Xx  -  Kx(3X/3z)‘ 


(10b) 


where  Pr-Kffl/K^  is  the  turbulent  Prandtl  number  and  Ri  is  the  gradient 


Richardson  number 


Ri  -  N  /S 


The  factor  S  is  the  square  of  the  vertical  gradient  of  the  vector  mean  wind 
and  N  is  the  Brunt-Vaisala  frequency,  which  in  the  atmosphere  is  approximately 
given  by 


"  -  (B/t)  ae/dz 


Using  (10)  we  can  obtain  the  ratio  of  the  scalar  dissipation  to  the  TKE 
dissipation 


*x/£  “  <W  Cx  /Cl  ”  [  (3X/3z)2/N2)  Ri/ (Pr-Ri)  (13) 


In  the  oceanographic  literature  (Gregg,  1987),  this  is  expressed  in  terms  of 


the  mixing  coefficient,  7, 


7X  “  XxN2/((3X/3z)2<] 


In  actively  turbulent  layers  with  Pr-1  and  Ri-0.25,  (13)  implies 
7-Ri/(Pr-Ri)«l/3.  Gossard  and  Frisch  (1987)  have  used  the  shear  budget 
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equations  to  show  that  ^x7/^ua,3/2  while  Gage  and  Nastrom  (1986)  have  shown 
that  equipartition  of  two-dimensional  TKE  and  potential  energy  due  to  gravity 
wave  displacements  in  a  stratified  fluid  implies  for  large 

scale  (anisotropic)  turbulence.  As  we  shall  discuss  later,  relations  derived 
from  (13)  have  been  used  to  compute  e  and  Kffl  from  clear  air  radar  data. 

C.  Clear  Air  Turbulence  Length  Scales 

The  field  of  turbulence  contains  a  bewildering  variety  of  length  scales 
which  we  will  not  attempt  to  discuss  here.  For  our  purposes,  we  focus  on  the 
smallest  scales  present,  the  Kolmogorov  microscale  L^,  and  the  scale  of  the 
dominant  vertical  motions  (eddy  overturning  scale)  Lq,  which  is  related  to 
the  integral  scale  calculated  from  the  vertical  velocity  autocorrelation 
function.  At  the  Kolmogorov  scale,  viscosity  is  rapidly  destroying  the 
turbulent  fluctuations;  the  spectrum  begins  to  deviate  from  the  -5/3  behavior 
at  sizes  an  order  of  magnitude  larger  than  this  (Hill  and  Clifford,  1978).  In 
stratified  turbulence,  the  overturning  scale  is  proportional  to  the  buoyancy 
length  scale,  (Gregg,  1987),  also  referred  to  as  the  Ozmidov  scale.  In 
terms  of  the  turbulent  parameters,  these  scales  are 

1^  -  (v/t3)1/U  (15a) 

-  («/N3)1/2  (15b) 

The  inertial  subrange  occurs  for  those  size  scales  smaller  than  and 
greater  than  L^.  If  then  the  energy  containing  vertical  motions 

are  rapidly  destroyed  by  viscosity;  if  L^-IOL^,  then  no  inertial  subrange 
is  anticipated. 

The  ratio  of  the  length  scales  forms  a  natural  activity  parameter  to 
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classify  the  turbulent  state  (Gibson,  1987) 


A  -  ( W 


4/3  -  €/(vN2) 


Gregg  (1987)  suggests  the  following  empirical  classifications 


A>200 


A>10000 


Decaying  turbulence,  <wu>— <w0>*O 

Isotropic 

Fully  developed 


The  physical  interpretation  of  the  mixing  coefficient  and  the  activity 
parameter  can  be  illustrated  by  noting  that  the  turbulent  diffusivity  can  be 
expressed  as  Kx“At7xu.  Thus,  when  the  product  of  the  activity  parameter 
and  the  mixing  coefficient  exceeds  one,  then  the  turbulent  mixing  processes 
are  more  efficient  than  molecular  diffusivity. 

The  concept  of  length  scales  is  also  used  to  eliminate  the  eddy  diffusion 
coefficients  from  (10)  by  invoking  mixing  length  theory  (Hinze,  1975) 


K  -  CjEJ 


where  c^  is  an  empirical  constant  and  the  square  it  of  E  represents  a 
velocity  scale.  The  dissipation  is  also  related  to  these  quantities 


«  -  c0EJ'  /L 


These  relations  can  be  used  in  (10)  to  yield 
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C2  -  2  £2/3-  2  (Cl/C^/3)  (3U/3z)2  L*/3  (19a) 

C2  -  2  (i9x/^u)  (cl/c2/3)/(Pr*Ri)  <3X/3z)2  L*/3  (19b) 


Notice  that  (19a)  can  be  manipulated  to  give 


,/»3  -  120  -  i» 


1/3 

The  most  common  convention  is  to  fix  the  ratio  c^/Cj  '  “1  and  to  assume 
that  in  actively  turbulent  layers  Ri“0.25.  This  implies  that  LQ-0.35 
(Hocking,  1982).  Because  the  constants  are  chosen  arbitrarily,  there  is  no 
physical  significance  to  this  particular  ratio  of  Lq  to  L^.  Since  is 
considered  to  be  the  outer  limit  of  the  validity  of  the  inertial  subrange,  it 
makes  sense  that  the  true  integral  scale  and  the  energy  containing  scale  are 
larger  than  L^.  It  is  also  clear  that  application  of  these  expressions  is 
likely  to  be  confused  by  our  inability  to  be  sure  of  the  values  of  Pr  and  Ri 
and  by  the  fact  that  the  assumptions  of  stationarity  and  negligible  transport 

will  not  be  valid  in  all  conditions. 

2 

D.  Statistical  Models  of  C 

n 

2 

Van  Zandt  et  al.  (1978,1981)  have  developed  a  model  of  based  on 

a  statistical  integration  of  simplified  forms  of  (19).  In  this  approach,  a 
smooth  mean  profile  with  a  vertical  resolution  roughly  equivalent  to  a 
rawinsonde  Is  the  input.  Velocity  and  temperature  (and,  therefore,  shear  and 
temperature  gradient)  fluctuations  are  allowed  (relative  to  the  mean  profile) 
with  probability  distributions  obtained  from  a  mix  of  empirical  analysis  and 
theory  about  gravity  wave  effects.  Active  turbulent  regions  are  assumed  to 
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exist  when  the  fluctuations  produce  RK0.25.  By  integrating  the  joint 

probability  distribution  over  shear,  temperature  gradient,  and  size  scale 

2 

space,  an  average  value,  <CR  >,  is  computed. 

2 

The  value  of  <C  >  obtained  in  this  manner  is  the  expected  value  at 
n 

some  altitude;  as  such,  it  does  not  actually  contain  any  information  about  the 

vertical  distribution.  We  can  see  that,  at  any  specific  time,  quite  different 
2 

results  for  Cn  can  be  obtained  from  measurements  depending  on  the 
vertical  resolution.  Suppose  that  the  vertical  resolution  (e.g.,  a  radar 
range  gate  of  several  hundred  meters)  is  much  greater  than  the  turbulent  patch 
thickness,  H.  In  this  case,  we  would  expect  that  the  radar  produces 
sufficient  averaging  to  be  consistent  with  model  assumptions .  A  high 
resolution  radar,  or  a  high  vertical  resolution  in  situ  measurement  (e.g.,  an 
aircraft  flying  at  constant  altitude)  is  likely  to  produce  a  measurement  that 
is  either  in  an  inactive  region  or  in  an  active  region.  Thus  a  high  resolution 
measurement  at  a  particular  altitude  may  require  averaging  over  a  long  period 
of  time  to  be  consistent  with  the  model.  The  time  required  would  be  many 
times  longer  than  the  typical  lifetime  of  an  actively  turbulent  layer,  which 
appears  to  be  on  the  order  of  a  few  hours  (Syrett,  1987;  also,  see  the 
discussion  in  Section  III-F). 

In  an  earlier  paper  Van  Zandt  et  al.  (1978)  examined  this  issue  with 
their  statistical  model  by  computing  the  probability  of  turbulence  occurring 
at  a  fixed  altitude,  F.  This  has  been  interpreted  by  some  (e.g.,  Weinstock, 
1981)  as  the  expected  fraction  of  a  range  gate  of  thickness,  2A,  that  is 
occupied  by  turbulence  of  average  scale,  Ljj“L0,  i.e.,  F“L^/(2A). 

Obviously,  this  assumes  that  the  patch  thickness  is  much  less  than  the  range 
gate  thickness  and  that  only  one  layer  is  likely  to  occur  per  range  gate. 

Thus,  average  measurements  produced  in  this  fashion  are  interpreted  as 


<«>  -  F  e 


(21a) 


«£>  -FC^  (21b) 

where  the  values  on  the  right  hand  side  represent  data  within  the  active 

layers  where  Ri«0.25.  While  the  value  of  F  depends  on  the  mean  conditions, 

1/3  2 

Gage  et  al.  (1980)  use  the  model  to  show  that  F  '  N  is  approximately 
constant  and  has  the  value  of  4*10*^  in  the  troposphere  and  8*10’^  in  the 
stratosphere. 

The  Van  Zandt  model  has  never  been  tested  in  detail.  It  has  been 

2 

evaluated  by  comparing  radar  measured  Cn  with  predictions  from 

ravinsondes.  In  a  few  cases,  optical  and  in  situ  data  have  also  been  used. 

2 

On  average,  the  model  does  quite  well  for  ;  this  is  not  shocking 

because  the  model  originally  contained  one  unspecified  constant  which  was 

selected  to  fit  radar  data.  The  internal  details  of  the  model,  such  as  the 

probability  distributions  for  shear,  N,  and  F,  have  not  been  evaluated.  Also, 

the  model  appears  to  overpredict  e  by  two  orders  of  magnitude  (Fairall  and 

2 

Markson,  1984) .  Even  the  model  predictions  of  Cn  have  never  been 
examined  for  a  large  data  set.  Furthermore,  the  model  is  incomplete  in  that 
it  provides  information  about  the  probability  of  observing  turbulence  (F)  and, 
perhaps,  the  typical  size  scale  of  the  turbulence  (<LQ>) ,  but  it  tells  us 
nothing  about  the  vertical  distribution  of  turbulent  regions  nor  does  it 
guarantee  that  the  patch  size  (H)  is  the  same  as  Lq.  Observations  in  the 
ocean  (Gregg,  private  communication)  and  in  the  atmosphere  (Barat  and  Bertin, 
1984)  show  that  H»L0. 

E.  Inactive  Regions 
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The  Van  Zandt  model  partitions  the  atmosphere  vertically  into  regions 

with  active  turbulence  (RK0.25)  and  regions  that  are  considered  to  be 

nonturbulent  (or  the  turbulence  makes  a  negligible  contribution  to  the  average 
2 

of  Cn  ).  By  inactive,  we  mean  that  there  is  no  production  of  TKE  or 
variance  which  implies  that  the  covariance  terms  in  (9)  must  be  zero.  This 
condition  is  met  when  the  activity  parameter  is  less  than  IS.  Thus,  we  can 
define  a  threshhold  value  (Gibson,  1987)  for  the  dissipation  rate,  t ,  such 
that  values  of  <  less  than  imply  decaying  turbulence, 

Ac-  «t/(uN2)  -  15  (22) 

-  8  2-3 

A  typical  value  in  the  free  troposphere  is  «t*5*10  ms  Even  if 
A>15  the  turbulence  may  be  decaying  if  the  destruction  terms  exceed  the 
production  terms.  Sometimes  the  term  'fossil’  turbulence  is  used  to  describe 
this  state  although  that  term  is  also  applied  to  residual  temperature 
structure  that  remains  after  the  velocity  turbulence  has  been  consumed. 

Decaying  turbulence  has  been  extensively  studied  in  the  laboratory  using 
flow  through  grids  (e.g.,  Itsweire  et  al.,  1986).  Atmospheric  studies  are 
quite  rare,  although  Nieuwstadt  and  Brost  (1986)  performed  a  large  eddy 
simulation  model  study  of  the  decay  of  convective  turbulence.  We  can  use  the 
TKE  budget  equation  to  analyze  decaying  turbulence,  where  we  drop  the 
production  terms 

3E/8t  -  -<  (23) 

However,  if  we  eliminate  E  with  (18),  then  we  obtain  a  differential  equation 


for  « 
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d£/e4/3  -  -(3/2)  <c2/L0)2/3  dt  (24) 

If  ve  assume  ChaC  the  size  scale  remains  constant  throughout  the  decay  process 
then  (24)  has  the  solution 

«  -  eo  (1  +  t/t0)*3  (25) 

where  eQ  is  the  value  of  e  at  the  beginning  of  the  decay  process  and  tQ  is 

tQ  -  2  [L2/(c2  £g)]1/3-  c*2/3  N*1  (26) 

Measurements  in  the  ocean  (Dillon,  1982)  suggest  toN“0.3,  which  implies  that 
an  inactive  region  has  a  turbulent  lifetime  that  is  only  on  the  order  of  a 
Brunt-Vaisala  period.  However,  this  is  the  value  of  Brunt-Vaisala  period 
within  the  turbulent  layer  itself,  which  may  be  much  longer  (because  of 
mixing)  than  the  value  associated  with  the  mean  background  temperature  profile 
(about  two  minutes  in  the  atmosphere) . 

F.  Active  Regions 

The  dynamics  and  time  scales  of  regions  of  active  turbulence  are  fairly 
complicated  and  very  little  complete  data  is  available  to  aid  in  the  analysis. 
It  is  generally  conceded  that  the  turbulence  is  initiated  by  Kelvin -Helmholtz 
instability  when  the  Richardson  number  decreases  to  a  value  near  0.25. 

Several  physical  processes  are  able  to  cause  this  decrease.  Nonturbulent 
portions  of  the  atmosphere  are  subjected  to  fluctuations  in  shear  and/or 
potential  temperature  gradient  caused  by  gravity  wave  disturbances.  These 
fluctuations  may  cause  turbulence  as  described  in  the  section  on  statistical 


-  B-23  - 


models  of  CAT.  Synoptic  and  mesoscale  dynamical  processes  also  cause 
evolutions  in  the  shear  and  lapse  rates  which  can  cause  turbulence.  A  classic 
example  of  this  is  the  persistent  regions  of  CAT  found  above  and  below  well 
developed  jet  streams.  If  we  take  the  vertical  derivatives  of  the  standard 
budget  equations  for  the  mean  wind  and  potential  temperature,  then  we  can 
create  budget  equations  for  the  evolution  of  mean  shear  and  mean  lapse  rate. 
These  can,  in  principle,  be  combined  to  yield  a  budget  equation  for  the 
Richardson  number-  a  'Richardsonnumberogenesis '  process,  to  paraphrase 
meteorological  jargon. 

Once  a  layer  of  thickness,  H,  becomes  turbulent,  it  is  of  interest  to 
ponder  the  temporal  duration  of  the  turbulent  event.  If  the  turbulence  is 
caused  by  a  gravity  wave  fluctuation,  then  we  do  not  expect  the  event  to  last 
more  than  a  fraction  of  an  inertial  period  (say,  a  few  hours).  If  the 
breakdown  is  due  to  synoptic  processes,  then  in  principle  the  turbulence  can 
endure  as  long  as  the  'Richardsonnumberogenesis'  can  maintain  the  instability 
against  turbulent  mixing.  Once  a  layer  is  turbulent,  the  mixing  process  tends 
to  cause  Ri  to  increase.  This  is  because  3U/3z  an  38/8 z  are  both  reduced  by 
mixing,  but  3U/3 ?  occurs  as  the  inverse  square  in  Ri.  Thus,  we  anticipate 
that  the  mixing  caused  by  the  turbulence  forces  Ri  to  increase  until  Ri»l.  At 
this  point,  shear  production  of  TKE  is  roughly  canceled  by  buoyant  destruction 
and  the  turbulence  begins  to  decay  due  to  dissipation.  This  process  occurs 
long  before  we  reach  et>  In  other  words,  shear  production  is  still  active, 
the  fluxes  are  still  nonzero,  but  buoyant  destruction  and  dissipation  exceed 
production.  The  layer  is  decaying  but  it  is  still  active. 

This  suggests  that  we  must  consider  two  more  time  scales:  the  time  scale 
for  the  mixing  process  to  occur,  r  ,  and  the  time  scale  for  the  active 
region  to  decay  to  «t<  We  assume  that  the  mixing  process  must  transform 
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kinetic  energy  into  potential  energy  by  destroying  the  ambient  potential 

temperature  gradient  for  a  layer  of  thickness.  Gregg  (1987)  has  shown  that 

2  2 

this  requires  an  amount  to  energy  equal  to  (1/12)N  H  ,  where  N  is  computed 
from  the  background  lapse  rate  when  turbulent  breakdown  occurs.  Suppose  we 
let  P  be  the  shear  production  of  TKE  and  B  the  buoyant  destruction  integrated 
of  the  entire  thickness  of  the  layer  (thus,  transport  terms  become  zero);  then 
in  a  state  of  slowly  evolving  equilibrium,  the  TKE  budget  equation  implies  the 
balance 

P  -  B  +  t  (27) 

If  we  substitute  the  flux  Richardson  number,  Rf-P/B,  we  get  a  relation  for  the 
mixing  time  constant 

% 

«rm  -  (l/Rf-l)Brm  -  (1/Rf -1)N2H2/12  -  [(Pr-Ri)/Ri]  N2H2/12  (28) 

Therefore,  we  can  write  r  as 

m 

rffl-  N2H2/(12  «  7)  (29) 

Aircraft  measurements  of  the  covariance  and  gradient  terms  (Kennedy  and 

Shapiro;  1975,  1980)  in  the  vicinity  of  the  Jet  stream  give  an  average  value 

of  7*0.40.  Representative  values  in  (29)  give  tb  on  the  order  of  1  hour  for 

-4  2-3 

a  100  m  thick  layer  with  €-1*10  m  s  .  Notice  that  (29)  implies  that 
thicker  layers  and  weaker  layers  (e  smaller)  will  persist  much  longer. 


Once  the  layer  has  been  mixed  as  discussed  above,  then  the  turbulence 
begins  to  decay  with  the  time  constant  t  described  in  III-D.  According  to 
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(25),  even  if  tQ  is  relatively  short  (a  few  tens  of  minutes),  it  will  still 
taxe  hours  for  «  to  decay  from  1C'4  to  the  transition  level  of  5*10'8. 
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IV.  Applications  of  VHF  and  UHF  Radar 

The  use  of  radar  data  in  optical  propagation  and  CAT  studies  requires 
consideration  of  several  factors:  radar  resolution,  range  and  measured 
variables,  the  relationship  ot  radar  and  optical  parameters,  and  the  types  of 
studies  that  are  appropriate  for  these  systems.  Details  are  provided  in  the 
literature  (e.g. ,  Hocking,  1982)  and  in  our  previous  proposal,  but  a  brief 
summary  of  the  main  points  will  be  presented  here. 

The  50  MHz  radars  have  vertical  range  gate  resolution  of  300  m  up  to  8  km 
and  900  m  resolution  up  to  18  km.  The  high  resolution  mode  requires  90  s  to 
obtain  one  profile,  the  low  resolution  3  minutes.  In  normal  operations,  a 
pair  of  resolution  profiles  is  obtained  in  5  minutes.  The  UHF  radar  has  100  m 
resolution  to  2.6  km,  300  m  resolution  to  8  km,  and  900  m  resolution  to  12  km. 
Both  systems  have  two  horizontal  and  one  vertical  beam. 

The  basic  raw  data  produced  by  the  radar  is  the  mean  Doppler  shift  (which 
is  used  to  compute  the  mean  wind  vector) ,  the  width  of  the  Doppler  spectrum 
(which  is  related  to  the  turbulence  associated  with  the  wind  variance  or  the 
mean  shear  and  can  be  used  to  estimate  e),  and  the  backscatter  intensity  of 

the  signal  (which  Is  related  to  the  radar  cross  section  and  is  used  to 

2  2 
calculate  radar  Cn  ) .  Details  concerning  the  computation  of  Cn  from 

the  radar  signal  are  discussed  in  Appendix  B.  Thus  the  radar  provides 

quantitative  measurements  of  the  profiles  of  U,  3V/dz,  aw,  and 
2 

Cr  .  Models  can  be  applied  to  this  data  to  estimate  c  (e.g.,  equation 

13),  K  (e.g.,  equation  10),  and  even  80/dz. 

The  radars  normally  operate  continuously  and  archive  data  in  one  hour 

average  blocks.  For  intensive  experiments,  we  archive  the  high  resolution 

(temporally  and  spatially)  data,  including  the  raw  Doppler  spectra  and  all 

2 

moment  tables.  An  example  of  a  time  series  of  Cn  from  two  range  gates  of 
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a  VHF  system  Is  shown  in  Fig.  1.  A  few  examples  of  studies  ideally  suited  to 

2 

these  systems  are:  vertical  distribution  of  Cn  'hot  spots',  comparison  of 

2 

daytime  versus  nighttime  Cn  at  different  times  of  the  year,  coincidence 
of  high  turbulence  levels  and  wind  shear/jet  stream  regions,  and  correlation 
with  synoptic  regimes. 


A .  Background 

This  project  is  concerned  with  the  relationship  of  meteorological 
conditions  to  parameters  and  processes  that  influence  the  optical  propagation 
properties  of  the  turbulent  atmosphere.  The  approach  is  centered  around  the 
establishment  of  a  climatology  of  refractive  index  structure  function 
parameter  as  measured  with  a  network  of  Doppler  radars.  The  relation  of  the 
atmospheric  turbulence  profile  to  the  synoptic  context  and  a  model  to  predict 
the  profile  using  standard  meteorological  profile  data  is  also  being 
investigated.  The  study  features  two  modes  of  data  archiving:  (1)  continuous 
archiving  of  1  hr  average  wind  profiles  and  turbulence  levels,  and  (2)  high 
time  resolution  measurements  in  association  with  other  measurements 
(ground-based  optical  scintillometers,  aircraft  or  radiosondes)  in  an 
intensive  'campaign'  given  the  acronym  EWAK. 

This  project  funded  one  graduate  assistant  (Syrett)  who  studied  the 
correlation  of  turbulence  with  synoptic  context  (in  relation  to  jet  streams). 
Two  Air  Force  graduate  students  have  also  been  working  on  the  project.  Capt. 
Michael  Moss  is  completing  a  Ph.D.  dissertation  involving  radar  data  and  the 
Van  Zandt  model  and  2nd  Lieut.  Elizabeth  Beecher  has  completed  a  comparison 
(M.S.  thesis)  of  meteorological  and  optical  microstructure  from  radar  and  in 
situ  (particularly  aircraft)  measurements. 

B.  The  EVAK  Experiment 

A  major  optical/meteorological  experiment  (acronamed  EWAK)  was  held  at 
Penn  State  from  15  April  to  15  May,  1986.  Rome  Air  Development  Center  (RADC) 
was  planning  a  combined  aircraft/optical  experiment  for  the  spring  and  we 
convinced  them  to  hold  the  field  program  at  PSU  in  order  to  take  advantage  of 
the  radar  data.  Given  the  scale  of  this  experiment,  optical  propagation 
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scientists  from  several  other  laboratories  were  also  invited  to  participate. 

The  enclosed  table  is  a  brief  summary  of  the  measurements  and  participants. 

All  optical  equipment  was  operated  during  an  intensive  period  from  30  April  to 

6  May  when  the  skies  were  clear  on  every  night  but  one.  Four  aircraft  flights 

were  made  (another  four  flights  had  been  made  earlier  to  calibrate  the  radar) 

and  approximately  35  thermosondes  were  launched  during  this  period.  High  time 

resolution  data  were  logged  on  the  VHF  radar  and  the  sodar  was  operated 

continuously  during  the  optical  measurements .  Synoptic  meteorological 

information  was  carefully  analyzed  and  archived.  A  meeting  was  held  at  AFGL 

in  early  September  to  discuss  the  preliminary  results.  One  of  the  Pi's  (CWF) 

and  three  graduate  students  attended  and  made  presentations. 

2 

The  standard  radar  equation  is  used  to  calculate  Cn  from  the 

received  power.  The  calibration  factor  for  the  radar  is  a  combination  of  a 

number  of  system  constants:  antenna  gain,  line  loss,  receiver  noise,  etc. 

However,  the  antenna  gain  for  the  colinear-coaxial  phased  arrays  used  with  the 

Doppler  radars  has  never  been  rigorously  measured.  Thus,  even  the  most 

careful  determination  of  the  other  calibration  factors  cannot  eliminate  all  of 

the  uncertainty.  Because  of  this,  since  the  Penn  State  radars  are  modeled 

after  the  NOAA/WPL  systems,  we  initially  decided  to  use  the  system  constants 

from  the  NOAA  radars.  The  analysis  of  the  EWAK  data  is  still  ongoing,  but  it 

2 

has  revealed  that  the  radar  calibration  is  quite  good;  the  values  of  Cn 

need  to  be  increased  by  about  a  factor  of  two  to  agree  with  the  thermosonde 

balloon  measurements.  A  sample  comparison  of  the  uncorrected  radar, 

2 

thermosonde  and  aircraft  profiles  of  optical  CR  Is  shown  in  Fig.  2. 
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Summary  of 

measurements  for  the  EWAK 

experiment 

Measurement 

Institution 

Surface  micromet 

PSU 

C. Fairall 

Sodar 

PSU* 

D. Thomson 

VHF1.VHF2 

PSU 

D . Thomson 

Thermosonde  (CnA2) 

AFGL 

J .  Brown 

Aircraft 

ARA 

R.Markson 

Optical  CnA2  profile#l 

RADC 

D. Stebbins 

Optical  CnA2  profile#2 

AFGL 

E . Murphy 

Optical  scintillometer(rO) 

NPS 

D. Walters 

Optical  isoplanometer 

AFWL 

J .Davidson 

Analysis  of  the  aircraft  data  has  been  completed  and  published  initially 

in  December  (Beecher,  1987)  in  the  form  of  an  M.S.  thesis.  An  initial 

scientific  publication  has  also  been  completed  (Fairall  et  al.,  1988).  We 

will  present  a  few  examples  from  this  work.  Fig.  3  shows  a  sample  aircraft 

profile  of  microturbulence  data  from  one  of  the  eight  flights  made  during 

o 

EWAK.  Fig.  4  shows  two  comparisons  of  optical  CR  obtained  from  the 
aircraft,  the  thermosonde,  the  radar,  the  Van  Zandt  model  and  optical 
scintillometer  data.  The  optical  data  is  highly  variable  over  the  one  hour 
period  but  is  in  agreement  with  the ' thermo sonde  and  clearly  lower  than  the 
other  two  independent  methods.  Finally,  a  conglomerate  analysis  of  the  small 
scale  turbulence  from  all  flights  is  shown  in  Fig.  5.  Here  the  relation  of 
thermal  and  velocity  microstructure  (as  in  equation  13)  is  examined.  Notice 
that  for  the  special  case  of  temperature 


(30) 


(8X/dz)2/ N2  -  (T/g)  86/8 z 

2  2 

By  plotting  Cj,  versus  (T/g)  86 /dz  Cu  in  log- log  format,  we  expect 
the  points  to  define  the  dependence  of  If  1$  is  constant, 

then  we  expect  a  straight  line  with  a  slope  of  one.  Clearly  this  is  not  what 
we  observed.  For  strong  turbulence  (clearly  active  and  nondecaying  regions 
with  <>10'^),  the  data  gives  a  value  of  7^  of  about  0.5,  which  is 
consistent  with  the  jet  stream  value  0.4  of  Kennedy  and  Shapiro  (1980)  and  the 
oceanic  value  0.3  of  Gregg  (1987).  For  the  weaker  turbulence  cases,  this 
value  increases  by  about  a  factor  of  10.  The  significance  of  this  is  not 
clear.  Presently  we  are  not  sure  of  the  effects  of  the  vertical  averaging 
necessary  to  produce  the  300  m  resolution  which  is  required  in  order  for  the 
data  to  be  compatible  with  that  from  the  radar. 

The  second  part  of  the  EWAK  study  is  an  in  depth  investigation  of  the 
o 

radar  Cn  and  wind  data  for  the  EWAK  period  by  Mike  Moss,  an  Air  Force 
graduate  student  enrolled  at  Florida  State  University.  This  work  involves  the 
Van  Zandt  model,  the  thermosondes ,  and  the  optical  data.  Since  this  is  a 
Ph.D.  thesis  project  (which  requires  three  or  four  years  of  work),  we  do  not 
anticipate  final  results  until  next  year.  A  preliminary  description  of  the 
work  in  progress  has  been  prepared  (Moss,  1986). 

C.  Climatological/Synoptic  Studies 

The  initial  steps  in  the  climatological  study  of  the  synoptic  context  of 

2 

Cn  have  been  completed.  We  have  been  routinely  archiving  winds  and  power 
at  one  hour  intervals.  During  a  one  year  period,  each  hour  was  assigned  a 
synoptic  classification  based  on  the  weather  charts.  A  grand  ensemble 
analysis  of  this  data  base  has  not  been  started.  For  the  purposes  of  this 
project,  we  have  completed  a  detailed  analysis  of  one  selected  synoptic 
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feature  that  is  of  particular  interest:  the  jet  stream  and  the  wind  shear  in 

its  vicinity.  This  work  has  now  been  published  initially  in  the  form  of  an 

M.S.  thesis  (Syrett,  1987).  A  scientific  publication  is  now  in  preparation 

(Syrett  and  Thomson,  1988).  Another  phase  of  this  work  involves  comparison  of 

the  radar  measured  winds  with  both  analyzed  and  predicted  winds  for  the  NWS 

Nested  Grid  Model  (NGM) .  A  paper  on  this  topic  has  been  presented  (Thomson  et 

2 

al.,  1988).  This  work  is  relevant  to  the  numerical  prediction  of  Cn  from 
NGM  products. 

Hourly  measurements  of  wind  speed  and  direction  were  examined  for  two 
prolonged  jet  stream  occurrences  over  western  Pennsylvania.  Data  from  two  of 
the  Penn  State  radars  (McElevys  Fort  and  Crown)  were  stratified  into 
categories  based  on  location  of  the  jet  axis  relative  to  the  site.  Low 
resolution  data  from  the  Crown  radar  were  also  compared  to  the  Pittsburgh 
rawinsonde.  Potential  temperature  profiles  were  obtained  using  isentropically 
interpolated  T  and  Td  soundings.  The  combination  of  measured  wind  and 
interpolated  temperature  profiles  allowed  low  resolution  Ri  profiles  to  be 
generated  for  the  profiler  sounding  volume.  Both  Ri  and  wind  shear  statistics 
were  examined  along  with  pilot  reports  of  turbulence  in  the  vicinity  of  the 
profiler. 

Two  cases  were  examined.  Case  1  lasted  from  7  November  to  14  November 
1986.  Jet  stream  case  2  lasted  from  January  IS  to  January  23,  1987.  A  sample 
time-height  cross  section  of  wind  speed  (Fig.  6)  illustrates  the  structure  and 
variability  of  the  jet  stream  at  the  Crown  site  (Just  north  of  Cook  Forest  in 
western  Pennsylvania).  The  data  were  separated  into  five  position  categories. 
The  average  wind  speed  and  Ri  profiles  for  case  2  are  shown  in  Fig.  7.  The 
vertical  resolution  of  the  profiler  data  is  slightly  less  than  1  km,  but  was 
interpolated  to  2S0  m  resolution  for  the  analysis.  The  potential  temperature 
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profiles  also  have  low  horizontal  resolution  because  they  are  extracted  from 
smoothed  fields  from  the  NWS  observing  network.  Notice  the  minima  in  Ri  in 
the  regions  of  maximum  shear  above  and  below  the  peak  of  the  jet.  Pilot 
reports  of  the  altitude  and  severity  of  turbulence  for  the  western 
Pennsylvania  region  were  compiled  for  the  period  of  the  event  (Fig.  8).  The 
PIREPS  were  most  closely  correlated  with  the  magnitude  of  the  shear  observed 
by  the  radar,  rather  than  the  occurrences  of  low  values  of  Ri  (Fig.  9).  This 
figure  unequivocally  shows  that  regions  of  strong  shear  are  almost  sure  to 
have  severe  turbulence.  This  is  consistent  with  (19a)  which  implies  the 
strongest  turbulence  in  the  thickest  layers  (largest  Lq)  and  the  regions  of 
greatest  wind  shear  (3U/3z) .  Also,  recall  that  (29)  implies  that  the  thickest 
layers  are  also  the  longest  lived  and,  therefore,  the  most  likely  to  be 
encountered  by  aircraft.  This  result  suggests  that,  as  far  as  severe 
turbulence  is  concerned,  it  is  probably  a  mistake  to  attempt  to  derive  the 
Richardson  number  from  low  resolution  sources,  partly  because  regions  of 
strong  shear  are  likely  to  have  small  Ri  anyway. 

Preliminary  results  from  the  comparisons  of  the  profiler  winds  with  the 
NGM  winds  have  been  surprising.  The  profiler  winds  have  been  compared  with 
analyzed  (0  hour  'forecast')  and  predicted  (48  hour  forecast)  NGM  winds.  The 
NGM  data  is  interpolated  to  the  profiler  location  (both  Crown  and  McElevys 
Fort  sites  were  used).  The  profiler  data  has  been  examined  with  several  time 
resolutions:  the  basic  one  hour  average  at  the  time  of  the  NGM  data,  a  seven 
hour  average  (plus  and  minus  three  hours  about  the  time  of  the  observation) , 
and  a  thirteen  hour  average.  The  comparison  is  best  for  the  longer  time 
average,  indicating  the  amount  of  spatial  averaging  inherent  in  the  regional 
scale  model  or  smoothing  of  extreme  observations  in  the  radar  data.  A  sample 
comparison  for  the  500  mb  wind  speed  and  direction  for  May,  1987,  is  shown  for 
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the  analysis  (Fig.  10)  and  the  48  hour  forecast  (Fig.  11).  The  initial 
analysis  is  typically  2.5  m/s  different  from  the  radar  data  and  there  appears 
to  be  no  significant  systematic  bias  in  the  wind  direction.  It  is  believed 
that  the  NGM  underestimates  the  winds  compared  to  the  radar  because  of 
smoothing  of  the  fields  in  the  analysis  procedure.  The  48  hour  forecast  is 
noticeably  more  inaccurate  than  the  analysis,  but  the  deterioration  is 
surprisingly  small. 
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VI.  Sueeested  Future  Research 

The  general  goal  of  this  research  is  Co  investigate  the  fundamental 

relationships  between  mean  atmospheric  structure  and  microturbulence 

2 

parameters  such  as  CR  .  We  wish  to  study  the  processes  that  lead  to 

2 

temporal/horizontal  variability  and  vertical  distribution  of  Cn  .  The 
next  obvious  step  in  this  investigation  is  to  acquire  and  analyze  data  from 
three  primary  sources: 

(1)  Existing  data.  Radar  data  from  VHF1  and  VHF2  archived  over  the  past 
two  years  will  be  available  for  completion  of  a  climatological 
analysis.  Raw  aircraft  turbulence  data,  which  is  archived  on  FM  tapes, 
and  raw  aircraft  mean  profile  data,  which  is  archived  on  digital  tapes 
with  2  second  resolution,  from  the  EWAK  experiment  could  be  reanalyzed 
with  finer  resolution. 

(2)  New  Profiler  data.  We  will  continue  to  archive  one  hour  average 
radar  data  from  VHF1,  VHF2,  and  VHF3  (which  is  to  be  in  operation  in 
September,  1988).  We  will  also  archive  one  hour  average  data  from  the 
UHF  radar,  which  is  to  begin  regular  operation  at  State  College  in 
July,  1988.  Also,  we  hope  to  begin  operation  of  the  mm-wave  radiometer 
system  in  State  College  about  July  of  1988.  This  system  will  archive 
two  minute  time  resolution  measurements  of  integrated  water 
vapor/liquid  and  temperature  profiles.  The  mesoscale  triangle  (about 
140  km  sides)  of  VHF  radar  combined  with  the  UHF  and  thermodynamic 
sounders  will  provide  an  unprecedented  look  at  atmospheric  turbulence 
and  dynamical  processes. 

(3)  Intensive  Turbulence  and  Optical  Propagation  Experiment  (EWAK  II). 
Sometime  in  the  near  future  it  makes  sense  to  have  another  field 
program  similar  to  EWAK.  We  would  try  to  work  with  the  same 
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institutions  as  before  for  the  optical  measurements.  Ue  could  also  try 
to  interest  and  invite  NASA  and  NOAA  groups  with  relevant  measurement 
systems  (e.g.,  Melfi's  raman  lidar) . 

The  variety  of  scientific  issues  of  interest  were  discussed  in  section 
III.  Ue  can  divide  the  issues  into  four  general  categories:  spatial 
structure,  turbulence  dynamics,  climatology/synoptic  correlations,  and 
comparisons  of  optical,  radar  and  in  situ  measurements. 

A.  Spatial  Structure 

The  vertical  distribution  of  the  active  turbulence  areas  can  be  studied 

with  combined  aircraft  and  radar  data.  Ue  are  interested  in  the  probability 

distributions  of  patch  size,  patch  lifetime,  turbulence  scales, and  the 

autocorrelation  functions  of  patch  distribution.  We  can  also  use  aircraft 

measurements  to  examine  the  probability  distributions  of  lapse  rate  and  shear 

(basic  parameterizations  in  the  Van  Zandt  model) .  The  horizontal  and  temporal 

persistence  of  individual  turbulent  layers  can  be  studied  by  cross  correlation 
o 

the  Cn  data  from  the  three  VHF  radars  and  by  examining  level  flight 
aircraft  data. 

B.  Turbulence  Dynamics 

Here  the  emphasis  is  on  the  microturbulence  scaling  relations  (e.g.,  the 
2  2 

ratio  CT  /Cu  ),  the  generation  and  decay  process,  and  the  turbulence 

size  scales  (L^  and  LQ) .  Aircraft,  radar  and  thermosonde  data  can  all  be 

used.  Ue  suggest  Investigating  the  use  of  the  radar  to  measure  t  from  the 

width  of  the  Doppler  spectrum.  Furthermore,  we  can  compare  «  derived  from  the 

width  of,  for  example,  a  300  m  range  resolution  spectrum  with  that  from  the 

mean  shear  derived  from  overlapping  100  m  range  resolution  data.  Thus,  we  can 

2  2 

examine  simultaneous  measurements  of  C,p  and  Cu  .  The  size  scales  of 
turbulence  can  be  investigated  in  several  ways .  The  operation  of  colocated 
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VHF  and  UHF  radar  will  also  provide  information  about  the  relative  intensity 
of  turbulence  size  scales  since  they  scatter  from  different  atmospheric  scales 
(A/2“3  m  and  0.3  m  respectively).  Profiles  of  mean  temperature  from 
aircraft,  rawinsonde,  thermosonde,  and  radiometers  would  give  an  indication  of 
N  which,  combined  with  «, gives  L^.  The  vertical  velocity  spectrum  from  the 
radar  can  also  give  an  indication  of  size  scales.  In  active  turbulent  layers 
it  should  reveal  Lo;  in  more  quiescent  layers  it  should  indicate  N.  In  the 
intensive  experiment,  a  great  deal  of  effort  could  be  devoted  to  examining 
these  issues. 

C.  Climatology/Synoptic 

The  mesoscale  triangle  of  VHF  radar  will  provide  high  quality  dynamical 

variables  usually  not  available  (e.g,  divergence  and  vorticity)  that  will  be 

2 

extremely  useful  in  interpreting  variability  in  .  For  example,  the  jet 
stream  study  referred  to  in  section  IV  would  have  greatly  benefited  from  a 
more  quantitative  classification  of  conditions  than  the  crude  'distance  from 
the  jet  axis'  categories  used.  The  temperature/humidity  radiometers  also 
provide  very  accurate  estimates  of  dynamic  height  (e.g.,  500  mb)  to  augment 
the  radar  information.  We  presently  plan  to  continue  the  large  scale  synoptic 
classifications  for  the  gross  climatological  analysis.  We  will  also  continue 
the  studies  of  the  comparisons  of  radar  data  and  NGM  analyzed  and  forecast 
fields. 

D.  Comparison  of  Optical  and  Meteorological  Sensors 

We  strongly  recommend  that  another  joint  meteorological  and  optical 
experiment  be  held  a  Penn  State.  EWAK  I  revealed  that  there  are  still 
substantial  differences  between  the  various  observing  systems  (often  an  order 
of  magnitude).  Some  of  this  is  due  to  differences  in  temporal  and  spatial 
averaging  techniques  but  a  great  deal  of  the  disagreement  must  be  due  to 
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undiagnosed  problems  with  the  sensors  and  data  processing.  A  new  intensive 
program  would  also  provide  us  with  other  sources  of  data  (e.g.,  thermosondes 
and  aircraft)  and  much  better  spatial  and  temporal  resolution  of  the  important 
variables.  This  is  necessary,  in  particular,  to  determine  the  importance  of 
the  averaging  process  on  the  interpretation  of  measurements  and  models. 


-  B-39  - 


VII .  References 

Barat,  J.  and  F.  Bertin,  1984:  Simultaneous  measurements  of  temperature  and 
velocity  fluctuations  within  clear  air  turbulence  layers:  Analysis  of 
the  estimate  of  dissipation  rate  by  remote  sensing  techniques.  2-  Atmos . 
Sci. ,  41,  1613-1619. 

Beecher,  E.A.,  1987:  Analysis  of  temperature  and  velocity  microturbulence 
parameters  from  aircraft  data  and  relation  to  atmospheric  refractior 
index  structure.  M.S.  thesis,  Pennsylvania  State  University,  ppl65. 

Carlson,  C.A. ,  1987:  Kinematic  quantities  derived  from  a  triangle  of  VHF 
Doppler  wind  profilers.  M.S.  thesis,  Pennsylvania  State  University. 

Dillon,  T.M. ,  1982:  Vertical  overturns:  A  comparison  of  Thorpe  and  Ozmidov 
length  scales.  Geophvs .  Res .  .  87,  9601-9613. 

Fairall,  C.W. ,  K.L.  Davidson  and  G.E.  Schacher,  1982:  Meteorological  models 
for  optical  properties  in  the  marine  atmospheric  boundary  layer.  Opt . 
EnEineer. .  21,  847-857. 

Fairall,  C.W.  and  R.  Markson,  1984:  Aircraft  measurements  of  temperature  and 
velocity  microturbulence  in  the  stably  stratified  free  troposphere. 

Proc.  Seventh  Symposium  on  Turbulence  and  Diffusion,  AMS,  Boulder,  CO. 

2 

Fairall,  C.W.,  1987:  A  top-down  and  bottom-up  diffusion  model  of  C^.  and 
2 

in  the  entraining  convective  boundary  layer,  i.  Atmos .  Sci . , 

44,  1009-1017. 

Fairall,  C.W.,  D.W.  Thomson,  and  R.  Markson,  1988:  An  aircraft  and  radar  study 
of  temperature  and  /elocity  microturbulence  in  the  stably  stratified  free 
troposphere.  Proc.  Eighth  Symposium  on  Turbulence  and  Diffusion,  AMS, 

San  Diego,  CA,  61-65. 

Fried,  D.L. ,  1966:  Optical  resolution  through  a  randomly  inhomogeneous  medium 
for  very  long  and  very  short  exposures,  i.  Opt .  Soc .  Am. ,  56 .  1372-1380. 


-  B-40  - 


Kennedy,  P.J.  and  M. A.  Shapiro,  1975:  The  energy  budget  in  a  clear  air 

turbulence  zone  as  observed  by  aircraft.  Mon.  Wea .  Rev. .  103 .  650-654. 

Kennedy,  P.J.  and  M.A.  Shapiro,  1980:  further  encounters  with  clear  air 
turbulence  in  research  aircraft.  2.  Atmos .  Scl . ■  12,  986-993. 

Knovlton,  L.W. ,  1987:  Kinematic  diagnoses  of  frontal  structure  and 

circulations  derived  from  two  and  three-station  VHF  Doppler  wind 
profilers  networks.  M.S.  thesis,  Pennsylvania  State  University. 

Livingston,  P.M.,  1972:  A  study  of  target  edge  response  viewed  through 
atmospheric  turbulence  over  water.  AppI .  Opt. .  H,  2352-2357. 

Moss,  M.T.,  1986:  Measurement  and  modeling  of  tropospheric/stratospheric 
refractive  index  structure  parameter.  PhD  thesis  prospectus,  Florida 
State  University. 

Nastrom,  G.D.,  B.B.  Balsley  ,  and  K.S.  gage,  1981:  Change  with  season  of 
2 

Cn  at  Poker  Flat,  Alaska,  from  MST  Doppler  radar  observations. 

Proc.  20th  Conf.  on  Radar  Meteorology,  AMS,  Boston,  MA. 

Neiman,  P.J.,  1987:  Wind  profiler  derived  temperature  gradients  and 
advections,  M.S.  thesis,  Pennsylvania  State  University. 

Nieuwstadt,  F.T.M. ,  and  R.A.  Brost,  1986:  The  decay  of  convective  turbulence. 

2.  Atmos .  Sci. .  43.  532-546. 

Syrett,  W.J.,  1987:  Some  applications  of  50  MHz  wind  profiler  data:  detailed 
observations  of  the  Jet  stream.  M.S.  thesis,  Pennsylvania  State 
University,  ppl35. 

Syrett,  W.J.  and  D.W.  Thomson,  1988:  Detailed  observations  of  the  jet  stream 
with  the  Penn  State  VHF  Doppler  radar.  J..  Appl.  Met.  ,  in 

preparation. 

Tatarski,  V.I.,  1961:  Wave  Propagation  in  a  Turbulent  Medium.  McGraw-Hill. 
New  York,  Chapter  7. 


-  B-41  - 


Thomson,  D.W, ,  W.J.  Syrett,  T.T.  Warner,  and  N.L.  Seaman,  1988:  Comparisons  of 
wind  Profiler  measurements  with  NMC  NGM  analyses  and  predictions.  Proc. 
8th  Conf.  on  Numerical  Weather  Prediction,  AMS,  Feb.  22-25,  Baltimore, 

MD ,  pp6 . 

Van  Zandt,  T.W. ,  J.L.  Green,  K.S.  Gage,  and  W.L.  Clark,  1978:  Vertical 

profiles  of  refractivity  turbulence  structure  constant:  comparison  of 
observations  by  the  sunset  radar  with  a  new  model.  Radio  Sci . , 

1,  819-829. 

Van  Zandt,  T.W.  ,  K.S.  Gage  and  J.M.  Wamock,  1981:  An  improved  model  for  the 

o 

calculation  of  profiles  of  Cn  and  «  in  the  free  atmosphere  form 
background  profiles  of  wind,  temperature  and  humidity.  Proc.  20th  Conf. 
on  Radar  Meteorology,  AMS,  Boston,  MA. 

Wesely,  M.L. ,  1976:  The  combined  effect  of  temperature  and  humidity 
fluctuations  on  refractive  index.  2.  AppI .  Meteor. .  15 .  43-49. 


-  5-42  - 


VIII. 


The  following  is  a  list  of  publications  and  presentations  at  least 
partially  funded  by  or  directly  concerned  with  this  research  project. 


A.  Publications 

Williams,  S.R.  and  D.V.  Thomson,  1986:  An  evaluation  of  errors  observed  in  the 


measurement  of  low  wind  velocities.  Handbook  for  MAP.  Vi 
secretariat,  Dep.  Elec.  Comp.  Eng.,  Univ.  IL. ,  256-262. 


,  SCOSTEP 


Williams,  S.R.,  and  R.  Peters,  1986:  The  Penn  State  Doppler  network  progress 


report.  Handbook  for  MAP.  Vol. 
Comp.  Eng.,  Univ.,  IL. ,  339-341. 


,  SCOSTEP  secretariat,  Dep.  Elec. 


Moss,  M.T. ,  1986:  Measurement  and  modeling  of  tropospheric/stratospheric 
refractive  index  structure  parameter.  PhD  thesis  prospectus,  Florida 
State  University. 

Beecher,  E.A. ,  1987:  Analysis  of  temperature  and  velocity  microturbulence 
parameters  from  aircraft  data  and  relation  to  atmospheric  refraction 
index  structure.  M.S.  Thesis,  Pennsylvania  State  University,  ppl65. 
Syrett,  W.J.,  1987:  Some  applications  of  50  MHz  wind  profiler  data:  detailed 
observations  of  the  jet  stream.  M.S.  thesis,  Pennsylvania  State 
University,  ppl35. 

2 

Fairall,  C.W. ,  1987:  A  top-down  and  bottom-up  diffusion  model  of  C^.  and 
2 

0^  in  the  entraining  convective  boundary  layer.  2.  Atmos .  Sci . . 

44,  1009-1017. 

Thomson,  D.W. ,  W.J.  Syrett,  T.T.  Warner,  and  N.L.  Seaman,  1988:  Proc.  8th 
Conf.  on  Numerical  Weather  Prediction,  AMS,  Feb.  22-25,  Baltimore,  MD, 


Fairall,  C.W.,  D.W.  Thomson  and  R.  Markson,  1988:  An  aircraft  and  radar  study 


-  B-43  - 


of  temperature  and  velocity  microturbulence  in  the  stably  stratified  free 
troposphere.  Proc.  8th  Symposium  on  Turbulence  and  Diffusion,  AMS,  April 
24-29, San  Diego,  CA, 61-65. 

Williams,  S.A.  and  D.W.  Thomson,  1988:  Comparisons  between  Wind  Profiler  and 
conventional  Loran/Omega  based  rawinsonde  wind  measurements.  In 
preparation. 

B.  Oral  presentations 

"Clear  air  Doppler  radar",  C.W.  Fairall,  RISO  National  Laboratory  (Denmark), 
Aug.  6,  1986. 

"Turbulence  measurements  with  Doppler  profilers",  D.W.  Thomson,  Naval 
Environmental  Research  and  Prediction  Facility,  June  15.  1986. 

"Wind  and  turbulence  measurements  with  a  clear  air  Doppler  radar",  C.W. 

Fairall,  Florida  State  University,  Aug.  28,  1986. 

"An  analysis  of  synoptic  conditions  during  EWAK" ,  W.  Syrett,  AFGL,  Sept.  4, 
1986. 

"Radar  measurements  of  Cn*2  during  EWAK",  M.  Moss,  AFGL,  Sept.  4,  1986. 

"Sodar  measurements  of  CnA2  during  EWAK",  T.  Messier,  AFGL,  Sept.  4,  1986. 
"Aircraft  measurements  of  microturbulence  during  EWAK",  C.  Fairall,  AFGL, 

Sept.  4,  1986. 

"Wind  Profilers",  D.W.  Thomson,  a  series  of  three  lectures,  Dept,  of 

Meteorology,  Naval  Postgraduate  School,  Monterey,  CA. ,  Jan.,  1987. 

2 

"Toward  a  physical  model  of  C  ",  C.  Fairall,  Workshop  on  the  Physics  of 

Directed  Energy  Propagation  in  the  Atmosphere,  Las  Cruces,  NM,  March  29, 
1988. 

"Atmospheric  radar  or  seeing  with  thick  air",  C. Fairall,  Dept,  of  Physics, 
Naval  Postgraduate  School,  Monterey,  CA,  April  20,  1988. 


OEWAA  RADAR  CN2 


OEWAA  RADAR  CN2 

uro2S«/04Z 

-16.4 

-te.e 
-16.0 
-17 
-17.2 
-17.4 

N 

§  -17.6 

|  -i7.a 

-18 
-ia.2 
-ia4 
-1»6 
-iaa 

-i* 


2 

Figure  1.  Sample  high  resolution  time  series  of  Cn  at  two  different 
range  resolutions.  The  upper  panel  is  with  0.9  km  vertical  resolution  and  3 
minute  time  resolution;  the  lower  panel  is  with  0.3  km  vertical  resolution  and 
90  second  time  resolution.  Each  panel  gives  a  time  series  for  two  different 
altitudes . 
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Figure  2.  Sample  profile  of  optical  C  2  from  the  EWAK  experiment.  The  * 
symbol  denotes  the  radar  estimate  (corrected  for  moisture  effects)  with 
horizontal  bars  to  indicate  the  standard  deviation  over  an  hour.  We  now 
believe  the  radar  data  should  be  increased  by  at  least  a  factor  of  two  from 
the  values  indicated  in  this  graph.  The  solid  line  is  AFGL  thermosonde  data 
degraded  to  the  radar  resolution  by  linear  averaging.  The  open  circles  are 
from  the  aircraft  measurements  using  the  spectral  method.  The  designator  in 
the  lower  left  corner  indicates  the  profile  is  for  May  4,  1986  at  0600  GMT. 
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Figure  4.  Comparison  of  various  measurements  of  optical  C  profiles  from 
EWAK.  The  solid  line  is  the  Van  Zandt  model,  the  dashed  line  is  the  radar, 
the  solid  line  with  horizontal  bars  is  the  averaged  scintillometer,  the  solid 
line  with  tic  marks  is  the  aircraft,  and  the  dotted  line  is  the  thermosonde. 
The  upper  panel  is  at  about  0000  local  time  on  U  May  1986;  the  lover  panel  is 
about  2300  local  time  on  5  May  1986. 
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Figure  5a.  A  comparison  of  simultaneous  values  of  temperature  and  velocity 
microturbulence  when  normalized  using  equation  13.  A  slope  of  1  would 
indicate  that  7^  was  a  constant. 
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Figure  5b.  A  composite  comparison  of  normalized  by  Cu  and 
potential  temperature  gradient  (as  in  Eqs.  13  and  30)  as  a  function  of  the 
activity  parameter,  A.  A  constant  value  would  indicate  that  7^  was  a 
constant. 


Figure  7a.  Average  wind  speed  profiles  with  the  measurement  site  within  100 
km  of  the  jet  axis.  The  horizontal  bars  indicate  the  standard  deviation  for 
the  period.  The  upper  panel  is  case  1;  the  lower  panel  is  case  2. 
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CASE  2  PILOT  REPORTS  OF  TURBULENCE  3X7  DEGREE  BOX 
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Figure  8.  Piloe  reports  of  turbulence  in  the  vicinity  of  the  radar  during 
case  2.  The  bottoa  panel  is  a  blowup  of  the  boxed  area  of  the  top  plot. 


Figure  9.  Surface  plots  of  v 
1987.  The  surface  defines  th 
The  shaded  areas  indicate  the 


-  B-54  - 


CROWN  500mb  WIND  DIRECTION,  MAY  1987 
NGM-00  vs  1-HR  PROFILER 
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Figure  10.  A  comparison  of  Che  NGM  interpolated  and  radar  derived  500  mb  wind 
direction  (upper  panel)  and  wind  speed  (lower  panel)  for  Che  month  of  May, 
1987.  The  NGM  data  represents  the  smoothed  analysis  field  used  to  initialize 
the  model  (0  forecast). 
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CROWN  WIND  DIRECTION,  MAY  1987 
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Figure  11.  As  in  figure  10  but  for  the  48  hour  NGM  forecast. 
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APPENDIX  A 

Profiler  System  Description 

Mesoscale  meteorological  measurements,  analysis  and  prediction  are  some 
of  the  principal  areas  of  research  in  the  Department  of  Meteorology  at  Penn 
State.  For  more  than  a  decade  those  members  of  the  faculty  concerned  with 
mesoscale  analysis,  numerical  modeling  and  forecasting  have  been  frustrated  by 
the  spatial  and  temporal  inadequacy  of  conventional  network  observations  for 
both  research  and  operational  applications.  For  more  than  five  years  the 
Department  had  sought  the  substantial  financial  resources  required  to  deploy  a 
network  of  VHF  Doppler  (ST)  radars  and  millimeter  wave  radiometers  for 
"operational"  test  and  evaluation  for  wind  and  thermodynamic  profiling. 
Construction  of  the  ST  radar  network  began  in  fall  of  1983  using  funding 
provided  by  the  Air  Force  Office  of  Scientific  Research  (through  the  DoD 
University  Research  Instrumentation  Program)  and  the  University.  In  1986  a 
second  AFOSR  URIP  grant  was  obtained  to  fund  the  acquisition  of  temperature 
(50-60  GHz)  and  humidity  (20-30  GHz)  profiling  radiometers.  This  document 
will  only  discuss  the  radar  systems. 

For  the  foreseeable  future  the  Penn  State  ST  radar  program  will  be 
focused  on  applications  rather  than  systems  development  research.  Deployment 
of  the  systems  would  not  have  been  possible  without  the  outstanding 
cooperation  provided  by  C.G.  Little  and  R.  Strauch  and  their  colleagues  at  the 
Wave  Propagation  Laboratory,  and  also  J.  Brosnahan  of  Tycho  Technology  from 
whom  we  bought  all  of  the  receivers,  transmitters  and  antennas.  With  regard 
to  the  other  major  systems  components,  we  have  assembled  in-house,  from  WPL 
documentation,  the  time-domain- integrator  and  computer  interfaces  and  have 
purchased  WPL  software -compatible  Data  General  Corp.  Eclipse  computers  for 
each  system. 
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The  Penn  State  network  consists  of  three  6  m  wavelength (VHF)  and  one  0.7 
m  wavelength(UHF)  radars.  Fig.  1  indicates  the  approximate  location  of  a 
mesoscale  triangle  formed  by  the  three  VHF  radar  sites  within  the  routine 
ravinsonde  network.  In  Table  1,  a  short  summary  of  the  specifications  for  the 
four  radars  is  given.  Experiences,  plans  and  improvements  for  the  PSU  network 
are  summarized  below. 

A.  VHF1  50-MHz  radar  located  15  km  south  of  State  College,  PA. 

1.  This  system  became  fully  operation  June  27,  1985.  The  primary  reason 
for  system  failures  since  the  onset  has  been  AC  outages  which  are  prevelant  in 
this  area.  Battery  back-up  and  computer-controlled  autorestart  of  the 
transmitters  has  circumvented  this  problem. 

2.  Initial  performance  statistics  done  by  Frisch  et  al.,  WPL  on  August 
1985  data  ,  indicate  very  good  performance  by  VHF1.  On  Beam  #1  the  next  to 
last  range  gate  (16.8  km  MSL)  was  able  to  make  a  wind  measurement  99%  of  the 
time,  while  Beam  #2  was  able  to  measure  the  wind  at  this  height  85%  of  the 
time.  The  difference  between  the  beams  is  probably  due  to  better  pickup  by 
beam  #2  of  computer/electronic  noise  from  the  building.  Samples  of  data  from 
the  low  altitude  range  of  VHF1  are  shown  in  Figs.  2  and  3. 

3.  A  vertical  beam  was  added  in  March,  1987,  and  the  latest  NOAA/VPL 
software  was  installed. 

B.  VHF2  50-MHz  radar  located  in  NW  Pennsylvania  near  Crown,  PA. 

1.  This  system  was  installed  May  1,  1986  and  has  operated  since. 
Performance  is  often  slightly  better  than  VHF1  but  a  source  of  interference 
(radiotelephones  from  a  local  trucking  company)  leads  to  periods  of  reduced 
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data  quality. 

2.  A  vertical  bean  was  added  in  January,  1987,  and  the  latest  NOAA/WPL 
software  was  installed. 

C.  VHF3  50 -MHz  radar  located  in  SV  Pennsylvania  near  Somerset,  PA. 

1.  All  hardware  is  ready  for  installation  pending  selection  and 
preparation  of  a  site. 

2.  A  three  beam  system  similar  to  VHF1  and  VHF2  will  be  used.  Some  of 
the  hardware  from  Tycho  Technology  has  been  slightly  updated. 

D.  UHF1  portable  405 -MHz  radar  to  be  semi -permanently  based  at  the  PSU 
Circleville  Farm  in  downtown  State  College,  PA. 

1.  This  system  has  been  used  on  three  major  field  deployments:  the 
SPACE/MI ST -COHMEX  experiment  in  Alabama  (June-July,  1986),  the  FIRE 
stratocumulus  IFO  on  San  Nicolas  Island,  CA  (July,  1987),  and  the  Arizona 
monsoon  experiment  (August,  1987). 

2.  Preparations  for  basing  at  Circleville  Farm  are  to  be  completed  in 


July,  1988. 
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Table  1:  Specifications  for  the  Penn  State  ST  Radar? 


Item 

VHF1.VHF2.VHF3 

UHF1 

Type 

Pulsed  Doppler,  3  beam 

Pulsed  Doppler,  3  beam 

Location 

1:  S.  State  College 

2 :  Crown ,  PA 

3:  Somerset, PA 

In  State  College,  trans 
portable 

Frequency  (non.) 

49.9  MHz 

404  MHz 

Bandwidth 

300,  100  KHz 

1,  0.3,0.07  MHz 

Peak  Power 

30  KW 

30  KW 

Pulsewidth 

3.67,  9.67  psec 

1,  4,  16  fisec 

Antenna : 

Type 

Dimensions 

Zenith  angles 

Phased  array  CoCo 

50  m  X  50  m 

90* (vert) ,75* 

Phased  array  CoCo 

9  m  X  9  m 

90* (vert) ,  75* 

One  site  computer 

Data  Gen.  Eclipse 

Data  Gen.  Eclipse 

On  site  processing* 
at  PW- 

3.67  ^s  9.67  ns 

1/js  4/xs  16/is 

Time  domain  ave 
Spectral  ave 

Height  spacing  (m) 
Spec .  resolut . (m/s) 
Maximum  absolute  rad. 
velocity  (m/s) 

400  125 

8  16 

290  870 

0.49  0.31 

15.7  19.6 

112  70  35 

16  32  64 

100  300  800 

0.29  0.29  0.29 

18.25  18.25  18.25 

♦Software  controllable 


FNT 


_e  1^  Locations  of  Penn  State  radar  sites  (open  circles)  and  standard 
rawinsonde  sites  (solid  circles). 
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figure  3  Two  dimension  plot  or  wind  speed  versus  Urne  and  altitude  for  the  same 
period  as  Fig.  2. 
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NOTES  CONCERNING  THE  USE  OF  CLEAR  AIR  DOPPLER  RADARS 
FOR  THE  MEASUREMENT  OF  C„ 2 

Robert  M.  Peters 

Department  of  Meteorology 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Recently  the  Penn  State  clear  air  radar  network  has  been 
increasingly  utilized  not  only  for  wind  profiling,  but  also  for  the 
measurement  of  the  refractive  structure  coefficient  associated  with 
isotropic  turbulence,  Co2 .  This  coefficient  is  related  to  received 
power  amplitude  of  active  indirect  sensing  systems  when  signal 
backscatter  is  caused  solely  by  refractive  index  variations  due  to 
isotropic  turbulence.  The  physical  mechanisms  associated  with  this 
backscatter  in  the  atmosphere  are  described  in  detail  by  WESELY 
(1976)  . 

The  present  Penn  State  clear  air  radar  network  consists  of  two 
VHF  ( 6M  wavelength)  Doppler  radars  and  a  transportable  UHF  (75  CM 
wavelength)  Doppler  radar.  The  systems  are  designed  primarily  for 
wind  profiling  of  the  troposphere  and  stratosphere.  The  Penn  State 
network  is  described  by  THOMSON  et  al .  (1984)  and  is  based  upon  the 

concepts,  signal  processing  hardware  designs  and  software  provided  by 
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the  NOAA/ERL  Wave  Propagation  Lab  (STRAUCH  et  al.,  1984).  An  ongoing 
experimental  research  effort  at  Penn  State  pertains  to  measurements  of 
turbulent  backscatter  from  radars,  optical  systems,  acoustic  systems 
and  direct  measurements  of  turbulence  via  aircraft  and  thermosondes. 
During  a  recent  experiment,  one  VHF  radar  was  used  for  this  purpose 
and  measurements  of  the  VHF  structure  coefficient  were  obtained  (MOSS, 
1986) .  It  has  become  evident  that  the  quantitative  measurement  of 
radar  system  power  parameters  that  are  required  for  the  determination 
of  CB 2  are  prone  to  errors  that  are  frequently  not  addressed  in 
designs  that  are  optimized  for  Doppler  measurements. 

In  order  for  backscattered  power  to  be  related  to  isotropic 
turbulence,  it  must  be  known  that  the  received  signal  amplitude  is  not 
modified  by  other  backscatter  sources,  reflections  or  interference 
from  other  emitters.  Continuous  interference  is  usually  man  made  from 
sources  such  as  arcing  power  lines,  radio  transmissions  or  broadband 
noise  from  the  digital  systems  of  the  radar.  Good  site  selection  and 
radar  system  engineering  practices  will  prevent  interference  from 
these  types  of  emission.  Hard  targets  will  cause  undesired 
reflections.  Reflections  due  to  ground  clutter  may  be  partially 
eliminated  in  spectral  processing  but  are  best  avoided  by  choosing 
radar  sites  that  have  no  large  targets  visible  at  distances  comparable 
to  the  radar's  range  response.  Hard  reflections  from  aircraft  will 
override  the  atmospheric  backscattered  signal.  The  consensus 
algorithm  used  for  averaging  in  the  Penn  State  radar  software  (STRAUCH 
et  al.,  1984)  will  eliminate  the  short  term  occurrence  of  aircraft 
reflections  from  an  averaged  data  set.  However,  the  location  of  major 
airways  should  be  a  factor  during  radar  site  selection.  Enhanced 
signal  returns  from  the  atmosphere  from  other  natural  sources  of 
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anisotropic  irregularities  will  occur  with  a  zenith  pointing  beam  as 
described  by  DOVIAK  and  ZRNIC'  (1983).  As  shown  by  DOVIAK  and  ZRNIC' 
(1983),  the  response  of  VHF  radars  to  anisotropic  irregularities  will 
be  insignificant  at  a  15  degree  off  zenith  pointing  angle,  assuming 
that  the  layered  structure  causing  the  signal  enhancement  is 
horizontal.  At  Penn  State,  one  of  the  15  degree  off-axis  beams  of  the 
VHF  radar  antenna  is  used  for  the  measurement  of  Cn2  as  well  as  for 
one  wind  velocity  component.  Antenna  beam  side  lobe  response  to 
enhanced  signal  returns  from  the  vertical  are  further  reduced  by 
having  a  zenith  null  designed  into  the  antenna  pattern. 

In  order  to  make  quantitative  measurements  of  received  signal 
amplitudes,  the  radar  systems 's  hardware  must  be  calibrated  or  its 
calibration  inferred  from  other  sources.  System  hardware  for 
calibration  purposes  are  the  antenna  system,  transmitters,  receivers 
and  analog  to  digital  converters  used  for  parameter  measurements. 

The  antenna  patterns  of  large  upward  pointing  VHF  phased 
arrays,  such  as  those  used  at  Penn  State,  are  difficult  to  measure  in 
their  far  fields.  Computer  antenna  pattern  models  are  quite  accurate 
with  regard  to  antenna  pattern  shape.  Experience  has  shown  that  an 
antenna  will  usually  have  a  gain  within  a  few  decibels  (dB)  of  that 
predicted  by  a  validated  antenna  model,  if  the  antenna  is  well 
engineered  and  installed.  Pattern  altering  metallic  objects  should 
not  be  placed  in  or  near  the  antenna  field  if  reliance  is  placed  upon 
previously  determined  antenna  gain  figures  that  are  to  be  used  in  the 
measurement  of  structure  parameters.  Astronomical  noise  sources  have 
been  used  to  confirm  antenna  pointing  maxima.  Galactic  sources  are 
useful  in  the  northern  hemisphere  (MOSS,  1986)  and  the  sun  can  be  used 
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in  equatorial  regions  (RIDDLE,  1985)  .  Precise  patterns  and  gains  have 
been  determined  for  the  Mu  radar  facility  in  Japan  by  integrating 
direct  satellite  measurments  as  described  by  SATO  et  al  (1985) .  The 
phase  and  power  distribution  provided  to  the  elements  of  the  VHF 
phased  array  may  be  directly  measured  to  provide  a  more  accurate  input 
to  modeled  antenna  patterns.  The  performance  of  an  uncalibrated 
antenna  system  may  be  statistically  compared  with  other  systems  of 
known  performance  and  its  gain  inferred  (MOSS,  1986),  as  long  as  other 
system  components  are  calibrated.  Computer  estimation  of  Cn 2  may  be 
derived  from  atmospheric  soundings  and  used  to  estimate  system 
calibration.  A  field  tested  computer  model  developed  by  VANZANDT  et 
al .  (1981)  could  possibly  be  used  for  this  purpose. 

Transmitter  parameters  such  as  power  and  pulse  width  must  be 
known.  The  Penn  State  network  incorporates  transmitters  which  have 
internal  microprocessor  control  and  monitoring  systems.  The  radar 
system  executive  control  computer  at  each  site  provides  a  means  to 
remotely  monitor  and  control  each  radar  transmitter  in  the  network. 
Transmitter  power  may  be  monitored  via  this  arrangement,  however,  this 
indirect  method  of  measurement  is  prone  to  error  if  the  associated 
electronic  circuits  are  not  calibrated.  With  Penn  State  systems,  the 
transmitted  power  level  is  derived  from  an  adjustable  RF  voltage 
monitoring  circuit  with  its  output  supplied  to  the  analog  to  digital 
converter  of  each  transmitter’s  processor.  Each  transmitter  includes 
an  integral  directional  coupler.  To  ensure  the  accuracy  of  the 
transmitter  power  monitor,  the  coupling  coefficient  of  each 
directional  coupler  should  be  calibrated  and  power  measurements  taken 
to  check  and  calibrate  the  microprocessor  based  power  monitoring 


system . 
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The  receiver  system  includes  not  only  the  radar  receivers  and 
preamplifiers,  but  also  the  analog  to  digital  converters.  For 
structure  parameter  measurements,  the  receiving  system  gain  must  be 
calibrated.  For  relative  power  measurements  of  received  signal 
amplitude,  the  calibration  must  be  done  with  test  instrumentation. 
Any  changes  to  the  receiving  system,  such  as  replacement  of  a  receiver 
preamplifier,  would  require  recalibration.  For  radar  systems  such  as 
the  Penn  State  UHF  Doppler  radar,  this  is  the  most  practical  method  of 
calibration  since  the  UHF  system  utilizes  one  receiver  time  shared 
among  several  antenna  pointing  angles .  The  minimun  discernable 
received  signal  power  of  the  UHF  system  is  limited  by  the  100  K  noise 
temperature  (approximate)  of  the  receiver  preamplifier.  This  type  of 
calibration  is  unnecessary  for  the  VHF  Doppler  radars.  For  these 
systems,  the  minimum  signal  detectabilty  is  limited  by  a  galactic 
noise  temperature  backround  in  the  order  of  several  thousand  Kelvins 
(KO,  1958).  For  this  reason,  the  galactic  noise  level  is  used  as  an 
automatic  calibration  for  all  receivers  at  Penn  State  VHF  radar  sites. 
This  is  accomplished  by  archiving  a  spectrally  derived  estimate  of 
consensus  averaged  signal  to  noise  ratio  (most  systems  presently 
archive  relative  power  on  a  routine  basis) .  Structure  parameter 
errors  at  close  ranges  using  VHF  radar  have  been  observe  during  recent 
experiments  (MOSS,  1986).  The  cause  for  this  class  of  error  has  not 
been  studied  but  appears  either  to  be  due  to  a  near  field  antenna 
response  error  or  receiver  dynamic  range  limiting. 

Less  than  ideal  receiver  low  pass  filter  stability, 
calibration  and  adjustment  has  been  the  most  significant  cause  of 
error  of  spectrally  derived  structure  parameters  from  the  VHF  radars. 
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Specifically,  the  presence  of  DC  offsets  and  gain  imbalance  between 
real  and  quadrature  channels  will  produce  a  host  of  errors  in  spectral 
signal  processing.  Each  received  signal  spectrum  is  actually  the 
result  of  several  averaged  spectra,  each  of  which  is  derived  from  a 
time  series  of  coherently  integrated  received  signal  amplitudes  at 
each  radar  altitude.  Each  signal  spectrum  is  derived  from  the  complex 
Fast  Fourier  Transform  of  the  receiver’s  range  gated  real  and 
quadrature  outputs.  Peaks  that  may  occur  in  each  spectrum  about  the 
zero  frequency  line  at  an  equal  amplitude  at  all  range  gates  are 
usually  the  result  of  the  presence  of  a  mean  or  trend  in  the  complex 
signal  data  set  that  is  not  completely  removable  by  the  signal 
processing  software.  Image  peaks  that  may  occur  opposite  the  spectral 
peak  associated  with  the  clear  air  signal  return  usually  are  a  result 
of  dynamic  range  limiting  of  the  receiving  system.  A  consistent  image 
peak  in  all  range  gates  is  usually  associated  with  a  gain  or 
quadrature  imbalance  between  the  I  and  Q  channels  of  the  receiving 
system.  Note  that  the  receiving  system  also  includes  the  analog  to 
digital  converters  of  the  signal  processing  system.  Imbalances 
between  the  I  and  Q  analog  to  digital  converters  will  produce  the  same 
errors.  Errors  produced  by  means,  trends  and  imbalances  produce 
errors  in  structure  parameter  data  that  are  not  as  easily  removed  as 
compared  to  Doppler  velocity  measurements.  The  amplitude  of  a 
spectral  image  will  always  be  lower  than  the  true  signal  so  that 
Doppler  peak  detecting  algorithms  will  never  choose  the  image  peak  for 
Doppler  processing.  Any  means  or  trends  will  be  have  an  effect  in 
Doppler  velocity  processing  only  if  the  zero  amplitude  peak  is  greater 
than  the  clear  air  peak.  If  the  zero  peak  is  larger,  then  the 
artificially  constant  velocity  produced  in  the  beam  component  (which 
is  frequently  blamed  on  ground  clutter)  can  be  easily  removed  from 
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wind  velocity  data  sets  by  computer  post  processing.  The  effect  of 
these  types  of  spectral  errors  in  obtaining  received  signal  power  or 
signal  to  noise  ratio  are  not  as  obvious.  The  spectral  processing 
algorithm  will  include  ground  clutter  peaks  as  part  of  the  returned 
backscattered  signal  at  low  radial  wind  velocities.  Any  image  peaks 
that  occur  will  also  bias  the  estimates  of  received  signal  power  and 
possibly  produce  a  velocity  bias  to  the  estimate.  The  presence  of 
spectral  images  will  contribute  to  the  estimated  noise  level  thereby 
producing  errors  in  the  signal  to  noise  ratio. 

To  reduce  the  errors  associated  with  the  receiver  system  for 
all  types  of  measurements,  the  receivers  must  be  carefully  adjusted  in 
the  lab  for  gain  balance  and  zero  offsets  after  a  period  of 
stabilization.  Software  solutions  can  contribute  to  the  reduction  of 
error.  The  present  software  used  at  Penn  State  (STRAUCH  et  al.,  1984) 
removes  the  mean  from  the  received  signal  data  set.  Trend  removal  is 
more  difficult  for  the  real  time  nature  of  our  signal  processing 
system  however  improvements  have  been  suggested  by  STRAUCH  (1986)  . 
The  ultimate  solution  is  to  treat  the  disease,  not  the  symptoms.  The 
filter  and  output  stages  of  VHF  and  UHF  Doppler  receivers  should  be 
designed  for  maximum  amplitude  stability  as  well  as  frequency 
stability.  The  manufacturer  of  the  radio  frequency  subsystems  of  the 
Penn  State  Doppler  radars  are  now  providing  improved  receivers  for  new 
systems  and  retrofits  for  present  systems  to  reduce  the  possibility  of 
these  types  of  problems.  The  Penn  State  UHF  system  will  have  a 
programmable  attenuator  added  to  the  receiver  system.  This  will  allow 
the  receiver's  gain  to  be  fixed  at  an  ideal  level  depending  upon 
operating  requirements.  This  feature,  along  with  improved  on-line 
spectral  display  software,  will  allow  field  optimization  of  radar 
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parameters  to  reduce  the  sources  of  error  described  herein.  An 
improved  method  of  spectral  signal  detection  has  also  been  developed 
(PETERS  and  WILLIAMS,  1988) .  The  new  method  provides  a  better  signal 
peak  selection  mode  in  the  presence  of  ground  clutter  or  in  low  signal 
to  noise  conditions.  Careful  attention  to  clear  air  radar  system 
performance  will  result  in  better  estimates  of  all  measurable  clear 
air  parameters. 
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ABSTRACT 


Due  to  its  inherently  turbulent  nature,  the  atmosphere  has  a 

temporally  and  spatially  variable  refractive  index,  which  degrades 

propagating  electromagnetic  radiation.  The  refractive  index  turbulen 

2 

structure  constant,  CR  ,  is  a  key  parameter  for  describing 

2 

refractive  variations.  Cn  can  be  inferred  from  optical 

turbulence  (scintillometer),  temperature  turbulence  (aircraft 

instrumentation  or  thermosonde)  or  radar  backscatter  (profiler) 

2 

measurements .  Vertical  integrals  of  give  transverse  coherence 

length,  rQ,  and  isoplanatic  angle,  8  These  can  also  be  measured 

remotely  by  an  rQ  scintillometer  and  isoplanometer ,  respectively. 

During  an  atmospheric  optics/meteorology  experiment  (acronym 

EWAK)  conducted  at  Penn  State  University  primarily  during  April  and 

May  of  1986 ,  data  was  collected  by  all  of  the  aforementioned 

instruments .  An  instrumented  research  aircraft  was  used  to  measure 

2 

vertical  profiles  of  temperature  and  velocity  turbulence  (C^,  and 
2 

Cu  )  and  other  meteorological  variables.  The  turbulence 
instrumentation  aboard  the  aircraft  consisted  of  cold-wire,  and 
hot-wire  sensors  and  FM  recording  apparatus.  The  taped  data  was 
processed  via  FFT  to  produce  one -dimensional  variance  spectra 
(wavenumber  range  0.01  to  10  m"^).  Flights  usually  produced  a  10  km 
vertical  profile;  the  data  was  processed  to  give  roughly  0.5  km 
vertical  resolution  (similar  to  that  of  the  radar  wind  profiler). 

Editing  was  based  on  percent  error  between  a  regression  analysis 
and  theoretical  -5/3  frequency  dependence  of  the  spectra.  A  majority 
of  spectra  showed  good  evidence  of  the  classic  inertial  subrange  with 
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2  2 

-5/3  slope.  and  Cu  were  calculated  from  the  regression 

fit  to  the  inertial  subrange  power  spectral  density.  The  race  of 

dissipation  of  turbulent  kinetic  energy,  e,  was  calculated  using  the 

2 

so-called  Corrsin  relation  for  velocity  turbulence.  Optical 

2 

was  calculated  from  . 

Considerable  interest  has  developed  in  models  chat  relate 

microturbulence  parameters  to  the  mean  gradients.  One  such  model, 

proposed  by  VanZandt  et  al.  (1978;  1981),  has  as  a  key  variable  the 

gradient  Richardson  number,  R^.  In  actively  turbulent  regions, 

2  2 

can  be  related  to  C^,  and  Cu  as  follows: 

4  /  Cl  "  1-6  V  <?r  *  V  (*/g)  ( M/dz ) 

2  2 

A  scatter  plot  of  C^,  versus  (6/ g  dd/dz)  in  the  free 

troposphere  showed  a  range  of  values  from  0.3  to  10  for  high  and  low 

2  2 
values  of  ,  respectively.  Scatter  plots  of  C^.  versus 

2 

Cu  showed  high  correlation  of  these  two  parameters,  with  log 

2  .2 

Cu  -  1.7  log  .  R^-N/S  and  plots  of  N  versus  S  roughly 

corresponded  to  the  curve  depicted  by  Fairall  and  Markson  (1985). 

2  2 

Plots  of  /  Cu  (0/g  89/8 z)  versus  a  turbulent  activity 
parameter  (Gregg  1987)  clearly  showed  agreement  with  the  suggested 
activity  levels  and  associated  values.  At  high  activity  levels,  the 
ordinate  value  approached  0.4  (equivalent  to  1.6  when 
approximating  Pr  as  1  and  R^  as  0.25). 

2 

Values  of  N  and  S  were  used  to  obtain  values  of  C  and  e 

n 

2 

with  the  model  of  VanZandt  et  al.  Comparisons  of  Cn  profiles 
measured  by  the  various  instruments  showed  got  .  agreement  between 
scintillometer  and  thermosonde,  and  between  the  VanZandt  et  al .  model 


-  C-6  - 


and  the  aircraft  data.  The  profiler  (converted  optical  values)  had 
some  agreement  with  the  thermosonde.  Average  ratios  of  each  profile 
to  the  thermosonde  profile  were:  aircraft  8,  model  6.8,  profiler  8.1. 
and  the  scintillometer  1.1. 

2 

Numerical  integration  of  Cn  gave  estimates  of  rQ  and  8  . 

Values  from  scintillometer,  thermosonde  and  aircraft  did  not  show 

clear  agreement.  Based  on  the  aircraft  values,  the  scintillometer  and 

thermosonde  8  values  differed  by  an  average  of  17%.  The 

thermosonde  rQ  values  differed  by  an  average  of  35%. 

2 

Despite  good  agreement  on  Cn  ,  the  aircraft  data  and  the 
VanZandt  et  al.  model  output  clearly  disagreed  on  vertical  profiles  of 
e,  in  some  cases  by  several  orders  of  magnitude. 

Overall,  EWAK  provided  a  chance  to  compare  methods  of  measuring 

2 

Cn  ,  as  well  as  an  opportunity  to  obtain  microturbulence  data  for 
examining  relations  to  the  mean  gradient  structure. 
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CHAPTER  1 


INTRODUCTION 

For  many  reasons,  there  is  considerable  interest  in  the  vertical 

2 

profile  of  the  refractive  index  structure  parameter,  C^  .  This 
parameter  quantifies  effects  of  refractive  fluctuations 

important  for  electromagnetic  radiation  systems  whose  signals 
propagate  through  the  turbulent  atmosphere.  Examples  of  such  systems 
are  satellite  communications,  remote  sensing  such  as  clear  air  radars, 
astronomical  observation,  and  certain  weapons.  Xu  situ  measurements 
are  a  necessary  first  step  to  set  up  a  database  for  verifying  physical 
and  numerical  models.  Both  are  required  to  facilitate  the 
interpretation  of  remote  sensing  measurements  and  to  achieve  the  goal 
of  prediction  based  on  readily  available  atmospheric  data  and  models. 

This  project  was  an  analysis  of  the  first  segment  of  the  EWAK 
(E  xperiment  W  ithout  ronym)  electro-optical/meteorological 
experiment  at  Penn  State  University.  The  purpose  of  EWAK,  conducted 
principally  in  April  and  May  of  1986,  was  to  simultaneously  operate  at 
the  same  location  several  different  instruments  to  measure  the 
refractive  index  structure  parameter,  both  in  situ  and  remotely, 
throughout  different  synoptic  and  diurnal  conditions. 

Specifically,  this  thesis  describes  a  detailed  analysis  of 

research  aircraft  measurements  of  temperature  and  velocity  turbulence 

taken  with  the  instrumented  Beechcraft  Baron  of  the  Airborne  Research 

Associates  of  Weston,  Massachusetts.  From  these  two  measurements  can 

2 

be  derived  the  temperature  and  velocity  structure  parameters,  C.^ 

2  2  2 
and  C  Optical  C  is  a  function  of  C-,  . 
u  n  I 


The  turbulent 
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kinetic  energy  dissipation  rate,  e,  is  related  to  Cu2 .  Another 

important  relationship  examined  was  that  between  C^2  and  C  2 

T  u  ' 
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CHAPTER  2 

BACKGROUND  ON  ATMOSPHERIC  PROPAGATION 

2 . 1  Refractive  Index 

A  turbulent  fluid  creates  three  general  effects:  it  imposes 
varying  forces  on  any  objects  embedded  in  the  fluid  or  in  the  fluid 
path,  it  generates  fluxes  of  fluid  properties  (e.g.,  temperature  or 
momentum) ,  and  it  creates  (by  density  fluxes)  variations  in  the 
refractive  index  of  the  fluid  (Panofsky  and  Dutton  1984).  These 
refractive  index  variations  will  affect  electromagnetic  (EM)  and  sonic 
radiation  propagating  through  the  fluid.  The  atmosphere,  being  a 
turbulent  fluid,  has  a  temporally  and  spatially  variable  refractive 
index  structure.  The  atmosphere’s  turbulent  refractive  index 
fluctuations  have  plagued  astronomers  since  some  of  the  earliest 
recorded  observations.  The  twinkling  of  stars  is  one  of  the  most 
obvious  visible  manifestations.  Today,  our  atmospheric  uses  of  EM 
radiation  are  far  more  diverse  than  merely  receiving  visible  starlight 
with  our  eyes  or  at  a  telescopic  aperture.  Hence,  the  need  exists  to 
characterize  and  predict  the  refractive  index  and  its  variations  over 
arbitrary  pathlengths  and  view  angles. 

Henceforth,  in  this  thesis,  unless  it  is  specifically  otherwise 
stated,  the  discussion  will  be  concerned  with  propagating 
electromagnetic  radiation.  Refraction  can  be  considered  as  a  type  of 
scattering  (i.e.,  deviation  from  the  forward  direction).  Turbulent, 
non-homogeneous  variations  of  the  atmospheric  refractive  index  are  on 
the  order  of  one  part  in  10^.  This  results  in  very  small  scattering 
angles,  for  which  depolarization  and  intensity  attenuation  generally 
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can  be  neglected  (Strohbehn  1978) .  The  real  part  of  the  complex 
refractive  index,  n,  of  a  medium  can  be  thought  of  as  a  measure  of  the 
medium's  effect  on  wave  propagation  as  compared  to  propagation  in  a 
vacuum:  n(medium)-c(vacuum)/c(medium) ,  where  c  is  the  phase  velocity 
of  the  wave.  Refractivity ,  usually  denoted  by  N,  is  (n-1)  X  106. 

The  phase  speed  of  the  wave  is  altered  (slowed  in  air  compared  to 
vacuum,  thus  n>l) ,  but  the  energy  E  and  frequency  i/  remain  the  same, 
E-hi/  (Hecht  1986).  Using  geometric  optics  theory,  this  can  be  shown  to 
result  in  a  deviation  of  the  beam's  path  which  is  referred  to  as  the 
phenomenon  of  refraction  (in  molecular  dipole  theory,  it  is  referred 
to  as  a  scattering  event) . 

Two  factors  determine  the  phase  speed  (and  thus  refractive  index) 
in  a  medium:  density  (number  of  atoms,  molecules  or  particles  per  unit 
volume)  and,  since  each  molecule  has  a  wavelength  dependent 
polarizability,  molecular  composition  (Hecht  1986).  At  optical  or 
shorter  wavelengths,  temperature  variations  effectively  dominate 
refractivity  changes  since  atmospheric  molecules  have  low  rotational 
polarizability  at  these  wavelengths  (Balsley  and  Gage  1980).  However, 
at  longer  wavelengths,  specifically  microwave  radar  wavelengths,  both 
temperature  and  the  relative  concentration  of  various  molecular 
species  with  permanent  dipole  moments  determine  the  refractive  index 
gradients  affecting  radiation  transmittance  (Wesely  1976;  Balsley  and 
Gage  1980).  In  the  boundary  layer,  water  vapor  content  is  often 
highly  variable,  and  thus,  may  dominate  refractivity  changes  (Balsley 
and  Gage  1980) . 

If  temperature  and  molecular  composition  varied  negligibly  from 
the  standard  atmosphere,  then  the  refractive  index  profile  could 
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easily  be  calculated.  The  importance  of  turbulence  is  its  resultant 
mixing  (Strohbehn  1978).  Mixing  due  to  turbulence  causes  the 
variations  in  space  and  time  of  the  refractive  index  profile.  It  is 
the  random  nature  of  turbulence  that  introduces  chaos  into  the  density 
profile . 

2.2  Propagation  Effects 

Variations  in  refractive  index  can  cause  at  least  six  types  of  EM 
propagation  degradation:  beam  steering,  image  dancing,  beam  spread, 
spatial  coherence  degradation,  temporal  coherence  degradation,  and 
scintillation  (Dewan  1980) . 

Beam  steering  is  deviation  from  the  line-of -sight  (i.e.  ,  the  most 
direct  path  from  source  to  receiver) .  The  transverse  deviation  of  a 
narrow  beam  increases  with  subsequent  distance  from  the  deviating 
medium.  If  a  turbulent  eddy  is  larger  than  the  width  of  the  beam,  it 
may  act  as  a  "prism."  This  effect  can  be  quite  significant  for 
earth- to-space  uplinks,  because  the  planetary  boundary  layer  (PBL)  is 
highly  turbulent  (Dewan  1980) . 

Image  dancing  is  the  modulation  of  the  arrival  angle  of  the  whole 
beam  wavefront.  It  can  cause  the  focal  plane  for  the  beam  to  be 
angled  to  a  receiver  normal  to  the  beam,  causing  photographic  blurring 
or  necessitating  larger  apertures  with  lower  signal -to -noise  ratios 
(Dewan  1980) . 

Beam  spreading  is  caused  by  scattering  at  small  angles  near  the 
edge  of  the  beam.  This  occurs  when  the  beamwidth  is  greater  than  the 
turbulent  eddies.  This  phenomenon  occurs  most  on  the  space- to-earth 
downlink  (Dewan  1980) . 
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At  large  distances,  a  spherical  wave  from  a  point  source  can  be 
approximated  as  a  plane  wave.  According  to  Huygen's  model,  a  single 
wavefront  can  be  thought  to  consist  of  an  infinite  number  of  points 
each  emitting  a  "wavelet."  The  wavelets  add  together  to  give  the 
total  wave.  If  the  wavelets  have  a  constant  phase  relation,  the 
wavefront  is  said  to  be  coherent.  When  phase  difference  randomly  and 
continuously  changes,  the  wavefront  becomes  incoherent  (Hecht  1986). 
When  a  pathlength  difference  is  introduced  across  a  coherent  wavefront 
(e.g.,  refraction),  a  Fraunhoffer  diffraction  pattern  may  be  seen  in 
which  maxima/minima  occur  where  the  amplitudes  of  wavelets 
contructively/destructively  interfere.  Thus  diffraction  is  one  method 
of  testing  coherence  (Hecht  1986) . 

A  coherent  source  can  be  viewed  as  emitting  "pulses , "  or  wave 
packets  (this  model  combines  wave  and  particle  concepts) .  No  source 
is  perfectly  monochromatic,  thus  each  wave  pulse  has  a  frequency 
bandwidth.  The  inverse  of  the  pulse  frequency  band  is  the  "period" 
band  or  coherence  period  (Hecht  1986).  Temporal  (or  longitudinal) 
coherence  occurs  within  the  coherence  period.  Spatial  (or  horizontal) 
coherence  occurs  across  a  wave  pulse  emitted  from  a  coherent  source. 
Spatial  degradation  is  loss  of  phase  coherence  across  a  wavefront, 
causing  spot  blurring  on  photographic  images.  It  also  distorts 
photomixing  processes  which  rely  on  mixing  of  an  incoming  to  a  local 
signal  to  form  a  "beat."  Temporal  degradation  is  a  loss  of  phase 
coherence  during  the  coherence  period,  thus  distorting  amplitude 
modulated  signals  (Dewan  1980) . 

Scintillation  is  a  coherence  effect  describing  amplitude  changes 
at  a  focal  plane  (interference  patterns).  Spatial  scintillation  being 
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fluctuations  across  a  beamwidth  (or  wavefront) ,  temporal  scintillation 
being  fluctuations  in  time  of  the  total  beam  amplitude  (signal  fading) 
(Dewan  1980) . 

2 

2.3  Propagation  Parameters  and  Cn  Measurements 


A  variety  of  theories  have  been  developed  for  these  degradation 

effects.  In  general,  an  atmospheric  turbulence  parameter,  the 

2 

refractive  index  structure  constant,  Cn  ,  must  be  specified  along 

the  propagation  path.  By  measuring  some  relevant  propagation 

parameter  (e.g.,  scintillation)  the  theory  can  be  "inverted"  to 
2 

determine  C 

n 

One  such  instrument,  the  NOAA  (National  Oceanic  and  Atmospheric 

2 

Administration)  scintillometer,  measures  Cn  -optical  remotely, 

based  on  starlight  amplitude  fluctuations.  As  one  of  the  six 

distortions  described  previously,  scintillation  has  been  quantified 

2 

using  an  integral  of  Cn  (z),  where  z  is  the  altitude.  Tatarskii 
has  shown  for  temporal  scintillations,  time  averaged  amplitude 
fluctuations  (Dewan  1980)  ,  that 


<( ln(A/A  )]2>  -  0.56  k7/6 


C2  (z)  z5/6dz 


(2.1) 


where  A  is  the  signal  amplitude,  Aq  is  the  mean  amplitude,  k  is  the 

optical  wave  number,  and  the  angle  brackets  denote  a  time  average.  A 

2 

stellar  scintillometer  measures  a ^  (d^),  the  normalized  variance 

of  the  spatially  filtered  signal  irradiance  for  3  spatial  wavelengths. 

2 

d^.  Then  a  profile  of  Cn  is  calculated  for  seven  predetermined 
levels,  based  on  seven  weighting  functions  (only  four  of  which  are 


independent) 


C-2 1 


C2  -  B  2  R.  K.8/3  a2  (d. ) 
n  i  i  li  i  l 

where  R^  weight  is  given  to  the  ith  spatial  wavenumber 
(K^-27r/d^)  and 

r« 

■  l 


(2.2) 


B  -  1.87  X  10*13  [ 


Wc(z)  dz  ] 


(2.3) 


where  Wc(z)  is  a  composite  path  weighting  function  (Ochs  et  al. 

1976).  The  seven  predetermined  levels  are  line-of -sight  distance  from 
the  instrument,  so  that  actual  height  above  ground  depends  on  the 
zenith  angle  of  the  observed  star.  Some  improvements  to  the  NOAA 
scintillometer  were  made  by  AFGL  (Air  Force  Geophysical  Lab)  to  their 
instrument.  They  are  described  in  more  detail  by  Murphy  and  Battles 


(1986)  . 


An  NPS  (Naval  Postgraduate  School)  rQ  scintillometer  optically 
measures  the  transverse  coherence  length,  rQ ,  of  the  image  in  the 
focal  plane  (Donald  Walters,  telephone  conversation,  October  1987). 
rQ  is  a  scalar  measure  (in  units  of  aperture  size)  of  spatial 
wavefront  coherence.  Crudely  stated,  the  resolution  of  an  image 
cannot  be  improved  by  increasing  a  single  aperture  beyond  rQ  (Dewan 
1980).  This  helps  delineate  the  aperture  size  for  adaptive  optics 
elements.  It  has  been  shown  that 


r  -  18.5  A1'2  ( 
o 


C2(z)  dz  )'3/5 


(2.4) 


Jo 

2  -  2/3 

with  rQ  in  cm,  wavelength  (A)  in  m,  Cn  in  m  ,  and  altitude 
in  m.  Using  this  integral,  rQ  can  also  be  calculated  from  profiles 
of  Cn2  (Fried  and  Mevers  1974). 
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An  AFWL  (Air  Force  Weapons  Lab)  isoplanometer  optically  measures 

SQ,  the  isoplanatic  angle.  $  is  a  measure  of  angular  coherence 

for  an  object  with  angular  size  (vs.  a  point  source)  that  delineates 

the  region  through  which  the  object  can  be  viewed  coherently  without 

adaptive  optics  (Loos  and  Hogge  1979;  Murphy  and  Battles  1986). 

2 

can  a^so  be  calculated  from  a  vertical  profile  of  (Murphy  and 
Battles  1986) 


6  -  58000  A1'2  [ 

o 


z5/3  C2(z)  dz  ]'3/5 


(2.5) 


with  9q  in  /*rad,  and  A,  z  and  in  the  same  units  as  for 

V 


Notice  that  the  three  parameters  discussed  here  have  different 

propagation  path  weighting  functions  (Oj2  as  z3^6,  rQ  as  1.0, 

5 /3 

$Q  as  z  '  ).  Thus,  for  ground- to -space  optical  paths,  rQ  tends 
to  be  dominated  by  the  boundary  layer,  SQ  by  the  stratosphere,  and 
(7j2  is  somewhere  in  between. 
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CHAPTER  3 

MICROTURBULENCE  THEORY 

3.1  The  Refractive  Index  Structure  Function 

The  atmosphere  is  predominantly  hydrostatically  stable  above 
the  planetary  boundary  layer  (PBL) .  Nevertheless,  observations,  such 
as  contrails  and  rocket  trails,  have  clearly  shown  that  discreet 
turbulent  layers  exist  (VanZandt  et  al.  1981).  Although  turbulence 
and  byproducts  of  turbulence  (fossil  eddies)  are  undoubtedly  present, 
they  exist  only  on  an  intermittent  basis. 

While  atmospheric  turbulence  is  inherently  anisotropic,  at  small 
size  scales,  where  viscous  dissipation  occurs,  turbulence  can  be 
considered  isotropic.  For  isotropic  and  homogeneous  turbulence,  the 
structure  function,  D,  depends  on  spatial  scalar  distance,  £,  and  is 
defined  as  the  mean  squared  difference  of  the  measured  passive 
parameter  (e.g.,  temperature)  at  two  points  in  space  (Hinze  19/ j). 

DT(f)  -  <(T(r)  -  T(r+£) )2>  (3.1) 

where  the  angle  brackets  indicate  an  averaged  quantity.  In 
anisotropic  turbulence,  D  also  depends  on  the  orientation  of  the 
separation  (Hinze  1975). 

The  autocorrelation  R  is  defined  as  the  average  of  the  product  of 
deviations  from  the  mean  (e.g.,  T')  at  two  points  in  space. 

R^f)  -  <T '  ( r )  T'(r+£)>  (3.2) 

In  the  limit,  as  £  approaches  infinity,  the  structure  function 
approaches  twice  the  variance.  Variance  is  equal  to  the 
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autocorrelation  function  at  £-0  (Tatarskii  1971).  Thus, 

2  0T(O  -  R-pCO)  -  R^O  (3.3) 

In  the  inertial  subrange,  it  can  be  shown  that  the  structure 
function  has  a  power  law  relation  to  separation  with  a  constant 
designated  C ^  (e.g.,  C^2)  (Hinze  1975). 

Dt<£)  -  C2  £2/3  (3.4) 

Similarly,  we  can  define  for  refractive  index, 

Dn(0  -  C2  £2/3  (3.5) 

The  structure  constant  is  only  a  "constant"  in  that  it  is  independent 

of  the  separation  distance;  _t  is  still  variable  over  time  and  space 
o 

[e.g.,  Cj  (x,y,z,t)].  The  optical  refractive  index  structure 
constant  can  be  related  to  the  temperature  and  humidity  structure 
constants  (Nastrom  et  al .  1981). 

3  -  *!  A*  2  £1  ^qT  *  A  4 

2 

where  f -  -  0.60  P/T  ,  and 

f2-  -  1.2  Pq/T2  6  -  77.6E-6  P/  T  6  (3.7; 

with  pressure  in  mb , and  where  q  is  specific  humidity.  Neglecting  the 
effects  of  humidity,  this  reduces  to  (Tatarskii  1971) 

C2  (optical)  -  (79  X  10'6  (P/T2))2  C2  (3.8) 

2 

During  the  EWAK  experiment,  the  aircraft-measured  allowed 

2 

computation  of  "dry"  Cr  -optical.  A  Lyman-alpha  detector  on  the 
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2 

aircraft  was  intended  to  measure  C  but  the  threshold  sensitivity 
was  too  high  for  the  vapor  concentrations  at  the  altitudes  of 
interest . 

2 

3.2  C  From  Radar  Measurements 
n 

Data  from  a  clear-air  radar  (often  referred  to  as  a  profiler)  can 

2 

be  used  to  calculate  Cn  -radar  based  on  radar  reflectivity.  The 

backscatter  cross  section  from  received  power,  a,  at  the  wavelength, 

2 

A,  can  be  characterized  by  (VanZandt  et  al.  1977)  as 

a  -  0.38  C2  A*1/3  (3.9) 

2 

Tatarskii  (1971)  derives  a  dimensional  relation  for  C  based 

n 

on  a  refractivity  gradient  across  a  turbulent  layer 

C2  -  a  M2  L4/3  (3.10) 

where  a  is  an  empirical  constant,  L  is  the  layer  thickness  and 

M2  -  -77.6  x  10’6  (P/T)  [  N/g( 1  +  15500q/T)  -  7750  q'/T] 

(3.11) 

where  P  is  pressure  in  mb,  T  is  absolute  temperature,  q  is  specific 

humidity  in  kg/kg,  q'  is  3q/3z,  g  is  acceleration  due  to  gravity,  and 

2 

N  is  the  Brunt-Vaisala  frequency  squared.  For  Cn  -optical,  q  and 

2 

q'  are  set  equal  to  zero  in  M  (Dewan  1980),  which  gives  the 
resulting  ratio  relating  optical  and  radar  values 

(C2-  radar)  /  (C2-  optical)  -  X2 


(3.12) 
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with  (Fairall  and  Thomson  1985) 

X2  -  [1-  7750  q'  (q/T) (g/N)  +  15500  q/T]2  (3.13) 

2  2 

3.3  Cj,  and  Cu  Relationships 


Dimensional  analysis  shows  that  in  the  inertial  subrange  of 
isotropic  turbulence,  the  variance  spectral  density  obeys  a  k*3//3 
wavenumber  dependence . 

^T(kx)  -  0T  *T  e‘1/3  k'5/3  (3.14) 


where  0j  is  an  empirical  constant,  x-p  is  the  rate  of  dissipation 
of  one  half  the  temperature  variance,  e  is  the  dissipation  rate  of  the 
turbulent  kinetic  energy,  and  k^  is  the  spatial  wavenumber  (k^  - 
2ir  A*^)  (Panofsky  and  Dutton  1984).  Since  the  one  dimensional 
variance  spectrum  is  defined  as  the  Fourier  transform  of  the 
autocorrelation  (Panofsky  and  Dutton  1984) , 


<^T(kl)  “  l 


RT(f1)  cos  (k^p)  d£ 


(3.15) 


it  can  be  shown  using  the  relation  of  the  structure  function  to 
variance  and  the  autocorrelation,  that 


^T(kp)  -  0.25  C2  k‘5/3 


(3.16) 


in  which  0.25  is  the  result  of  combining  several  constants.  Mow,  the 
so-called  Corrsin  relation  is  obtained  from  3.14  and  3.16  (Tatarskii 
1971)  . 

<4  -  4  0T  xT  «'1/3  (3.17) 


0  has  been  determined  for  temperature,  0^-0 . 8 ,  and  velocity, 
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0U~O . 5  (Champagne  et  al.  1977;  Panofsky  and  Dutton  1984).  is 

assumed  to  be  equal  to  Pj.  Since  Xu“€  for  velocity,  a  simplified 
Corrsin  relation  is  obtained 

C2  -  4  p  £2/3  -  2  £2/3  (3.18) 

u  u 

Using  this  relation,  it  is  simple  to  calculate  e,  the  turbulent 
kinetic  energy  dissipation  rate. 

e  is  an  important  term  in  the  variance  budget  equations  for  a 
turbulent  layer.  These  equations,  neglecting  transport  (Ottersten 
1969)  ,  are  (Fairall  and  Markson  1984) 


de/dt  +  u'w'  du/dz  -  g /8  v'9'  -  -«  (3.19) 

2  do2/dt  +  w’T'  39/9 z  -  *XT  “  -cx  «1/3 /  ^T)  (3.20) 


where  primed  quantities  are  fluctuations,  barred  quantities  are  means. 
2  2  2 

e-l/2(ou  +<xv  +crw  )  is  the  turbulent  kinetic  energy,  g  is 
the  acceleration  due  to  gravity,  z  is  altitude,  o g  is  variance  of 
potential  temperature,  Xj  is  the  rate  of  dissipation  of  one  half  the 
thermal  variance,  and  u,  v,  and  w  are  the  horizontal  and  vertical 
velocity  components  respectively.  If  the  flux  Richardson  number  is 
defined  as 


Rf  -  [ (V  8 '  g/9  )/(  u'w'  3u/3z) ] 


(3.21) 


and  the  eddy  diffusion  coefficient  is  invoked 


w'  9  '  -  -Ky  <59 /dz  (3.22) 


u'w'  —  -K  3u/3z 
m 


(3.23) 
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chen  the  following  ratio  results  when  using  the  Corrsin  relations  in 
a  ratio  of  equation  3.20  to  3.19.  (Fairall  and  Markson  1984) 


C2  /  C2  -  1.6  (9/g)  39/dz  Rf/(1  -  Rf )  (3.24) 

C2  /  C2  -  1.6  (9/e)  89/dz  V(Pr*  Ri}  (3'25) 

where  R^.  is  replaced  by  the  gradient  Richardson  number, 

Rt  -  g  31n0/3z  /  (3u/3z) 2  (3.26) 

and  the  turbulent  Prandtl  number, 

Pr-  Km/KH  (3.27) 


In  an  actively  turbulent  layer  in  the  free  troposphere , 

Rj^  =*  Ricrit“0-25  and  Pr»l  (Dewan  1980),  then 

C2  /  C2  «  1.6  (9/e)  39/dz  (3.28) 

Gradient  Richardson  number  is  the  basis  for  a  model  developed  by 

2 

VanZandt  et  al.  (1978,  1981).  This  model  predicts  -optical  and 

2 

c  based  on  measured  values  of  N  (-g /9  39/d z)  and  S  (-  du/dz  + 

2 

3v/3z  ),  and  input  values  of  L,  the  outer  scale  or  overturning  size 

of  the  turbulent  eddies.  Note  that  the  gradient  Richardson  number  is 

equal  to  N/S .  Based  on  equation  3.28,  this  model  can  also  be  related 

2  2 

to  the  ratio  of  and  Cu  .  The  model  is  based  on 
Tatarskii's  dimensional  relationships 

C2  -  a  M2  LV3  (3.29) 

n 


e  -  b  S3/2  L2 


(3.30) 
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where  a  and  b  are  constants.  A  statistical  integration  is  performed 

assuming  that  turbulence  occurs  in  discreet  layers  with  R^<  0.25. 

Theoretical  probability  density  functions  are  assumed  for  shear  and 

temperature  gradients.  The  internal  model  variable  L  is  adjusted  to 

2 

give  best  match  to  radar-derived  Cn  profiles.  The  model  is 
inherently  based  on  averaging  over  a  length  much  greater  than  the 
turbulent  scale,  of  which  only  a  fraction,  F,  is  actively  turbulent 
(VanZandt  et  al.  1981).  Thus  high  resolution  in  situ  measurements 
(such  as  measurements  by  aircraft)  may  differ  from  model  predictions. 

3.4  The  Turbulent  Kinetic  Energy  Dissipation  Rate 

2 

Few  in  situ  measurements  of  or  e  have  been  made  above  the 

2  2 

PBL.  Barat  and  Bertin  (1984)  measured  CT  and  Cu  with 

stratospheric  balloons,  and  Fairall  and  Markson  (1984)  with  an 

aircraft.  Kennedy  and  Shapiro  (1980)  calculated  various  turbulence 

parameters  from  research  aircraft  observations  of  CAT  episodes 

associated  with  jetstreams. 

2 

Substituting  for  Cu  in  equation  3.22  gives  another 
relationship  for  « 

C2  -  3.2  (0/g)  dd/dz  R£  «2/3  (3.31) 

2 

This  can  be  used  to  calculate  e  from  measurements,  and  from 

high  altitude  (where  water  vapor  contribution  is  negligible)  profiler 

2 

measurements  of  Cn  -optical  (Tatarskii  1971)  with 

C2  -  [79  X  10'6  P/T2)]2  3.2  (d/g)  dd/dz  R.^  f2/3 


(3.32) 
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or  by  assuming  a  Corrsin  relation  for  C  2 


C2  -  3  2  v  f'1/3 
Cn  *n  £ 


and  assuming 


*n  “  Kn  <a"/az)2 
KH  “  Ri  N"1  € 

Kn  "  ^ 

M2  -  (3n/3z)2 


then 


(3.33) 


(3.34) 

(3.35) 

(3.36) 

(3.37) 


C2  -  3.2  R.  M2  N*1  c2/3  (3.38) 

n  i 

where  N  is  Brunt-Vaisala  frequency  squared  and  M  is  Tatarskii's 
refractivity  gradient  (Barat  and  Bertin  1984,  Panofsky  and  Dutton 
1984).  Notice  that  equation  3.28  implies  the  factor  should 
actually  be  R^/(Pr  *  R^)- 


3.5  Turbulence  Size  Scales 


For  oceanic  turbulence,  Gregg  (1987)  defines  a  turbulence 

3/2  1/2 

bandwidth  based  on  the  ratio  of  Ozraidov  scale  Lq-[«/N  ]  (a 

buoyancy  scale  at  which  buoyancy  equals  inertial  force)  and  the 
Kolmogorov  scale,  ]  /  (the  inner  scale  or  cutoff 

wavelength  for  the  inertial  subrange  where  viscous  equal  inertial 
forces) . 

L0/n  - 


(3.39) 
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where  1/  is  the  kinematic  viscosity.  The  eddy  overturning  scale  is 
proportional  to  the  buoyancy  scale  in  actively  turbulent  layers. 
Thus,  the  inertial  subrange  occurs  for  size  scales  less  than  the 
buoyancy  scale  and  greater  than  the  inner  scale.  Gregg  provides  the 
following  empirical  guideline  to  interpret  the  bandwidth  ratio 
parameter  as  a  turbulence  activity  parameter. 


Value  of  e/i/N 
<15 
>200 
>10000 


Turbulence  State 

Decaying  turbulence,  fluxes  negligible 
Fully  isotropic 
Fully  developed  flow 


A  spectrum  begins  to  deviate  from  -5/3  slope  at  size  scales  one  order 
of  magnitude  greater  that  the  Kolmogorov  microscale  (Hill  and  Clifford 
1978).  For  Lo=»10rj ,  no  inertial  subrange  should  be  apparent. 
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CHAPTER  4 

THE  EWAK  EXPERIMENT 


4 . 1  Overview 


Two  segments  of  the  EWAK  experiment  were  analyzed  in  this  thesis 
study.  The  first  and  major  segment  was  performed  from  14  April  1986 
to  6  May  1986.  This  segment  was  conducted  at  the  Meteorological  Field 
Site,  locate J  several  miles  southwest  of  State  College,  PA,  on  the 
Rock  Springs  Agronomy  Research  Farm.  The  relevant  instruments 
operating  at  various  times  during  the  experiment  were  an  Air  Force 
Geophysical  Lab  (AFGL)  scintillometer,  Rome  Air  Development  Center 
(RADC)  scintillometer,  an  AFGL  thermosonde,  an  Air  Force  Weapons  Lab 
(AFWL)  isoplanometer,  a  Naval  Postgraduate  School  (NPS)  rQ 
scintillometer,  a  Penn  State  University  (PSU)  Doppler  radar  profiler, 
a  PSU  Doppler  sodar,  and  the  Airborne  Research  Associates  (ARA) 
aircraft.  A  summary  of  the  database  relevant  to  this  study  is  given 
in  Table  1. 

The  aircraft  was  based  at  University  Park  Airport,  several  miles 
northwest  of  State  College,  PA.  Data  gathering  flights  were  made  over 
the  field  site  and  insofar  as  possible  within  a  radius  of  10  miles. 
Eight  flights  were  made,  six  yielded  useable  results:  flight  1  was  a 
test  flight  and  flight  5  had  a  total  failure  of  the  instrument 
recorder.  Each  flight  consisted  of  a  fairly  rapid  climb  to  maximum 
altitude  (approximately  30  kft  or  10  km)  and  then  a  slow  descent  (500 
fpm) .  All  descents  were  analyzed  (because  of  the  greater  resolution) , 
but  for  comparison  purposes  several  ascents  were  also  analyzed.  Level 
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flight  data  of  return  trips  by  the  aircraft  to  Boston  was  also  added 
to  several  flights . 

The  second  segment  of  EWAK  was  conducted  at  Griffiss  Air  Force 
Base,  Rome,  NY,  from  11  August  through  15  August  1986.  Both  ARA 
aircraft  and  RADC  scintillometer  measurements  were  taken,  however, 
scintillometer  data  taken  during  these  flights  was  unavailable  for 
this  thesis.  One  additional  flight  was  made  on  20  October  1986.  For 
these  flights,  the  aircraft  was  based  at  Oneida  County  Airport  and 
performed  similar  flight  patterns  centered  around  RADC's  field  site  in 
Verona,  NY. 

4.2  Aircraft  Instrumentation 

A  summary  of  the  instrumentation  available  on  the  aircraft  is 
given  by  Fairall  and  Markson  (1984;  1987).  Thus  only  the  relevant 
turbulence  measurements  are  discussed  here.  Temperature  and  velocity 
signals  were  amplified,  filtered,  and  recorded  on  FM  tape.  Recorder 
gain  settings  were  correlated  to  tape  counters  and  recorded  manually. 
Filters  were  set  for  a  bandpass  between  2  and  200  Hz  for  the 
temperature  signal,  and  5  and  200  Hz  for  the  velocity  signal.  The  FM 
cape  recorder  was  set  at  the  15/16  ips  speed,  which  provided  an 
additional  low  pass  filtering  at  256  Hz. 

Temperature  fluctuations  were  measured  with  a  fast  response 
resistance  bridge  driving  microthermal  sensor  probes  as  the 
temperature  sensitive  resistance  elements.  The  bridge  used  was  a 
Thermo  Systems,  Inc.  (TSI)  type  1044,  dc  Wheatstone  bridge  with  a 
frequency  response  of  800  hz.  The  microthermal  probes  were  TSI  type 
1210  with  T1 . 5  configuration  (W  wire,  U.5fim  dia.,  ice  point  resistance 
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Ri-50) .  A  wingtip  boom  carried  matched  sensors  that  were  connected 
to  the  two  bridge  arms.  One  sensor  was  covered  so  that  bridge  signal 
output  was  proportional  to  temperature  fluctuations  sensed  by  the 
exposed  probe. 

Velocity  fluctuations  were  measured  with  hot-wire  anemometry  in  a 
constant  temperature  circuit.  The  bridge  used  was  a  TSI  model  1050. 
Probes  used  were  the  same  TSI  type  1210  T1.5  tungsten  sensors  with  a 
50%  overheat. 

Other  measurements  taken  that  were  used  in  this  analysis  were  air 
temperature  (measured  with  a  standard  Rosemont  sensor  and  bridge 
circuit),  dew-point  temperature  (measured  with  a  Cambridge  Systems, 
Inc.  aircraft  chilled  mirror  instrument),  pressure  altitude  (measured 
with  a  Rosemont  sensor) ,  and  radar  altitude  (measured  with  a  Bonzair 
Mark  10  altimeter) . 

4.3  Additional  Instrumentation 

Both  the  RADC  and  the  AFGL  scintillometers  are  instruments  that 
measure  stellar  scintillation  amplitude  fluctuations  based  on  the  NOAA 
model  II  instrument  described  in  detail  by  Ochs  et  al .  (1977). 
Additional  modifications  were  made  to  the  AFGL  scintillometer  as 
described  by  Murphy  and  Battles  (1986).  The  AFWL  isoplanometer  and 
NPS  rQ  scintillometer  are  instruments  also  based  on  measuring 
stellar  scintillations  and  may  be  regarded  to  be  derivatives  of  the 
NOAA  instrument.  These  instruments  are  described  in  some  detail  by 
Eaton  et  al.  (1985),  and  Stevens  (1985),  respectively. 
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The  AFGL  thermosonde  is  essentially  a  special  temperature 

turbulence  measuring  package  attached  to  a  radiosonde.  The  instrument 

package  is  described  in  detail  by  Brown  et  ai.  (1982). 

Penn  State's  Doppler  radar  profiler  is,  basically,  of  the  type 

described  by  Strauch  et  al.  (1986 j.  It  is  a  VHF  (6  m  wavelength) 

pulsed  Doppler  radar  with  a  phased  array  antenna  located  near  McAlevys 

Fort,  PA,  on  Shantytown  Road,  approximately  10  miles  south  of  State 

Collage,  PA.  With  the  profiler,  vertical  profiles  of  wind  speed  and 

direction  are  obtained  as  well  as  radar  reflectivity  (backscatter 

cross  section) .  Radar  volume  reflectivity  can  be  used  to  calculate 
2 

Cn  values  at  each  of  the  radar  gates  (VanZandt  et  al.  1978).  The 
profiler  produces  both  high  and  low  resolution  profiles.  High 
resolution  is  -290  m  from  -1  km  to  8  km  MSL.  Low  resolution  is  -870  m 
from  -1.5  km  to  17.5  km  MSL. 

PSU  Doppler  sodar  data  was  not  used  in  this  study. 
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CHAPTER  5 
DATA  PROCESSING 


The  relation  of  the  structure  constant  to  the  variance  spectral 

density  as  a  function  of  wavenumber  (Equation  3.1b)  allows 

2  2 

determination  of  C^,  and  Cu  from  spectral  analysis  of 
turbulence  data.  The  aircraft  data  was  a  time  series  record  of 
voltage  fluctuations,  which  could  be  transformed  into  frequency  . 
spectra.  Taylor's  frozen  turbulence  hypothesis  allows  for  translation 
of  the  spatial  spectrum  <£(k^)  to  a  frequency  spectrum  <£'(f) 

(Panofsky  and  Dutton  1984) ,  by  the  relation 

kx  -  2jrf/u  (5.1) 

where  k^  is  the  spatial  wavenumber,  u  is  the  relative  speed  of  the 
passing  eddies  (in  this  case,  u-60  m/s,  the  aircraft  speed),  and  f  is 
the  frequency  in  Hz.  Since  total  spectral  intensity  must  be  invariant 
(Panofsky  and  Dutton  1984),  temporal  and  spatial  spectra  can  be 
related  by 

^(k1)dk1  -  *'(f)df  (5.2) 

Substitution  gives 

*'(f)  -  ^(2wf/u)dk1  /df  -  <(>( 2*f/u)  2w/u  (5.3) 

and  then 


^(f)  -  2jt/u  ^T(kL) 

(5.4) 

^(f)  -  2jt/u  (0.25  Cj  k‘5/3  ) 

(5.5) 
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^(f)  -  (2*/u)'2/3  (0.25  C2  f'5/3) 


(5.6) 


Similarly  for  velocity 

^(f)  -  (2*/u)'2/3  (0.25  C2  f‘5/3) 


(5.7) 


Rearranging  gives  the  equation (s)  relating  the  structure 
functions  to  the  measured  spectra 


CT,u  “  4  <2*/u)2/3  *x,u(f)  fV3 


(5.8) 


where  <fr' ^  u(f)  is  a  representative  value  in  the  inertial  subrange 
(Fairall  and  Markson  1984) . 

Another  method  to  determine  the  structure  constants  (the  single 
probe  RMS  method)  is  suggested  by  the  relation  of  total  variance  to 
spectral  density  (Fairall  and  Markson  1984) 


Variance  -  mean  square 


*'(f>  df 


(5.9) 


where  fu~upper  frequency  limit,  and  f^-lower  frequency  limit. 

Thus  the  filtered  RMS  of  the  signal  can  be  used  to  directly  obtain  the 
structure  constant  by  integrating  5.6  or  5.7 


Variance  -  (RMS)2-  (2*/u)  '2/30 . 25  C2  u(-  \  )(fu'2/3-  t{2/Z) 

(5.10) 


Rearranging  gives 


CT,u  "  3  <2*/u)2/3  (^S)2  <f{2/3  -  ^2/3  yl  (5.11) 

2 

The  thermosonde  instrument  calculated  C^,  directly  using  the 
double  probe  RMS  method.  If  two  probes  with  a  known  separation  £  are 
used,  then  with  an  RMS  of  the  difference  of  the  signal 
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(Brown  et  al.  1982) 

C2  -  <  (T1  ‘  T2)2  >  /  ?2/3  “  (MS)2/  £2/3  (5.12) 

The  spectral  method  has  the  advantage  of  greater  accuracy  over 
the  RMS  method  because  a  region  of  the  spectrum  can  be  selected  which 
is  relatively  free  of  noise.  The  RMS  variance  includes  all  noise  in 
the  bandwidth  analyzed  (Fairall  and  Markson  1984).  However,  the 
disadvantage  to  the  spectral  method  is  the  greatly  increased 
computation  time.  RMS  methods  can  be  done  on  a  virtually  real  time 
basis . 

The  aircraft  temperature  turbulence  data  was  taken  as  voltage 
fluctuations  in  bridge  output  due  to  resistance  fluctuations  of  the 
probes.  Temperature  fluctuations  cause  these  resistance  changes. 

This  relation  is  given  by 

dV/dt  -  i  Gfi  dR/dt  (5.13) 

where  the  current  to  the  probe  sensor  is  i,  and  bridge  amplification 
or  gain  is  Gg(«4000)  (Fairall  and  Markson  1984) . 

The  resistance- temperature  relation  is  linear  for  tungsten 
sensors  over  a  broad  temperature  range,  giving 

R  -  R.[l  +  a(T  -  Ti)]  (5.14) 

where  and  T^  are  values  at  the  ice  point  (Ri-5fl) ,  and  a  is 
the  effective  thermal  resistance  coefficient  (a-. 0026  K*3-),  which 
includes  the  loss  of  sensitivity  due  to  probe  support  effects  (Fairall 
and  Markson  1984).  Thus, 


dR/dt  -  Rl  a  dT/dt 


(5.15) 
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and 


dV/dt  -  i  Gg  G  (Rt  a  dT/dt)  (5.16) 

where  the  added  G  term  accounts  for  any  additional  signal  conditioning 
or  gains  applied  (usually  G-10) .  The  voltage  variance  spectrum  is  the 
Fourier  transform  of  the  square  of  voltage  fluctuations  which  then  can 
be  related  to  the  temperature  variance  spectrum  (Fairall  and  Markson 
1984) , 

^  -  (i  G  Gg  R1  a)2  ^(f)  (5.17) 

Including  noise  sources ,  this  would  be 

4$  -  (i  G  Gg  Rt  a)2  (^  +  Nu  +  NT  +  Nv/(iR.a)2  )  (5.18) 

where  Nv  is  broadband  voltage  noise  in  the  system,  is 
temperature  noise  detected  by  the  sensors  but  not  due  to  turbulence 
(e.g.,  sonic  adiabatic  compression),  and  Nu  is  velocity 
contamination  due  to  cooling  rate  variations  associated  with  velocity 
fluctuations.  An  indepth  study  of  signal- to -noise  ratio  was  done  by 
Fairall  and  Markson  (1984).  It  was  concluded  that,  in  flight,  most 
noise  was  temperature  noise.  Although  the  source  of  this  noise  could 
not  be  identified,  it  could  be  lessened  by  changing  sensor  location. 
Optimal  sensor  location  prov  -a  t.>  be  on  the  wingtip.  Sensor  currents 
were  kept  at  a  maximum  just  below  the  velocity  contamination 
threshhold  (i-5  mA) . 

Velocity  turbulence  data  was  also  taken  as  a  time  series  of 
voltage  fluctuations.  King's  Law  voltage  relationship  (Hinze  1975) 


was  used 


-  C-42  - 


V2  -  V2  +  B  u1/2  (5.19) 

o 

2 

where  V  and  B  are  calibration  constants  (Fairall  and  Schacher 
o 

2  2 

1977)  computed  for  50%  overheat  at  STP  as  Vq  -2.7  Volts  and 
2  -11/2 

Bq— 1 . 0  Volts  /(ms  )  '  .  Correcting  for  varying  temperature 
and  pressure  gives 

B  -  Bq  (P/1013)1/2  (288/T)1/4  [ (453 -T)/288 ]  (5.20) 

(Fairall  and  Markson  1984) .  Velocity  fluctuations  cause  bridge  outpu 
voltage  variations  as 

dV/dt-  [  B  G/(4  V  u1/2)]  du/dt  (5.21) 

with  G  additional  signal  conditioning  gains  (usually  G-100)  and  V  the 
mean  bridge  voltage  (about  3  Volts  for  nominal  aircraft  speeds) . 

Relating  the  square  of  fluctuations  to  the  variance  spectrum 
gives 

dy  -  (B  G/  4  V  7u)2  (5.22) 

2  2 

The  final  forms  used  for  calculating  and  from  the 
aircraft  data  (Fairall  and  Markson  1984)  were 

C2  -  4  (2jt/u) 2/3  (i  G  Gb  a)*2  ^(T)(f>  f5/3  (5.23) 

C2  -  4  (2*)2/3  u1/3  (4  V/  B  G)2  ^(u)(f)  fV3  (5.24) 

Flight  12  gain  settings  were  not  recorded.  They  were  assumed  to 
be  the  same  as  flight  11  for  processing,  but  for  that  flight,  this 
should  be  noted. 
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Initially,  the  temperature  and  velocity  channels  of  the  taped 
data  were  converted  from  analog  to  digital  with  an  Infotek  A/D  board. 
Approximately  two  minutes  of  analog  data  were  digitized.  A  128  point 
Fast  Fourier  Transform  (FFT)  was  done  and  a  final  spectrum  produced 
from  an  average  of  400  spectral  blocks.  This  gave  a  vertical 

resolution  of  0.3  to  0.5  km,  about  the  same  as  the  radar  profile. 

2 

The  log  of  power  spectral  density  (Volts  /Hz)  was  plotted  as 
2 

decibels  [10*log  (V  /Hz)]  against  log  of  frequency  (Hz)  for  4  to  256 
Hz.  A  key  frequency  was  chosen  that  consistently  appeared  in  the 
inertial  subrange  (in  a  region  of  high  signal  to  noise)  and  the 
spectral  density  at  this  point  was  used  to  plot  a  *5/3  slope  through 
che  spectrum  to  show  goodness-of-fit  to  the  Kolmogorov  turbulence 
theory.  Due  to  the  small  number  of  points,  the  spectra  obtained  with 
this  method  were  somewhat  crude  and  very  difficult  to  interpret. 

Figure  1  shows  three  different  spectra  from  flight  2  as  examples 
(indicated  as  EWAK2BT-rec#) .  Record  1  has  a  high  signal-to-noise 
ratio,  record  31  is  very  noisy  with  little  signal,  and  record  45  is  an 
unusual  steeply  sloped  spectra. 

A  second  method  gave  far  more  satisfactory  results.  Roughly  two 
minutes  of  data  were  analyzed  by  a  dual  channel  Hewlett-Packard  3562A 
dynamic  signal  analyzer  by  linearly  averaging  55  times  with  an  overlap 
of  75*.  A  1600  point  FFT,  using  a  Hanning  window,  produced  a  power 
spectrum  from  .125  to  100  Hz.  This  frequency  range  corresponds  to  a 
wavenumber  range  of  0.01  to  10  m*^.  Resolution  of  the  spectra  was 
greatly  increased.  The  same  three  records  from  flight  2  are  shown  in 
Figure  2 . 
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Figure  1.  Voltage  variance  spectra  from  first  FFT  program  (128 
points) .  For  each  record,  the  top  spectrum  (channel  »*1)  is 
temperature  variance  and  the  bottom  spectrum  (channel  #2)  is  velocity 
variance.  Voltage  power  spectral  density  (dB)  is  plotted  against 
frequency  (log  Hz).  The  slanted  line  represents  a  -5/3  slope  for 
comparison.  Coded  flight  name  (e.g.,  EWAK2BT)  and  record  number 
(e.g.,  -  1),  tape  footage,  and  altitude  are  shown  interior  to  each 
channel  #1  graph.  Three  different  dual-channel  records  in  flight  2 
are  shown:  1,  31,  and  45. 
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Figure  1. 


(continued) 
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Figure  1.  (continued) 
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Figure  2.  Voltage  variance  spectra  produced  by  the  signal  analyzer's 
1600  point  FFT.  For  each  record,  Ovl  (top)  is  temperature  variance 
and  Ov2  (bottom)  is  velocity  variance.  Voltage  power  spectral  density 
(dB)  is  plotted  against  frequency  (log  Hz).  Processing  parameters  are 
stated  above  each  spectrum.  Coordinates  of  the  cursor  are  given  in 
the  upper  left  corner  of  each  spectrum.  Vertical  axis  range  has  been 
fixed  as  stated  in  the  lower  left  of  each  spectrum.  Coded  flight  name 
(e.g.,  EAIR)  and  record  number  (e.g.,  -  1),  tape  footage,  and  altitude 
are  annotated  above  each  Ovl  graph.  The  same  three  dual -channel 
records  in  flight  2  are  shown  as  in  Figure  1. 
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CHAPTER  6 
ANALYSIS  METHODS 

6.1  Aircraft  Turbulence  Data  Analysis 

An  indepth  analysis  on  aircraft  instrumentation/experimertal 
error  was  doie  by  Fairall  and  Markson  (1984).  Uncertainty  due  to 
measurement  error  will  therefore  not  be  discussed  here.  However,  for 
this  thesis,  much  effort  was  put  into  detailed  analysis  of  the 
aircraft  turbulence  data  in  attempts  to  improve  analysis  precision. 
These  processes  are  described  in  this  chapter. 

To  quantify  goodness-of -fit,  a  linear  regression  analysis  was 
performed  on  the  log- log  spectral  data  and  percent  error  was 
calculated  for  deviation  from  -5/3  slope.  Most  spectra  exhibited  the 

-5/3  dependence  out  to  the  broadband  noise  level  of  the  instruments, 

-7  2  -6  2 

*3  X  10*  Km  for  temperature  and  «1  X  10*  (m/s)  m  for 

velocity  (Fairall  and  Markson  1985) .  Because  of  this  background  noise 

level,  notable  in  the  high  frequency  range  (over  20  Hz),  it  was 

determined  that  the  best  fit  would  be  obtained  at  the  lowest  frequency 

possible  within  the  inertial  subrange.  Unfortunately,  the  temperature 

data  had  been  low  pass  filtered  at  2  Hz  and  the  velocity  data  at  5  Hz. 

Filters  of  the  type  used  on  the  aircraft  (Thermo  Systems,  Inc.  model 

1057  signal  conditioner  with  TSI  model  1051-2  monitor  and  power 

supply)  were  analyzed  with  white  noise  input  on  the  signal  analyzer. 

Then  a  power  spectrum  ratio  of  unfiltered  noise  to  filtered  noise  was 

calculated  (Figures  3a  and  3b) .  This  function  could  then  be 

multiplied  point  by  point  to  essentially  "defilter"  the  data. 
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Spectral  drop-off  was  still  noted  in  the  low  frequency  end  «1 
Hz)  of  the  spectra  after  this  filter  function  was  applied.  Several 
reasons  for  this  drop-off  were  proposed:  a  real  phenomenon  (i.e.,  the 
outer  scale  of  the  inertial  subrange  or  the  production  size  of  the 
turbulence),  a  non-linear  characteristic  of  the  filter  when  the  signal 
is  log-distributed  and  not  linearly  distributed  like  the  white  noise, 
or  individual  error  in  the  filter  (two  identical  filters  were  tested 
and  found  to  have  different  signatures).  Assuming  f  is  1  Hz  and  u  is 
60  m/s,  the  frozen  turbulence  theory  (A*^-2jrf/u)  gives  A  of  10  m  as 
the  production  scale.  For  the  altitudes  flown,  this  is  within  an 
order  of  magnitude  of  typical  estimates  (VanZandt  et  al.  1981;  Gregg 
1987;  Barat  and  Bertin  1984).  For  comparison  to  the  characteristics 
of  the  measured  filter  function,  an  atmospheric  filter  function  was 
derived  to  compensate  for  non-linear  filter  effects  and  any  unknown 
atmospheric  effects  (e.g.,  variable  production  size  scale).  Several 
sample  spectra  from  the  boundary  layer  with  exceptionally  high 
signal- to-noise  ratio  and  -5/3  slope  were  used  to  derive  a  function 
based  on  -5/3  slope  (Figures  4a  and  4b). 

All  spectra  were  then  run  through  a  three  pass  progressive 
editing  program.  If  the  percent  error  on  a  regression  fit  from  5  to 
18  Hz  was  greater  than  20%,  the  spectrum  was  rejected  on  pass  1. 

Figure  5a  again  shows  records  1,  31,  and  45  of  flight  2  for  pass  1. 

On  the  second  pass,  the  atmospheric  filter  was  used  to  defilter  the 
data  from  0.1  to  5.6  Hz  (the  cutoff  frequency  where  both  2  and  5  Hz 
filters  were  approximately  equal  to  one).  Figure  5b  shows  the  same 
three  records  for  pass  2.  If  percent  error  was  greater  than  15%  on  a 
second  regression  fit  from  2  to  18  Hz,  then  the  spectrum  was  run 
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Figure  5a.  Voltage  variance  spectra  showing  editing  procedure.  These 
spectra  are  filtered  spectra  derived  from  the  taped  data.  This 
represents  editing  pass  1.  Voltage  power  spectral  density  (dB)  is 
plotted  against  frequency  (log  Hz).  For  each  record,  channel  1  (top) 
is  temperature  variance  and  channel  2  (bottom)  is  velocity  variance. 
Channel  number,  pass  number,  and  coded  flight  name  and  record  number 
are  shown  interior  to  each  graph.  The  slanted  line  represents  the 
regression  fit.  Regression  slope,  correlation  coefficient,  and 
percent  error  are  given  in  the  upper  left  corner  of  each  spectrum. 

The  same  three  dual -channel  records  in  flight  2  are  shown  as  in  Figure  1. 


VOLTAGE  SPECTRAL  DENSITY  <dB)  VOLTAGE  SPECTRAL  DENSITY  CdB) 


-  C-55  - 


Figure  5a.  (continued) 
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Figure  5b.  Voltage  variance  spectra  showing  editing  procedure.  These 
spectra  have  been  defiltered  with  the  atmospheric  filter  function. 

This  represents  editing  pass  2.  Voltage  power  spectral  density  (dB) 
is  plotted  against  frequency  (log  Hr).  For  each  record,  channel  1 
(cop)  is  temperature  variance  and  channel  2  (bottom)  is  velocity 
variance.  Channel  number,  pass  number,  and  coded  flight  name  and 
record  number  are  shown  interior  to  each  graph.  The  slanted  line 
represents  the  regression  fit.  Regression  slope,  correlation 
coefficient,  and  percent  error  are  given  in  the  upper  left  corner  of 
each  spectrum.  The  same  three  dual-channel  records  in  flight  2  are 
shown  as  in  Figure  1. 
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Figure  5b.  (continued) 


VOLTAGE  SPECTRAL  DENSITY  (dB)  VOLTAGE  SPECTRAL  DENSITY  (dB) 


-  C-59  - 


Figure  5b.  (continued) 
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Figure  5c.  Voltage  variance  spectra  showing  editing  procedure. 

Average  noise  level  has  been  subtracted  from  each  defiltered  spectrum. 
This  represents  editing  pass  3.  Voltage  power  spectral  density  (dB) 
is  plotted  against  frequency  (log  Hz).  For  each  record,  channel  1 
(top)  is  temperature  variance  and  channel  2  (bottom)  is  velocity 
variance.  Channel  number,  pass  number,  and  coded  flight  name  and 
record  number  are  shown  interior  to  each  graph.  The  slanted  line 
represents  the  regression  fit.  Regression  slope,  correlation 
coefficient,  and  percent  error  are  given  in  the  upper  left  corner  of 
each  spectrum.  The  same  three  dual -channel  records  in  flight  2  are 
shown  as  in  Figure  1. 
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Figure  5c.  (continued) 
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Figure  5c.  (continued) 
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through  pass  3.  An  average  noise  level  was  calculated  from  the  twenty 

lowest  spectral  density  values  and  subtracted  from  the  entire 

spectrum.  Figure  5c  shows  the  same  three  records  for  pass  3.  If  the 

percent  error  on  a  new  regression  fit  from  2  to  18  Hz  was  still 

greater  than  15%,  the  spectrum  was  rejected  on  pass  3  and  flagged. 

Spectra  accepted  on  one  of  the  above  three  passes  were  coded  as  "1". 

Flagged  spectra  were  then  manually  examined.  If  a  spectrum  contained 

a  visible  portion  of  the  inertial  subrange  that  was  represented  by  the 

regression  fit,  it  was  coded  "2".  Spectra  with  too  small  a 

signal- to -noise  ratio  (no  visible  inertial  subrange)  were  coded  "3". 

Unreal  spectra  were  coded  as  "4”.  These  unreal  spectra  were  steeply 

sloped  in  log-log  space.  No  theory  predicts  such  a  spectrum.  Thus, 

these  spectra  were  thought  to  be  byproducts  of  measurement  error, 

either  overdriving  the  filter  (i.e.,  amplitude  fluctuations  so  great 

that  filter  saturation  occurred)  or  noise  distortion  (e.g.,  change  of 

aircraft  power  setting  or  a  radio  transmission) .  Code  4  spectra  were 

eliminated  from  the  analysis.  All  other  analysis  was  performed  on 

code  1  and  2  spectra,  unless  otherwise  noted.  Table  2  gives  a 

representative  breakdown  of  spectra  by  editing  codes. 

Based  on  the  regression  analysis,  the  following  formulas  were 

used  to  derive  the  spectral  density  necessary  to  calculate  the 

2  2 

structure  parameters,  and  Cy  .  'X'  was  assigned  to 

2 

frequency  and  'Y'  to  Volts  /Hz  (spectral  density) .  Using  the 
power- law  equation, 


Y  -  a  Xm 
x5/3 


(6.1) 


Y  -  a  X 


(m  +  5/3) 


(6.2) 
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Table  2. 

Representative  editing  results  for  each  flight  and  in  total.  The 
breakdown  of  spectra  in  each  pass  for  code  1  is  shown,  as  well  as  the 
number  of  spectra  in  the  remaining  codes. 
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Table  2. 
(continued) 


Velocity  Spectra 
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TOTALS 

Temperature  Spectra 
Code  1  :  Pass  1  2 

222  25 


1  Code  2 
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Velocity  Spectra 
Code  1  :  Pass  121 
328  8  33 
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Q  -  X5/3  Y 

(6.3) 

a  -  Y  /  x“ 

(6.4) 

Q  -  10(<1°6  ^  '  m  <loS  10(«  +  5/3)  <log  X> 

(6.5) 

Q  _  10<<loS  Y>  +  5/3  Clog  X» 

(6.6) 

Q  _  10[(Z  log  Y)/N  +  1.667  (2  log  X)/N] 

(6.7) 

where  Q  is  the  spectral  density  times  frequency  to  the  5/3  and  is 

based  on  the  spectral  density  and  frequency  values  midrange  in  the 

2 

regression  fit.  Q  is  a  necessary  input  parameter  for  the  and 

2 

equations  (equations  5.23  and  5.24),  calculated  from  the 
temperature  and  velocity  turbulence  spectra,  respectively. 

6.2  Mean  Aircraft  Meteorological  Profiles 

The  aircraft  temperature  and  humidity  profiles  were  digitized  at 
levels  corresponding  to  the  profiler  high  resolution  range  gate 
heights.  Dewpoint  temperature  was  converted  to  vapor  pressure,  e, 
over  water  using  the  following  formula  (Iribarne  and  Godson  1981). 

log  e  -  -(2937.4  /  T)  -  4.9283  log  T  +  23.5471  (6.8) 

where  e  is  vapor  pressure  in  mb  and  T  is  in  Kelvin.  Vapor  pressure 
was  then  converted  to  specific  humidity,  q,  using  the  approximation 
e  *  (P  q)/0.622.  Potential  temperature  and  temperature  have  the 
defined  relationship 


9  -  T  (1000  mb/P) 


0.286 


(6.9) 
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with  pressure  from  the  standard  atmosphere  relationship 

P  -  PQ  e*(z/8)  (6.10) 

where  PQ  is  1013.2  mb,  altitude  z  is  in  km  and  scale  height  H  is  =*8 
km  (Iribarne  and  Godson  1981;  Wallace  and  Hobbs  1977).  Aircraft 
potential  temperature  profiles  are  shown  for  all  flights  in  Figure  6 . 
Aircraft  specific  humidity  profiles  for  only  the  PSU  flights  are  shown 
in  Figure  7. 

6.3  Other  Data  Sources 

Profiler  data  was  logged  at  one  hour  intervals  denoted  by  the  GMT 

time  (zulu  time,  "z")  at  the  end  of  the  average  period.  Height  is 

2 

given  in  m,  wind  speed  in  m/s,  direction  in  degrees,  and 
(radar)  in  m (Mike  Moss,  letter  to  author,  July  1987;  Mike  Moss, 
letter  to  author,  November  1987).  Radar  wind  profiles  for  the  PSU 
flights  are  shown  in  Figure  8. 

AFGL  thermosonde  data  was  available  at  20  m  resolution.  The 

following  information  was  available  for  each  level:  altitude  (km), 

pressure  (mb),  temperature  (°C) ,  relative  humidity  (%), 

-2/3 

(m  '  ),  wind  speed  (m/s),  and  direction  (degrees)  (Robert  Beland, 

letter  to  William  Syrett,  September  1986).  AFGL  thermosonde  potential 

temperature  profiles  corresponding  to  aircraft  flights  are  shown  in 

Figure  9,  Corresponding  AFGL  thermosonde  specific  humidity,  mean 

2 

wind,  and  calculated  raw  profiles  are  shown  in  Figures  10,  11, 
and  12,  respectively. 

AFGL  scintillometer  data  was  provided  as  seven- level  profiles 
2  -2/3 

with  Cn  in  m  '  ,  and  height  (corrected  for  zenith  angle)  in  km 
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Figure  6.  (continued) 
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Figure  6.  (continued) 
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Figure  6.  (continued) 
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Figure  6.  (continued) 
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Figure  7.  Specific  humidity  vertical  profiles  as  measured  by 
aircraft  instrumentation.  Altitude  is  MSL.  Only  flights  during  which 
the  PSU  profiler  operated  are  shown,  Flights  2-8  (labeled  above  each 
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Figure  7.  (continued) 
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Figure  7.  (continued) 
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RADAR  PROFILE  14  APRIL  07  15-16  Z 
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Figure  8.  Mean  wind  speed  (solid  line)  and  direction  (broken  line) 
measured  by  the  profiler.  Date  and  the  hour  interval  averaged  are 
given  above  each  graph.  Altitude  is  MSL. 
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Figure  11.  AFGL  the nnosonde -measured  vertical  profiles  of  mean  wind 
speed  (solid  line)  and  direction  (broken  line).  Altitude  is  MSL. 
Thermosonde  flight  number  is  given  above  each  graph. 
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(Murphy  and  Battles  1986) .  AFGL  scintillometer  profiles  corresponding 

to  aircraft  flights  are  shown  in  Figure  13. 

RADC  scintillometer  data  was  provided  as  seven- level  profiles 
2  -2/3 

with  C  in  m  '  at  the  standard  kernel  levels  (AGL  altitude 
n 

assuming  zenith  angle  of  zero)  (Donald  Stebbins,  letter  to  author, 
September  1986).  Only  one  data  run  corresponded  to  the  aircraft 
flight  times.  It  is  shown  in  Figure  14. 

NPS  rQ  scintillometer  data  was  provided  as  rQ  in  cm  as  a 
function  of  zulu  time  (Donald  Walters,  telephone  conversation, 

November  1987) .  Data  taken  during  aircraft  flight  was  limited  to  one 
night,  4  May  1986.  The  average  value  of  rQ  during  this  time 
segment,  0400-0540  2,  was  5  cm. 

AFWL  isoplanometer  data  was  unavailable  for  use  in  this 
comparison. 
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CHAPTER  7 

RESULTS  AND  DISCUSSION 
7.1  Analysis  of  C^/C^ 

2  2 

The  vertical  profiles  of  C^  and  C^  calculated  from  the 

aircraft  data  show  strong  correlation  for  flights  2,  6 (ascent  and 

descent),  8(ascent  and  descent),  9(ascent),  lO(descent) ,  11,  and  12 

(Figure  15).  The  other  flights  also  show  some  correlation,  but  not  at 

all  levels.  Correlation  of  higher  values  of  these  turbulence 

parameters  indicates  regions  of  the  atmosphere  that  contain  turbulent 

2  2 

layers  (Ottersten  1969).  Cj.  and  Cn  are  directly  correlated 

2 

at  all  levels  for  all  flights,  as  is  to  be  expected  since  Cn  is  a 
2 

function  of  Gj.  . 

A  rough  estimate  of  boundary  layer  height  can  be  determined  by 

spectral  signal- to-noise  ratio  for  these  profiles,  since  the  signal  is 

much  stronger  in  the  turbulent  boundary  layer.  For  this  purpose 

spectra  were  examined  in  detail  manually.  Table  3  gives  an  estimate 

2  2 

of  PBL  height  for  each  flight.  If  C^  versus  Cu  is  graphed 
(for  values  above  the  boundary  layer  where  the  TKE  and  variance  budget 
approximations  discussed  in  Chapter  3  are  valid),  a  power-law  relation 
is  evident  (Figures  16  and  17) .  This  ratio  is  proportional  to 
Richardson  number  and  temperature  gradient  (right  hand  side  of 
equation  3.28).  All  flights  except  6(ascent  and  descent),  7(descent), 
and  ll(ascent)  show  distinct  sloped  correlations.  All  flights  were 
plotted  simultaneously  and  a  regression  analysis  was  done  to  determine 
slope  and  goodness -of- fit  (Figure  18).  The  plot  of  flight  6  had  high 
scatter.  Removal  of  this  flight  decreased  the  scatter  and  thus 


CuA2  CT*2  (C *2/«*2/3> 


Figure  15.  Vertical  profiles  of  C  ,  ,  and  C  .  T^e 

top  graph  shows  the  C  profile  calculated  from  the 

profile.  (solid  line)  and  Cu  (broken  line)  are  shown  on 

the  bottom  graph  versus  MSL  altitude.  These  values  were  calculated 

from  the  aircraft  spectral  data.  The  flight  number  and  ascent/descent 

designation  for  the  profile  is  given  interior  to  the  top  graph.  Dual 

graphs  are  shown  for  flights  2-12. 
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Figure  15.  (continued) 
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Figure  15.  (continued) 
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Figure  15.  (continued) 
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Figure  15.  (continued) 


ALTITUDE  (km)  ALTITUDE  (k») 


ALTITUDE  <ka)  ALTITUDE  (ka> 


Figure  15.  (continued) 
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Figure  15.  (continued) 
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Figure  15.  (continued) 
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Figure  15.  (continued) 
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Table  3. 

Estimates  of  the  height  of  the  PBL  based  on  the  strength  of  the 
signal -to -noise  ratio  exhibited  in  the  spectra. 


night 

2 

PBL  Height 

cmsl  ian) 
1.3 

3 

2.3 

4 

1.5 

6 

2.2 

7 

2.2 

8 

2.4 

9 

2.2 

10 

2.8 

11 

2 

12 

1 
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Figure  16.  (continued) 
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Figure  16.  (continued) 
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Figure  16.  (continued) 
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Figure  L6 .  (continued) 
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Figure  16.  (continued) 
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Figure  16.  (continued) 


CT~2  <K^2  m^-2/3) 


-  C-110  - 


Figure  16.  (continued) 
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Figure  17.  (concinued) 
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Figure  18.  (continued) 
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increased  correlation  of  the  entire  data  base  (Figure  18).  The  slope 

for  all  the  vertical  profiles  was  close  to  3/5.  Correlation  and  slope 

for  all  the  level  flights  showed  approximately  the  same  scatter 

(correlation)  and  slope  as  the  vertical  profiles  (Figure  19). 

2  2 

In  Figure  20,  was  plotted  versus  Cu  multiplied  by 

9/ g  (39 /3z)  -  T/g  (3T/3z  +  D  (7.1) 

for  all  vertical  flights  (having  associated  temperature  profiles). 
Again,  flight  6  had  high  scatter  and  flight  7  showed  little  slope.  A 
combined  plot  of  all  flights  showed  slightly  higher  correlation 
without  flight  6  and  again  roughly  a  slope  of  3/5  (Figure  20) .  This 
ratio  corresponds  to  values  of  1.6  Chat  range  from 

0.3  for  high  values  of  to  as  much  as  10  for  low  values  of 


2  2 

As  an  additional  comparison,  /  (Cu  9/ g  39/3 z)  was 
plotted  against  a  turbulence  activity  parameter  (Gregg  1987)  ,  e/i/N 
(Figure  21) .  N  was  calculated  as  g  31n9/3z  from  a  polynomial  fit  to 
the  aircraft  potential  temperature  profile.  The  formula  used  to 
derive  kinematic  viscosity  was  taken  from  the  U.  £.  Standard 
Atmosphere  Supplements ,  1966  (Environmental  Science  Services 
Administration  et  al.  1966). 

H  -  (0  T3/2)/(T  +  S)  (7.2) 

where  m  is  the  dynamic  viscosity,  0  -  1.458E-6  kg  sec  ^K  . 

and  S  -  110.4  Kelvin.  Kinematic  viscosity,  v ,  equals  n/p.  Density 

.  3 

was  calculated  as  p  -  pQ  exp(-Z/8  km)  where  pQ  is  1.225  kg  m 

2  2 

T  was  taken  from  the  aircraft  profiles.  Gp  ,  Cu  ,  e ,  and  N 
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Figure  19.  A  composite  comparison  of  C^  versus  Cu  for  all 
level  profiles.  The  broken  line  represents  a  regression  analysis 
Slope  and  correlation  coefficient  are  given  below  each  graph. 
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2  7 

Figure  20.  A  composite  comparison  of  C^  versus  C 
normalized  by  the  expression  [T/g  ( dT/az  +  D ] .  T^e  broken  line 
represents  a  regression  analysis.  Slope  and  correlation  coeffic 
are  given  below  each  graph. 
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Figure  21.  (continued) 
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Figure  21.  (continued) 


CW  ✓  [Cu''2#(T/'g  )J»<d7VdZ  *  r>3  CW  X  CCu~2#nVg)*<dTVdZ  +  r)3 


-  C- 122  - 

FLIGHT  8  ASCENT 


Figure  21.  (continued) 
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Figure  21.  (continued) 
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Figure  21.  (continued) 
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were  all  calculated  from  the  aircraft  data.  Figure  22  is  a  composite 

of  all  vertical  flights.  This  plot  clearly  shows  that  as  turbulence 

2  2 

activity  goes  down,  the  ratio  of  to  normalized  Cu  goes 

2  2'  2 

up.  For  activity  values  above  10  ,  the  /Cu  (norm)  ratio 
becomes  constant  at  about  0.4. 

7.2  Mean  Gradients  and  Richardson  Number 

Next,  the  Richardson  number  was  evaluated.  Since  VanZandt  et  ai. 

(1981)  defined  R^  as  N/S ,  N  versus  S  was  plotted,  thus  showing 

Richardson  number  as  a  slope  (Figure  23).  N  was  calculated  by  a 

finite  difference  derivative  between  each  level  in  the  digitized 

aircraft  potential  temperature  profiles,  while  S  was  calculated  by 

finite  difference  from  the  profiler  hourly  average  taken  during  the 

time  the  aircraft  was  collecting  data  (Table  1) .  S  was  calculated  as 
2  2 

(du/dz)  +  (5v/3z)  .  Figure  24  is  a  composite  of  all  flights 

shown  in  Figure  23.  Most  values  show  R^>  0.25  above  the  broken 

line,  with  the  shear  term  on  the  order  of  10  ^  and  the  buoyancy  term 
-  5  -4 

ranging  from  10  to  4  X  10  Richardson  number  profiles  for 
each  flight  were  directly  calculated  (Table  4) ,  and  then  averaged 
(Table  5)  . 

These  values  of  N  and  S  are  the  necessary  input  parameters  for 

the  VanZandt  et  al.  model  which  was  used  to  predict  vertical  profiles 
2 

of  Cn  and  e.  This  was  done  for  each  flight  conducted  at  Penn 

State  where  profiler  data  was  available  (Figure  25).  The  model  outtu: 

2 

is  normalized  optical  C .  Standard  pressure  calculated  from  the 
geometric  height  and  the  aircraft  temperature  profile  were  used  to 
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Figure  23.  Comparisons  of  Che  VanZandc  et  al.  parameters  N.  the 
Brunt-Vaisala  frequency  squared,  and  S,  Che  shear  squared.  Compariso 
is  shown  for  each  vertical  aircraft  flight  during  which  the  radar 
operated  (labeled  above  each  graph).  N  was  calculated  using  aircraft 
potential  temperature  profiles  and  S  was  calculated  using  radar  wind 
profiles.  Note  units.  Broken  line  represents  R.  -  0.25. 
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Figure  23.  (continued) 
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Figure  24.  A  composite  comparison  of  N  versus  S  for  all  flights  shov 
in  Figure  23.  Note  units.  Broken  line  represents  R,.  -  0.25." 
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Table  4. 

Calculated  vertical  profiles  of  Richardson  number  based  on  R.-N/S 
and  using  N  and  S  values  shown  in  Figure  17.  Altitude  shown1 is  MS 
and  is  the  average  altitude  between  each  vertical  radar  observacio 
level . 


Flight  2 

Flight  3 

Alt  (km) 

Alt  (km) 

Si 

1.4 

1.2 

1.7 

9.04 

1.7 

2.61 

2 

20.99 

2 

3.92 

2.3 

5.18 

2.3 

31.98 

2.6 

6.99 

2.6 

6.73 

2.9 

2.45 

2.9 

.19 

3.2 

88.53 

3.2 

.56 

3.4 

1.58 

3.4 

18.94 

3.7 

11.59 

3.7 

74.74 

4 

0 

4 

1.4 

4.3 

.07 

4.3 

.21 

4.6 

4.51 

4.6 

1.68 

4.9 

10.18 

5.2 

1.15 

5.5 

11.4 

5.8 

1.17 

6.1 

2.56 

6 . 4 

1.83 

6.6 

8.96 

6.9 

2.39 

7.2 

0 

7.5 

2.02 
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Table  4. 
(continued) 


Flight  4 
AlC  (kjn) 

Si 

Flight  6 
Alt  (km) 

1.4 

2.45 

1.4 

-  .03 

1.7 

34.24 

1.7 

142.19 

2 

86.56 

2 

.78 

2.3 

5.99 

2.3 

21.71 

2.6 

5.87 

2.6 

8.65 

2.9 

1.09 

2.9 

179.84 

3.2 

-.02 

3.2 

1.6 

3.4 

39.34 

3.4 

3.2 

3.7 

3.23 

3.7 

181.64 

4 

6.27 

4 

7.6 

4.3 

6.21 

4.3 

2.24 

4.6 

10.47 

4.6 

9.6 

4.9 

6.66 

4.9 

27.89 

5.2 

3.32 

5.2 

4.13 

5.5 

5.11 

5.5 

1.92 

5.8 

1.84 

5.8 

2.19 

6.1 

9.16 

6.1 

.05 

6.4 

31.57 

6.4 

-  .02 

6.6 

138.56 

6.6 

-  .02 

6.9 

2.49 

6.9 

.19 

7.2 

1.53 

7.2 

.25 

7.5 

1.6 

7.5 

0 

7.8 

-1.53 
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Table  4. 
(continued) 


Flight  7 

Flight  8 

Gsm) 

Alt  (km) 

1.4 

-  .05 

1.4 

.02 

1.7 

-  .21 

1.7 

0 

2 

1.82 

2 

2.31 

2.3 

1.81 

2.3 

8.51 

2.6 

1.28 

2.6 

8.05 

2.9 

.82 

2.9 

8.64 

3.2 

3.21 

3.2 

8.93 

3.4 

1.67 

3.4 

23.48 

3.7 

32.94 

3.7 

17.85 

4 

.13 

4 

-1.22 

4.3 

.53 

4.3 

1.22 

4.6 

3.81 

4.6 

8.08 

4.9 

4.54 

4.9 

67.55 

5.2 

4.94 

5.2 

3.82 

5.5 

2.02 

5.5 

6.6 

5.8 

10.74 

5.8 

5.3 

6.1 

11.11 

6.1 

7.18 

6.4 

1.56 

6.4 

11.93 

6.6 

4.67 

6.6 

23.27 

6.9 

6.6 

6.9 

-8.82 

7.2 

1.9 

7.2 

3.68 

7.5 

1.1 

7.5 

7.03 

7.8 

1.99 

7.8 

3.3 
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Table  5. 

The  arithmetically  averaged  Richardson  number  at  each  height  for  al 
vertical  flights. 


Average  Richardson  Number 


(las) 

1.4 

0.72 

1.7 

31 

2 

19 

2.3 

12 

2.6 

6 

2.9 

32 

3.2 

17 

3.4 

15 

3.7 

54 

4 

2.36 

4.3 

1.75 

4.6 

6 

4.9 

23 

5.2 

3.67 

5.5 

5 

5.8 

4.25 

6.1 

6 

6.4 

9 

6.6 

35 

6.9 

0.57 

7.2 

1.47 

7.5 

2.35 

7.8 

1.25 
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Figure  25.  Comparison  of  che  VanZandc  et  al.  model  vertical 
prediction  of  -optical  (solid  line)  and  the  aircraft 
CT,  *derlved  vertical  profile  of  Cn  -optical  (broken  line). 

PjU  flights  are  shown  (labeled  interior  to  each  graph)'. 
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"de-normalize"  Che  model-derived  parameter.  Since  only  the  high 
resolution  profiler  range  gates  were  used  to  calculate  S  at  a 
resolution  of  290  m,  the  prediction  was  only  calculated  from  1410  m 
MSL  to  7515  m  MSL.  An  average  MSL  altitude  for  State  College,  Pa,  is 
300  m.  Except  for  flight  6(ascent)  the  model  and  aircraft  profiles 
agreed  to  within  an  order  of  magnitude.  This  particular  discrepancy 
in  flight  6  is  thought  to  be  the  result  of  errors  in  gain  records  for 
the  flight. 

7.3  Profiles  of  C  ^ 
n 

2 

The  VanZandt  et  al.  model  can  also  be  used  to  give  radar  . 

2 

The  relationship  between  profiler  measured  Cn  -radar  and 
2 

-optical  is  a  function  of  temperature  and  humidity  (Tatarskii 

2 

1971).  This  parameter,  previously  referred  to  as  X  ,  has  been 

calculated  using  the  aircraft  specific  humidity  and  temperature  data. 

2  2  2  2 
X  is  the  ratio  of  C  -radar  to  C  -optical.  X  at  each 

n  n  r 

aircraft  average  observation  level  was  calculated  with  a  polynomial 

fit  to  the  temperature  profile  and  a  finite  difference  derivative  from 

the  humidity  profile.  Smoothing  of  the  vertical  profiles  was  achieved 

by  setting  negative  or  zero  temperature  derivatives  equal  to  0.00001 

K/m  and  positive  specific  humidity  profiles  equal  to  zero.  This 

2 

calculation  is  given  in  Table  6.  X  at  each  average  profiler  high 

resolution  gate  height  was  also  calculated  with  finite  difference 

derivatives  from  both  temperature  and  humidity  digitized  profiles. 

Smooching  was  performed  as  before.  This  calculation  is  given  in  Table 

2 

7.  An  arithmetically  averaged  X  profile  from  Table  7  is  given  in 
Table  8.  High  humidity  in  the  boundary  layer  gives  extremely  high 
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Table  6. 

vertical  profiles  for  each  flight  during  which  the  profiler 
operated,  calculated  at  each  aircraft  observation  level  in  MSL 
altitude  with  aircraft-measured  temperature  and  specific  humidity 
profiles . 


Flight  2 

Descent 

Flight  3 

Descent 

Alt  (km) 

z1 

(km) 

x2 

4.4 

1 

9 

1 

3.8 

1 

8.7 

1 

3 

1 

8.4 

1.6 

2.3 

1.2 

8 

1.3 

1.8 

1.1 

7.6 

1.3 

1.4 

42.8 

7.2 

1.2 

1 

25.5 

6.6 

1.4 

.7 

12.5 

6.3 

1.6 

.5 

148.6 

5.8 

1.9 

5.4 

1 

5 

1 

4 . 6 

1.5 

4.1 

1.8 

3.7 

2.3 

3.4 

1.7 

3 

1.9 

2.6 

1.5 

1.7 

12.4 

1.4 

2.6 

1 

2.3 

.6 

1914 

.4 

13019 

.8 

3717 

.3 

12907 
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Table  6. 
(continued) 


Flight  4 

Descent 

Flight  6 

Descent 

Ale  (knj) 

X2 

Alt  (km) 

£2 

9 

76189 

9 

1 

8.7 

106301 

8.5 

1.3 

8.4 

104241 

8.2 

1.2 

8.2 

232623 

7.9 

6.4 

7.4 

1 

7.5 

1.1 

6.8 

1 

7 

1.8 

6.1 

1 

6.6 

18.9 

5.6 

1 

6.3 

4.6 

5.1 

4 . 6 

5.9 

1 

4.8 

1 

5.4 

1 

4.3 

3.8 

5 

1.9 

4 

1.1 

4.4 

12.8 

3.5 

1 

3.8 

2.4 

3 

1.8 

3.4 

1 

2.7 

1.6 

2.9 

13.4 

2.2 

1.8 

2.4 

3.2 

1.9 

3.4 

2 

13.9 

1.5 

227.8 

1.4 

7.7 

.8 

1.4  . 

1.2 

7.7 

.7 

2.33E+6 

.8 

1.3 

.6 

2.32E+6 

.7 

4.2 

.5 

2.3E+6 

.6 

3.7 

.5 

3.2 
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Table  6. 
(continued) 


Flight  7 

Descent 

Flight  8 

Descent 

till  (Jsm) 

Alt  (kin) 

8.5 

1 

8.4 

1.9 

8.3 

1 

8.1 

3.5 

7.9 

2.6 

7.7 

3.3 

7.6 

3.5 

7.3 

10.8 

7.2 

13.4 

7 

11.6 

6.8 

4.1 

6.7 

3.8 

6.4 

39.2 

6.3 

21.7 

5.9 

3.3 

5.8 

15.3 

5.7 

17.8 

5.3 

6.1 

5 

24.6 

4.9 

1.2 

4.6 

5.8 

4 . 6 

54.8 

4.2 

4 

4.1 

1.2 

3.8 

1.3 

3.8 

1.1 

3.3 

1 

3.3 

87.5 

2.9 

1 

2.8 

6.6 

2.3 

21.4 

2.6 

6.8 

1.9 

81.9 

2.2 

101.9 

1.6 

1.9 

1.7 

5.8 

1 

351269 

1.4 

2.2 

.8 

34447 

1 

1.6 

.6 

16 

.8 

1.6 

.5 

3.2 

.6 

1.5 

.5 

1.5 
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2  Table  7 . 

Yr  vertical  profiles  for  each  flight  during  which  the  profiler 
operated,  calculated  at  each  profiler  observation  level  in  MSL 
altitude  with  aircraft-measured  temperature  and  specific  humidi 
protiles  . 


Flight  2 

Flight  3 

9 

Alt  (km) 

£ 

Alt  (Km) 

x2 

1.3 

46.9 

1.3 

5.7 

1.6 

3.5 

1.6 

2.5 

1.9 

1.1 

1.9 

11.9 

2.1 

1.3 

2.1 

138.2 

2.4 

1.2 

2.4 

3.2 

2.7 

1.1 

2.7 

1.5 

3 

1 

3 

1.9 

3.3 

3.3 

1.6 

3.6 

1 

3.6 

2.2 

? .  9 

1 

3.9 

2.2 

4.2 

1 

4.2 

1.8 

4.5 

1 

4.5 

1.5 

4.8 

1.2 

5 

1 

5.3 

l 

5.6 

25.5 

5.9 

1.9 

6.2 

1.6 

6.5 

1.4 

6.8 

1.4 

7.1 

1.2 

7.4 

1.2 

7.7 

1.3 

8 

1.3 
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Table  7 . 
(continued) 


Flight  4 

X1 

Flight  6 

„2 

Alt  Qsd) 

(ism) 

Jr 

1.3 

34.9 

1.3 

7.8 

1.6 

210.3 

1.6 

2.9 

1.9 

3.4 

1.9 

14.9 

2.1 

1.8 

2.1 

56.2 

2.4 

1.9 

2.4 

3.2 

2.7 

1.6 

2.7 

1 

3 

1.8 

3 

13.2 

3.3 

1.3 

3.3 

1 

3.6 

1 

3.6 

1.2 

3.9 

1.1 

3.9 

2.5 

4.2 

3.8 

4.2 

1.1 

4.5 

6.9 

4.5 

13.9 

4.8 

1 

4.8 

4.2 

5 

4.6 

5 

2 

5.3 

2.8 

5.3 

1 

5.6 

1 

5.6 

5.9 

5.9 

1 

5.9 

1 

6.2 

1 

6.2 

4.2 

6.5 

1.9 

6.5 

20.4 

6.8 

1 

6.8 

6.1 

7.1 

1 

7.1 

1.7 

7.4 

1 

7.4 

1.1 

7.7 

1.5 

7.7 

2.4 

8 

1.1 

8 

6.2 
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Flight  7 

Al£  (Jsn) 

1.3 

1.6 

1.9 

2.1 

2.4 

2.7 
3 

3.3 

3.6 

3.9 

4.2 

4.5 

4.8 
5 

5.3 

5.6 

5.9 
6.2 
6.5 

6.8 
7.1 

7.4 

7.7 
8 


Table  7. 
(continued) 


0 

Flight  8 

X* 

Alt  (km) 

1.2 

1.3 

2.3 

1.9 

1.6 

1.6 

12.9 

1.9 

5.4 

27.4 

2.1 

107.4 

6.5 

2.4 

7.7 

1 

2.7 

6.7 

1 

3 

112.1 

1 

3.3 

87.6 

1.2 

3.6 

9.6 

1.3 

3.9 

1.1 

4 

4.2 

1.2 

5.5 

4.5 

50.9 

1.2 

4.8 

1.2 

25.6 

5 

1.2 

1.5 

5.3 

6.1 

16.3 

5.6 

25.8 

3.3 

5.9 

15.3 

4.7 

6.2 

22.4 

41.3 

6.5 

10.3 

4.1 

6.8 

3.7 

14.8 

7.1 

11.1 

6.5 

7.4 

10.8 

3.5 

7.7 

3.2 

2.6 

8 

3.1 
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Table  8 . 

Arithmetically  averaged  X^  vertical  profiles  from  the  Table  7 
profiles . 

2 

Average  Radar  X  Profile 


All  (kn) 

1.3  16.5 

1.6  37.1 

1.9  23.3 

2.1  55.4 

2.4  3.9 

2.7  2.1 

3  21.8 

3.3  15.6 

3.6  2.7 

3.9  1.5 

4.2  2.2 

4.5  13.3 

4.8  1.8 

5  6.9 

5.3  2.5 

5.6  14.9 

5.9  4.5 

6.2  6.8 

6.5  15 

6.8  3.3 

7.1  6 

7.4  4.1 

7.7  2.4 

8  2.9 
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2 

values  of  C  -radar  at  these  levels.  Aircraft  profiles  of 
n 

2  2  2 
Cn  -optical  adjusted  with  X  to  give  CR  -radar  show 

2 

agreement  with  the  profiler  derived  Cn  -radar  profiles  (Figure 
26). 

For  the  two  night  flights  during  which  multiple  instruments  took 

2 

data,  comparisons  of  Cn  -optical  vertical  profiles  are  shown.  The 
2 

level  at  which  C  was  measured  by  the  AFGL  scintillometer  varied 
n 

because  the  instrument  was  measuring  at  different  zenith  angles.  To 

aid  in  graphical  interpretation,  for  each  night,  all  the  AFGL 

2 

scintillometer  profile  Cn  values  were  log  averaged  and  the 

heights  arithmetically  averaged,  at  each  level.  RADC  data  (Figure  14) 

was  available  for  one  night  only  and  was  omitted  for  clarity.  Figure 

2 

27  shows  a  comparison  of  Cn  -optical  from  the  averaged  AFGL 

scintillometer  profiles,  the  aircraft,  the  VanZandt  et  al.  model 

predictions,  the  AFGL  thermosonde,  and  the  profiler  (unadjusted  raw 
2 

Cn  -radar).  The  ascent  profile  for  the  aircraft  on  flight  6  is 
thought  to  have  been  an  incorrect  gain  setting  which  leads  to  the 
unrealistically  high  values  in  the  lower  part  of  the  profile.  Figure 

28  is  the  same  comparison  except  that  the  radar  profile  has  been 

2  2 
converted  to  -optical  using  the  calculated  X  profiles  of 

Table  7.  These  plots  show  again  the  agreement  between  the  model  and 

the  aircraft,  while  the  thermosonde  has  good  agreement  with  the 

2 

scintillometer.  The  profiler  CR  -optical  values  (from  the  hourly 
average  observation  during  the  aircraft  flight)  show  some  agreement 


with  the  thermosonde . 


10 


I 


9  Aar- ived  from  profile 

t.  26-  comparison  •'  .ulcipU**  by  =he 

in  line)  end  t«»  aircr.l t  5?r,dar  (i„ken  line)  •  ?SU 

^labeled  incSrios  to  each  Staph,- 


ALTITUDE  (ka)  ALTITUDE  (ka> 


5 


Figure  26.  (continued) 


ALTITUDE  (km)  ALTITUDE  (k*> 


C-152 


Figure  26.  (continued) 


Cn*2  optical  <»*-2/3> 


2 

Figure  27.  Comparison  of  vertical  profiles  of  C  -optical. 
Altitude  is  MSL.  Bars  represent  standard  deviation  of  AFGL 
scintillometer  data.  Profiler  values  are  raw  C  -radar  (not 
converted  to  Cn  -optical).  Aircraft  flight  designations  are 
labeled  interior  to  each  graph. 


Log  averaged  scintillometer  - j - 

VarvZandt  et  al.  model  - - 

Arithmetically  smoothed  thermosonde 


Aircraft 

Profiler 


56  - 


FLIGHT  6 
RSCENT  PROFILE 


•19  ia-18  ia-17 


10-16  10‘13  10M 


Cn*2  optics!  (**-2/3) 


FLIGHT  6 


PROFILE 


Cn*2  optics!  (■—2/3) 


Figure  28.  Comparison  of  vertical  profiles  of  -optical. 
Altitude  is  MSL.  Bars  represent  standard  deviation  of  AFGL 
scintillometer  data.  Profiler  values  have  been  converted  to 
C  -optical.  Aircraft  flight  designations  are  labeled  interior 
each  graph. 


Log  averaged  scintillometer  - 1 - 

VanZandt  et  al.  model  - 

Arithmetically  smoothed  thermosonde 


Aircraft 

Profiler 
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7.4  Propagation  Parameters 


Other  optical  parameters  of  interest  are  the  coherence  length 

rQ  and  isoplanatic  angle  SQ.  These  are  related  to  vertical 

2 

integrals  of  C (Murphy  and  Battles  1986) 


r 


o 


o 


-3/5 


(7.3) 


8 


o 


ctf 


Cn(*> 


z5/3dz  f3/5 


(7.4) 


1  2 

where  kc  and  kr  are  functions  of  wavelength  (A  *  dependence) . 

Assuming  a  wavelength,  and  a  linear  relation  between  observations  of 

Cn  -optical  at  each  level,  the  integrals  can  be  evaluated 

numerically.  But  contributions  from  above  and  below  the  measured 

profile  must  be  calculated  based  on  theoretical  assumptions.  In  the 

2 

boundary  layer,  Murphy  and  Battles  (1986)  assumed  Cn  to  go  as 


<£<*>  "  CJ<zb)  <*/*b>'2/3 


(7.5) 


where  z^  is  the  lowest  level  of  the  measured  profile.  In  the  upper 
2 

atmosphere,  Cn  was  modeled  by  an  exponential  decrease 


«n<*> 


r2f~  1  -a(z-zr> 
"  C  (z  )  e  c 


where  a  -  (In  10)  DR,  and  DR  (the  atmospheric  drop  off  rate)  -  1.3  X 
10'4.  With  these  relationships,  the  integral  portions  for 
above  and  below  become  (Murphy  and  Battles  1986) 


Above  -  z3/3  (C3(zc)/a)  [1  +  (5/3)(azc) 


(10/9) 


(7.7) 
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Below-  \  z ®/3  C^(Zb)  (7.8) 

and  for  r 

o 

Above  -  C^(zc)/a  (7.9) 

Below  -  3  zfe  C^(zb)  [1-  (6/zb)1/3]  (7.10) 

As  part  of  Che  AFGL  thermosonde  flights,  rQ  and  SQ  were 

calculated  using  equations  7.3  and  7.4  linearly  integrated  for  just 

the  thermosonde  profile,  assuming  \  -  .3  n  (Robert  Beland,  letter  to 

William  Syrett,  September  1986;  Robert  Beland,  telephone  conversation. 

November  1987).  These  values  are  shown  in  Table  9.  r  and  9 

o  o 

were  also  calculated  by  Murphy  and  Battles  (1986)  with  all  the  above 
integral  equations  and  theoretical  assumptions  for  the  atmosphere 
above  and  below,  for  the  AFGL  scintillometer  data,  using  0.5  m-  An 
average  value  for  each  night  is  given  in  Table  10.  No  rQ  values  are 
shown  for  the  AFGL  scintillometer  because  the  boundary  layer  has  a 
large  contribution  to  the  tQ  integral.  The  lowest  scintillometer 
measurement  level  is  roughly  2  km  AGL.  Comparisons  of  contributions 
to  the  total  integral  showed  that  the  theoretical  extrapolations  belov 
the  measured  profile  were  contributing  the  major  percentage  to  the 
computed  value.  Calculated  rQ  values  were  thus  considered 
unreliable  (Murphy  and  Battles  1986).  To  facilitate  comparison,  this 
same  method  and  assumptions  were  used  to  compute  rQ  and  from 
the  aircraft  data  at  0.3  (j  and  0.5  n  wavelength  (Table  11).  Both  the 
aircraft  and  thermosonde  profiles  start  at  much  lower  AGL  altitudes 
and  thus  the  measured  profile  is  a  significant  contributor  to  the  rQ 
calculations  for  both  these  instruments  (Rcoert  Beland,  telephone 
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Table  9. 

AFGL-calculated  optical  parameters  from  thermosonde  da 

AFGL  Thermosonde  Optical  Data 
(calculated  with  A-.3 m) 

Flight  ro  (SJS)  *0  (M£ad) 

L4033  4.3  3.2 


L4019 


4.4 


4.5 
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Table  10. 

Average  of  AFGL-calculated  optical  parameters  from  scintillomete 
data . 


AFGL  Scintillometer  Data 
(calculated  wich  A-.5 fi) 


Date 

SB I  lime 

e0  (Miai) 

Percent 

:  Contribution 

5/4/86 

0342-0513 

8.02 

a 

56.5 

b 

42.2 

c 

1.3 

5/6/86 

0317-0420 

10.29 

a 

61.4 

b 

36.8 

c 

1.8 

a-  Scintillometer  Measurement 
b-  Atmosphere  Above 
c-  Atmosphere  Below 
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Table  11. 

Optical  parameters  from  aircraft  data  calculated  with  the  Murphy  ar.c 
Battles  assumptions. 


Aircraft  Data 
(calculated  with  A-.5m) 


Date 

flight; 

90  (Miai) 

Percent  Contribution 

4  May 

6 

8.3 

a 

84.8 

b 

15 

c 

0.2 

6  May 

8 

7.9 

a 

49.2 

b 

50.7 

c 

0 

Aircraft  Data 

(calculated  with  A-, 

3m) 

Date 

Flight; 

oq  (Miah) 

Percent  Contribution 

4  May 

6 

4.5 

a 

84.8 

b 

15 

c 

0.2 

6  May 

a 

4.3 

a 

49.2 

b 

50.7 

c 

0 

Date 

Flight; 

rQ  (fiffl) 

Percent  Contribution 

4  May 

6 

2.4 

a 

68.5 

b 

0.5 

c 

31 

6  May 

8 

4.8 

a 

82.2 

b 

6 . 4 

c 

11.3 

a-  Therraosonde  Measurement 
b-  Atmosphere  Above 
c-  Atmosphere  Below 
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conversation,  November  1987).  The  9  integral,  inversely,  weights 

the  upper  atmosphere  most  heavily.  Since  the  aircraft  is  limited  in 

altitude  to  “10  km,  the  contribution  to  0  is  limited.  However,  due 

2 

to  the  integral  approximation,  profiles  in  which  Cn  drops  off 

sharply  with  altitude  will  weight  the  measured  layer  more  heavily. 

This  is  seen  in  the  profile  of  flight  6.  0Q  values  from  the 

scintillometer  and  the  aircraft  (at  A  -  0.5  /*)  were  similar  on  4  May, 

but  disagreed  on  6  May.  The  thermosonde  and  aircraft  6q  values  (at 

A  -  0.3  n)  disagreed  on  4  May,  but  agreed  on  6  May. 

rQ  was  also  measured  at  optical  wavelengths  during  the  EWAK 

experiment  by  the  NPS  rQ  scintillometer.  The  measured  rQ  average 

for  the  aircraft  flight  tim*  on  the  night  of  4  May  1986  was  5  cm.  The 

thermosonde  and  aircraft  rQ  values  (at  A  -  0.5  m)  differed  by  a 

factor  of  two  on  4  May.  The  night  of  6  May,  the  thermosonde  and 

aircraft  r  values  were  similar, 
o 

7.5  VanZandt  et  al.  Model  for  e 

The  VanZandt  et  al.  model  e  predictions  are  compared  to  e 

profiles  from  the  aircraft  data  (Figure  29) .  The  hot  wire  failed 

during  flight  4,  so  that  only  a  small  amount  of  data  is  available. 

The  model  shows  predictions  generally  ranging  in  the  area  of  5  X  10--* 

for  all  flights.  There  is  poor  agreement  between  the  aircraft  and 

model,  the  aircraft  (except  for  flight  3)  often  giving  values  three  or 

four  orders  of  magnitude  smaller.  Aircraft  values  also  vary  much 

more,  ranging  from  10'^  to  10'^. 

A  crude  average  profile  of  rj,  the  microscale,  based  on  the 

3  25 

aircraft  t  profiles  is  shown  in  Figure  30.  »>  -  (w  /«)* 


was 


Figure  29.  Comparison  of  vertical  profiles  of  epsilon  (f),  the 
turbulent  kinetic  energy  dissipation  rate,  predicted  by  the  VanZandt 
et  al.  model  (solid  line)  and  measured  by  the  aircraft  (broken  line) 
Altitude  is  MSL.  Aircraft  flight  designations  are  given  interior  to 
each  graph. 
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Figure  29.  (continued) 
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Figure  30.  Average  vertical  profile  of  eta  (tj)  ,  the  turbulence  inne 
size  scale,  from  aircraft-measured  epsilon  profiles.  Altitude  is  MS 
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computed  for  each  t  value  using  the  previously  mentioned  formula  for 
with  the  standard  atmosphere  temperature  profile,  and  then  averaged 
1  km  layers.  The  average  value  of  r?  spans  one  order  of  magnitude, 
from  1  mm  in  the  boundary  layer  to  an  average  of  1  cm  above  the 
boundary  layer. 


-  C— 170  - 


CHAPTER  8 

SUMMARY  AND  CONCLUSIONS 


The  atmosphere  has  a  temporally  and  spatially  variable  refractive 

index  due  to  its  inherently  turbulent  nature.  This  can  affect 

electromagnetic  beam  propagation  by  introducing  beam  steering,  image 

dancing,  beam  spread,  spatial  coherence  degradation,  temporal 

coherence  degradation,  and  scintillation  (Dewan  1980) .  Refractivity 

is  wavelength  dependent.  At  microwave  radar  wavelengths,  molecular 

composition  (e.g.,  water  content)  dominates  refractivity  changes.  At 

visible  wavelengths,  density  (parameterized  by  temperature)  dominates 

(Balsley  and  Gage  1980).  For  many  quantification  purposes,  the 

2 

refractive  index  turbulent  structure  constant,  ,  is  a  key 

parameter.  This  parameter  can  be  inferred  from  optical  turbulence 

(scintillometer) ,  temperature  turbulence  (aircraft  instrumentation  or 

thermosonde)  or  radar  backscatter  (profiler)  measurements.  However, 

2 

Cn  from  a  radar  profiler  will  differ  because  of  the  wavelength 

dependence  of  refractivity  Theoretically,  this  can  be  compensated 

for,  if  the  humidity  profile  is  known  (VanZandt  et  al.  1981). 

2 

Vertical  integrals  of  Cn  give  transverse  coherence  length,  rQ , 
and  isoplanatic  angle,  8  .  These  can  also  be  measured  remotely 
(based  on  optical  turbulence)  by  an  rQ  scintillometer  and 
isoplanometer ,  respectively. 

During  an  atmospheric  optics/meteorology  experiment  (acronym 
EWAK)  conducted  at  Penn  State  University  primarily  during  April  and 
May  of  1986,  data  was  collected  by  the  aforementioned  instruments. 
AFGL  and  RADC  scintillometers  produced  vertical  profiles  of 
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2 

Cn  -optical.  An  AFGL  thermosonde  produced  vertical  profiles  of 
2 

Cj  and  other  meteorological  variables.  A  PSU  profiler  produced 

2 

vertical  profiles  of  wind  direction  and  speed,  and  Cn  -radar.  The 

instrumented  ARA  research  aircraft  produced  vertical  profiles  of 

2  2 

temperature  and  velocity  turbulence  (C^  and  Cu  )  and  other 
meteorological  variables. 

An  indepth  analysis  of  the  aircraft  turbulence  data  was 

performed.  The  turbulence  instrumentation  aboard  consisted  of 

cold-wire  and  hot-wire  sensors  (for  temperature  and  velocity 

variations,  respectively)  and  FM  recording  apparatus.  The  taped  data 

was  processed  via  FFT  to  produce  one -dimensional  variance  spectra 

(wavenumber  range  0.01  to  10  m"^).  Flights  usually  produced  a  10  km 

vertical  profile;  the  data  was  processed  to  give  roughly  0.5  km 

vertical  resolution,  similar  to  that  of  the  profiler. 

An  atmospheric  filter  function  was  derived,  based  on  PBL  spectra 

and  -5/3  inertial  subrange  slope.  This  function  was  applied  to 

compensate  for  low  frequency  dropoff,  as  part  of  a  spectral  editing 

program.  Editing  was  based  on  percent  error  between  a  regression 

analysis  and  theoretical  -5/3  slope.  Noise  subtraction  was  performed 

in  some  low  signal  cases.  In  general,  PBL  spectra  clearly  showed  the 

classical  -5/3  inertial  subrange  slope,  out  to  the  broadband 

instrument  noise  level.  A  majority  of  spectra  in  the  free  troposphere 

showed  good  evidence  of  the  inertia...  subrange  with  -5/3  slope. 

Further  data  collection  should  emphasize  the  low  frequencies,  since 

this  appears  to  be  the  low  noise  region.  Occasional  anomalous  and  low 

signal- to-noise  spectra  were  encountered  above  the  PBL.  Low  signal 

2  2 

spectra  represent  quiescent  layers .  Cj  and  Cu  were 
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calculated  from  the  regression  fit  to  the  inertial  subrange  power 
spectral  density.  The  rate  of  dissipation  of  turbulent  kinetic 

energy,  f,  can  be  calculated  by  using  the  so-called  Corrsin  relation 

2  2 
for  velocity  turbulence.  Cn  can  be  calculated  from  C^.  . 

Considerable  interest  has  developed  in  models  that  relate 

microturbulence  parameters  to  the  mean  gradients.  One  such  model, 

proposed  by  VanZandt  et  al.  <1978;  1981),  has  as  a  key  variable  the 

gradient  Richardson  number,  R^.  In  actively  turbulent  regions,  R; 

2  2 

can  be  related  to  and  as  follows: 

-  1.6  R./  (Pr  -  Rt)  <0/g)  (30/3z)  (8.1) 

2  2 

A  scatter  plot  of  versus  (0/g  30/3z)  Cu  in  the 

free  troposphere  has  a  slope  that  corresponds  to  1.6  Ri/(Pr  * 

R^) .  This  plot  showed  a  range  of  values  from  0.3  to  10  for  high  and 

2  2 
low  values  of  C^.  ,  respectively.  Scatter  plots  of  versus 

2 

showed  high  correlation  of  these  two  parameters,  with  log 
2  2 

Cu  -  1.7  log  Cj,  .  This  implies  that  the  mean  gradient 

structure  of  the  free  troposphere  would  have  a  probability  density 

maximum  (Fairall  and  Markson  1985).  This  could  be  illustrated  as  a 

single  curve  in  N-S  space,  where  N  is  the  Brunt-Vaisala  frequency 

squared,  and  S  is  the  shear  squared  and  R^  -  N/S .  Plots  of  this 

type  were  done  using  aircraft  temperature  data  and  profiler  wind  data. 

A  distribution  was  obtained  roughly  corresponding  to  the  curve 

depicted  by  Fairall  and  Markson  (1985).  Generally,  regions  of  high 

shear  are  at  or  below  the  R.-  0.25  line,  as  expected. 

2  2 

Plots  of  /  Cu  (0/g  30/3z)  versus  a  turbulent  activity 
parameter  (Gregg  1987)  clearly  show  agreement  with  the  suggested 
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activity  levels  and  associated  values .  As  the  ordinate  value 
(corresponding  to  1.6  R^/(Pr  *  Rj.)]  decreases,  turbulent 
activity  increases.  At  high  activity  levels,  the  ordinate  value 
approaches  0.4  (equivalent  to  1.6  when  approximating  P^  as  1 
and  R^  as  0.25). 

2 

Values  of  N  and  S  were  used  to  obtain  values  of  C  and  « 

n 

2 

with  the  model  of  VanZandt  et  al.  Comparisons  of  Cn  profiles 

measured  by  the  various  instruments  showed  good  agreement  between 

scintillometer  and  thermosonde,  and  between  the  VanZandt  et  al.  model 

and  the  aircraft  data.  The  profiler  (converted  optical  values)  had 

2 

some  agreement  with  the  thermosonde.  The  difference  in  Cn 

profiles  between  methods  varied.  An  average  ratio  of  each  profile 

versus  the  thermosonde  profile  gave  the  following  factors  for  flight 

6:  aircraft  11.7  (standard  deviation  23.4),  model  11.5  (stand,  dev. 

18.5),  profiler  4.3  (stand,  dev.  5.3),  and  the  scintillometer  1.6 

(stand,  dev.  0.5).  In  the  lower  levels  of  the  flight  6  ascent,  the 

aircraft  profile  is  thought  to  have  an  error  in  gain  setting.  For 

flight  8,  the  factors  were:  aircraft  4.2  (stand,  dev.  4.5),  model  2. 

(stand,  dev.  2),  profiler  12  (stand,  dev.  34),  and  the  scincillometer 

0.7  (stand,  dev.  0.6).  Similar  differences  between  aircraft  and 

thermosonde  were  noted  by  Brown  and  Good  (1984) .  It  was  later  though 

that  this  difference  might  have  been  due  to  spectral  analysis 

technique.  However,  the  increased  effort  to  improve  spectral  analysi 

in  this  paper  apparently  seemed  to  have  little  effect.  Structure 

constant  values  obtained  by  the  two  spectral  analysis  methods  examine 

in  this  work  (FFT  board  and  the  signal  analyzer)  differed  on  the 

2 

average  by  about  a  factor  of  two  for  :  however  individual 
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values  differed  by  as  much  as  Cvo  orders  of  magnitude. 

2 

To  account  for  the  difference  in  Cn  profiles  between 

instrument  systems,  further  research  might  seek  to  determine  the 

difference  in  values  obtained  by  thermosonde  and  aircraft 

instrumentation,  since  they  actually  measure  some  of  the  same 

2 

variables  (e.g.,  )  in  slightly  different  ways.  Since  the 

scintillometer  tended  to  agree  with  the  thermosonde  and  the  VanZandt 
et  al.  model  tended  to  agree  with  the  aircraft,  no  clear  choice  is 
evident.  The  radar  does  not  have  an  absolute  calibration,  so  there  is 
no  significance  to  its  agreement  with  either  system  (Christopher 
Fairall,  personal  communication,  December  1987). 

Another  noticeable  difference  between  the  aircraft  and 
thermosonde  is  the  measured  vertical  profiles  of  potential 
temperature.  Since  potential  temperature  and  its  gradient  are  input 
parameters  in  the  VanZandt  et  al.  model,  this  difference  could 
obviously  affect  model  output.  Barat  and  Bertin  (1984)  state  the  need 
for  accurate  temperature  profiles  and  suggest  changes  in  Richardson 
number  depend  more  on  temperature  gradient  than  on  shear.  This  study 
only  used  aircraft  data  as  model  input.  Comparing  model  output  from 
thermosonde  and  aircraft  input  could  be  another  interesting  comparison 
for  further  study.  There  is  also  an  obvious  need  for  additional 
simultaneous  data  collection.  Logistical  and  weather  problems 
contributed  to  make  only  two  nights  of  simultaneous  operation  possible 
during  EWAK. 

2 

Numerical  integration  of  C^  gave  estimates  of  rQ  and  9  . 

Values  from  scintillometer,  thermosonde  and  aircraft  did  not  show 
consistent  agreement.  Based  on  the  aircraft  values,  both  the 
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scintillometer  and  thermosonde  differed  by  an  average  of  17%  for  8  , 

and  the  thermosonde  differed  by  an  average  of  35%  for  r  .  With  such 

a  limited  data  set  it  is  not  possible  at  this  time  to  conclude  whether 
2 

profiles  of  can  be  used  to  estimate  these  optical  parameters. 
Acceptable  approximations  for  unmeasured  segments  of  the  atmosphere 
must  also  be  further  studied. 

2 

Despite  good  agreement  on  Cn  ,  the  aircraft  data  and  the 

VanZandt  et  al.  model  output  clearly  disagreed  on  vertical  profiles  of 

t.  This  difference  was  also  noted  by  Fairall  and  Markson  (1985). 

They  suggest  this  effect  is  due  to  the  adjustment  of  the  model 

parameter  L,  turbulent  layer  thickness,  to  produce  the  best  fit  of 

2 

model  output  to  profiler  Cn  data.  Since  few  measurements  of  e 

have  been  done,  there  is  a  definite  need  for  a  larger  data  base. 

Other  areas  of  study  might  include  adjusting  L  to  best  fit  aircraft  e 

profiles.  However,  this  would  require  either  adjusting  the  constants 

(a  and  b)  in  the  model  equations  or  modifying  the  theory  in  order  to 

2 

preserve  the  good  predictions  for  .  Raw  high  speed  aircraft 
data  could  also  be  used  to  study  turbulent  "episodes"  individually 
without  vertical  averaging.  The  fraction  of  the  profile  that  is 
actually  turbulent  could  be  computed  and  the  microturbulence  ratio 
could  be  examined  for  each  active  layer. 

Overall,  EWAK  provided  a  chance  to  compare  methods  of  measuring 

2 

Cn  ,  as  well  as  an  opportunity  to  obtain  microturbulence  data  that 
has  some  relation  to  the  mean  gradient  structure.  Hopefully,  further 
research  will  provide  a  relationship  that  could  be  exploited  in  future 
predictive  models. 
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RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (CSI)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  C5I  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /maintainability  and  compatibility. 


